Showing 1–50 of 55 results for author: Waldmann, I P

The effect of spectroscopic binning on atmospheric retrievals

Authors: Jack J. Davey, Kai Hou Yip, Ahmed F. Al-Refaie, Ingo P. Waldmann

Abstract: With the James Webb Space Telescope (JWST) offering higher resolution data in space-based transmission spectroscopy, understanding the capabilities of our current atmospheric retrieval pipelines is essential. These new data cover wider wavelength ranges and at much higher spectral resolution than previous instruments have been able to offer. Therefore, it is often appealing to bin spectra to fewer… ▽ More With the James Webb Space Telescope (JWST) offering higher resolution data in space-based transmission spectroscopy, understanding the capabilities of our current atmospheric retrieval pipelines is essential. These new data cover wider wavelength ranges and at much higher spectral resolution than previous instruments have been able to offer. Therefore, it is often appealing to bin spectra to fewer points, better constrained in their transit depth, before using them as inputs for atmospheric retrievals. However, little quantitative analysis of the trade-off between spectral resolution and signal-to-noise ratio has been conducted thus far. As such, we produce a simulation replicating the observations of WASP-39b by the NIRSpec PRISM instrument on board JWST and assess the accuracy and consistency of retrievals while varying resolution and the average photometric error. While this probes a specific case we also plot `binning paths' in the resulting sensitivity maps to demonstrate the best attainable atmospheric parameter estimations starting from the position of the real JWST Early Release Science observation. We repeat this analysis on three different simulation setups where each includes an opaque cloud layer at a different height in the atmosphere. We find that a much greater resolution is needed in the case of a high cloud deck since features are already heavily muted by the presence of the clouds. In the other two cases, there are large `safe zones' in the parameter space. If these maps can be generalised, binning paths could inform future observations on how to achieve the most accurate retrieval results. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: Submitted to Monthly Notices of the Royal Astronomical Society, 25 pages, 27 figures

Enhancing 3D Planetary Atmosphere Simulations with a Surrogate Radiative Transfer Model

Authors: Tara P. A. Tahseen, João M. Mendonça, Kai Hou Yip, Ingo P. Waldmann

Abstract: This work introduces an approach to enhancing the computational efficiency of 3D atmospheric simulations by integrating a machine-learned surrogate model into the OASIS global circulation model (GCM). Traditional GCMs, which are based on repeatedly numerically integrating physical equations governing atmospheric processes across a series of time-steps, are time-intensive, leading to compromises in… ▽ More This work introduces an approach to enhancing the computational efficiency of 3D atmospheric simulations by integrating a machine-learned surrogate model into the OASIS global circulation model (GCM). Traditional GCMs, which are based on repeatedly numerically integrating physical equations governing atmospheric processes across a series of time-steps, are time-intensive, leading to compromises in spatial and temporal resolution of simulations. This research improves upon this limitation, enabling higher resolution simulations within practical timeframes. Speeding up 3D simulations holds significant implications in multiple domains. Firstly, it facilitates the integration of 3D models into exoplanet inference pipelines, allowing for robust characterisation of exoplanets from a previously unseen wealth of data anticipated from JWST and post-JWST instruments. Secondly, acceleration of 3D models will enable higher resolution atmospheric simulations of Earth and Solar System planets, enabling more detailed insights into their atmospheric physics and chemistry. Our method replaces the radiative transfer module in OASIS with a recurrent neural network-based model trained on simulation inputs and outputs. Radiative transfer is typically one of the slowest components of a GCM, thus providing the largest scope for overall model speed-up. The surrogate model was trained and tested on the specific test case of the Venusian atmosphere, to benchmark the utility of this approach in the case of non-terrestrial atmospheres. This approach yields promising results, with the surrogate-integrated GCM demonstrating above 99.0% accuracy and 101 factor GPU speed-up of the entire simulation compared to using the matched original GCM under Venus-like conditions. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 17 pages, 11 figures

Is the atmosphere of the ultra-hot Jupiter WASP-121b variable?

Authors: Quentin Changeat, Jack W. Skinner, James Y-K. Cho, Joonas Nättilä, Ingo P. Waldmann, Ahmed F. Al-Refaie, Achrène Dyrek, Billy Edwards, Thomas Mikal-Evans, Max Joshua, Giuseppe Morello, Nour Skaf, Angelos Tsiaras, Olivia Venot, Kai Hou Yip

Abstract: We present a comprehensive analysis of the Hubble Space Telescope observations of the atmosphere of WASP-121 b, a 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 cl… ▽ More We present a comprehensive analysis of the Hubble Space Telescope observations of the atmosphere of WASP-121 b, a 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 ~0.1 bar on the dayside, solar to subsolar metallicity Z (i.e., -0.77 < log(Z) < 0.05), and super-solar 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 atmosphere at high-resolution, we show that the identified signals are consistent with quasi-periodic 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. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: Accepted for publication in ApJS. 43 pages, 31 figures, 2 animations (available online at the journal)

Retrieval study of cool, directly imaged exoplanet 51 Eri b

Authors: Niall Whiteford, Alistair Glasse, Katy L. Chubb, Daniel Kitzmann, Shrishmoy Ray, Mark W. Phillips, Beth A. Biller, Paul I. Palmer, Ken Rice, Ingo P. Waldmann, Quentin Changeat, Nour Skaf, Jason Wang, Billy Edwards, Ahmed Al-Refaie

Abstract: Retrieval methods are a powerful analysis technique for modelling exoplanetary atmospheres by estimating the bulk physical and chemical properties that combine in a forward model to best-fit an observed spectrum, and they are increasingly being applied to observations of directly-imaged exoplanets. We have adapted TauREx3, the Bayesian retrieval suite, for the analysis of near-infrared spectrophot… ▽ More Retrieval methods are a powerful analysis technique for modelling exoplanetary atmospheres by estimating the bulk physical and chemical properties that combine in a forward model to best-fit an observed spectrum, and they are increasingly being applied to observations of directly-imaged exoplanets. We have adapted TauREx3, the Bayesian retrieval suite, for the analysis of near-infrared spectrophotometry from directly-imaged gas giant exoplanets and brown dwarfs. We demonstrate TauREx3's applicability to sub-stellar atmospheres by presenting results for brown dwarf benchmark GJ 570D which are consistent with previous retrieval studies, whilst also exhibiting systematic biases associated with the presence of alkali lines. We also present results for the cool exoplanet 51 Eri b, the first application of a free chemistry retrieval analysis to this object, using spectroscopic observations from GPI and SPHERE. While our retrieval analysis is able to explain spectroscopic and photometric observations without employing cloud extinction, we conclude this may be a result of employing a flexible temperature-pressure profile which is able to mimic the presence of clouds. We present Bayesian evidence for an ammonia detection with a 2.7$σ$ confidence, the first indication of ammonia in an exoplanetary atmosphere. This is consistent with this molecule being present in brown dwarfs of a similar spectral type. We demonstrate the chemical similarities between 51 Eri b and GJ 570D in relation to their retrieved molecular abundances. Finally, we show that overall retrieval conclusions for 51 Eri b can vary when employing different spectral data and modelling components, such as temperature-pressure and cloud structures. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Identification of carbon dioxide in an exoplanet atmosphere

Authors: The JWST Transiting Exoplanet Community Early Release Science Team, Eva-Maria Ahrer, Lili Alderson, Natalie M. Batalha, Natasha E. Batalha, Jacob L. Bean, Thomas G. Beatty, Taylor J. Bell, Björn Benneke, Zachory K. Berta-Thompson, Aarynn L. Carter, Ian J. M. Crossfield, Néstor Espinoza, Adina D. Feinstein, Jonathan J. Fortney, Neale P. Gibson, Jayesh M. Goyal, Eliza M. -R. Kempton, James Kirk, Laura Kreidberg, Mercedes López-Morales, Michael R. Line, Joshua D. Lothringer, Sarah E. Moran, Sagnick Mukherjee , et al. (107 additional authors not shown)

Abstract: Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres… ▽ More Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2 but have not yielded definitive detections due to the lack of unambiguous spectroscopic identification. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science Program (ERS). The data used in this study span 3.0 to 5.5 μm in wavelength and show a prominent CO2 absorption feature at 4.3 μm (26σ significance). The overall spectrum is well matched by one-dimensional, 10x solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide, and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 μm that is not reproduced by these models. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

ESA-Ariel Data Challenge NeurIPS 2022: Inferring Physical Properties of Exoplanets From Next-Generation Telescopes

Authors: Kai Hou Yip, Ingo P. Waldmann, Quentin Changeat, Mario Morvan, Ahmed F. Al-Refaie, Billy Edwards, Nikolaos Nikolaou, Angelos Tsiaras, Catarina Alves de Oliveira, Pierre-Olivier Lagage, Clare Jenner, James Y-K. Cho, Jeyan Thiyagalingam, Giovanna Tinetti

Abstract: The study of extra-solar planets, or simply, exoplanets, planets outside our own Solar System, is fundamentally a grand quest to understand our place in the Universe. Discoveries in the last two decades have re-defined our understanding of planets, and helped us comprehend the uniqueness of our very own Earth. In recent years the focus has shifted from planet detection to planet characterisation,… ▽ More The study of extra-solar planets, or simply, exoplanets, planets outside our own Solar System, is fundamentally a grand quest to understand our place in the Universe. Discoveries in the last two decades have re-defined our understanding of planets, and helped us comprehend the uniqueness of our very own Earth. In recent years the focus has shifted from planet detection to planet characterisation, where key planetary properties are inferred from telescope observations using Monte Carlo-based methods. However, the efficiency of sampling-based methodologies is put under strain by the high-resolution observational data from next generation telescopes, such as the James Webb Space Telescope and the Ariel Space Mission. We are delighted to announce the acceptance of the Ariel ML Data Challenge 2022 as part of the NeurIPS competition track. The goal of this challenge is to identify a reliable and scalable method to perform planetary characterisation. Depending on the chosen track, participants are tasked to provide either quartile estimates or the approximate distribution of key planetary properties. To this end, a synthetic spectroscopic dataset has been generated from the official simulators for the ESA Ariel Space Mission. The aims of the competition are three-fold. 1) To offer a challenging application for comparing and advancing conditional density estimation methods. 2) To provide a valuable contribution towards reliable and efficient analysis of spectroscopic data, enabling astronomers to build a better picture of planetary demographics, and 3) To promote the interaction between ML and exoplanetary science. The competition is open from 15th June and will run until early October, participants of all skill levels are more than welcomed! △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: 13 pages, Accepted in the Conference on Neural Information Processing Systems 2022 (NeurIPS 2022), Competition hosted on https://www.ariel-datachallenge.space/

Five key exoplanet questions answered via the analysis of 25 hot Jupiter atmospheres in eclipse

Authors: Quentin Changeat, Billy Edwards, Ahmed F. Al-Refaie, Angelos Tsiaras, Jack W. Skinner, James Y-K Cho, Kai H. Yip, Lara Anisman, Masahiro Ikoma, Michelle F. Bieger, Olivia Venot, Sho Shibata, Ingo P. Waldmann, Giovanna Tinetti

Abstract: Population studies of exoplanets are key to unlocking their statistical properties. So far the inferred properties have been mostly limited to planetary, orbital and stellar parameters extracted from, e.g., Kepler, radial velocity, and GAIA data. More recently an increasing number of exoplanet atmospheres have been observed in detail from space and the ground. Generally, however, these atmospheric… ▽ More Population studies of exoplanets are key to unlocking their statistical properties. So far the inferred properties have been mostly limited to planetary, orbital and stellar parameters extracted from, e.g., Kepler, radial velocity, and GAIA data. More recently an increasing number of exoplanet atmospheres have been observed in detail from space and the ground. Generally, however, these atmospheric studies have focused on individual planets, with the exception of a couple of works which have detected the presence of water vapor and clouds in populations of gaseous planets via transmission spectroscopy. Here, using a suite of retrieval tools, we analyse spectroscopic and photometric data of 25 hot Jupiters, obtained with the Hubble and Spitzer Space Telescopes via the eclipse technique. By applying the tools uniformly across the entire set of 25 planets, we extract robust trends in the thermal structure and chemical properties of hot Jupiters not obtained in past studies. With the recent launch of JWST and the upcoming missions Twinkle, and Ariel, population based studies of exoplanet atmospheres, such as the one presented here, will be a key approach to understanding planet characteristics, formation, and evolution in our galaxy. △ Less

Submitted 13 May, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: 66 pages, 23 figures, 7 tables. Published in The Astrophysical Journal Supplement Series

Journal ref: ApJS 260 3 (2022)

A retrieval challenge exercise for the Ariel mission

Authors: Joanna K. Barstow, Quentin Changeat, Katy L. Chubb, Patricio E. Cubillos, Billy Edwards, Ryan J. MacDonald, Michiel Min, Ingo P. Waldmann

Abstract: The Ariel mission, due to launch in 2029, will obtain spectroscopic information for 1000 exoplanets, providing an unprecedented opportunity for comparative exoplanetology. Retrieval codes - parameteric atmospheric models coupled with an inversion algorithm - represent the tool of choice for interpreting Ariel data. Ensuring that reliable and consistent results can be produced by these tools is a c… ▽ More The Ariel mission, due to launch in 2029, will obtain spectroscopic information for 1000 exoplanets, providing an unprecedented opportunity for comparative exoplanetology. Retrieval codes - parameteric atmospheric models coupled with an inversion algorithm - represent the tool of choice for interpreting Ariel data. Ensuring that reliable and consistent results can be produced by these tools is a critical preparatory step for the mission. Here, we present the results of a retrieval challenge. We use five different exoplanet retrieval codes to analyse the same synthetic datasets, and test a) the ability of each to recover the correct input solution and b) the consistency of the results. We find that generally there is very good agreement between the five codes, and in the majority of cases the correct solutions are recovered. This demonstrates the reproducibility of retrievals for transit spectra of exoplanets, even when codes are not previously benchmarked against each other. △ Less

Submitted 15 March, 2022; v1 submitted 1 March, 2022; originally announced March 2022.

Comments: 28 pages, 14 figures. Accepted in Experimental Astronomy (2022)

A comparison of chemical models of exoplanet atmospheres enabled by TauREx 3.1

Authors: Ahmed F. Al-Refaie, Quentin. Changeat, Olivia Venot, Ingo P. Waldmann, Giovanna Tinetti

Abstract: Thermochemical equilibrium is one of the most commonly used assumptions in current exoplanet retrievals. As the James Webb Space Telescope (JWST) and Ariel launch draw near, assessing the underlying biases and assumptions made when applying self-consistent chemistry into spectral retrievals is crucial. Here we use the flexibility of TauREx 3.1 to cross-compare three state of the art chemical equil… ▽ More Thermochemical equilibrium is one of the most commonly used assumptions in current exoplanet retrievals. As the James Webb Space Telescope (JWST) and Ariel launch draw near, assessing the underlying biases and assumptions made when applying self-consistent chemistry into spectral retrievals is crucial. Here we use the flexibility of TauREx 3.1 to cross-compare three state of the art chemical equilibrium codes: ACE, FastChem and GGchem. We simulate JWST spectra for ACE, FastChem, GGchem and GGchem+condensation containing only C, H, O, N elements and spectra for FastChem, GGchem and GGchem+condensation with a more extensive range of elements giving seven simulated JWST spectra in total and then cross-retrieve giving a total of 49 retrievals. Our analysis demonstrates that like-for-like, all chemical codes retrieve the correct parameters to <1% of the truth. However, retrievals, where the contained elements do not match the truth, parameters such as metallicity deviate by 20% while maintaining extremely low uncertainties <1% giving false confidence. This point is of major importance for future analyses on JWST and Ariel, highlighting that self-consistent chemical scheme that do not employ the proper assumptions (missing species, fixed elemental ratios, condensation) are at risk of confidently biasing interpretations. Free chemistry retrievals employing parametric descriptions of the chemical profiles can provide alternative unbiased explorations. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: 31 pages, 15 figures

Ariel: Enabling planetary science across light-years

Authors: Giovanna Tinetti, Paul Eccleston, Carole Haswell, Pierre-Olivier Lagage, Jérémy Leconte, Theresa Lüftinger, Giusi Micela, Michel Min, Göran Pilbratt, Ludovic Puig, Mark Swain, Leonardo Testi, Diego Turrini, Bart Vandenbussche, Maria Rosa Zapatero Osorio, Anna Aret, Jean-Philippe Beaulieu, Lars Buchhave, Martin Ferus, Matt Griffin, Manuel Guedel, Paul Hartogh, Pedro Machado, Giuseppe Malaguti, Enric Pallé , et al. (293 additional authors not shown)

Abstract: Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths.… ▽ More Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

Report number: ESA/SCI(2020)1

An exploration of model degeneracies with a unified phase curve retrieval analysis: The light and dark sides of WASP-43 b

Authors: Quentin Changeat, Ahmed F. Al-Refaie, Billy Edwards, Ingo P. Waldmann, Giovanna Tinetti

Abstract: The analysis of exoplanetary atmospheres often relies upon the observation of transit or eclipse events. While very powerful, these snapshots provide mainly 1-dimensional information on the planet structure and do not easily allow precise latitude-longitude characterisations. The phase curve technique, which consists of measuring the planet emission throughout its entire orbit, can break this limi… ▽ More The analysis of exoplanetary atmospheres often relies upon the observation of transit or eclipse events. While very powerful, these snapshots provide mainly 1-dimensional information on the planet structure and do not easily allow precise latitude-longitude characterisations. The phase curve technique, which consists of measuring the planet emission throughout its entire orbit, can break this limitation and provide useful 2-dimensional thermal and chemical constraints on the atmosphere. As of today however, computing performances have limited our ability to perform unified retrieval studies on the full set of observed spectra from phase curve observations at the same time. Here, we present a new phase curve model that enables fast, unified retrieval capabilities. We apply our technique to the combined phase curve data from the Hubble and Spitzer space telescopes of the hot-Jupiter WASP-43 b. We tested different scenarios and discussed the dependence of our solution to different assumptions in the model. Our more comprehensive approach suggests that multiple interpretation of this dataset are possible but our more complex model is consistent with the presence of thermal inversions and a metal rich atmosphere, contrasting with previous data analyses, although this likely depends on the Spitzer data reduction. The detailed constraints extracted here demonstrate the importance of developing and understanding advanced phase curve techniques, which we believe will unlock access to a richer picture of exoplanet atmospheres. △ Less

Submitted 13 May, 2021; v1 submitted 26 March, 2021; originally announced March 2021.

Comments: Accepted in ApJ, 57 pages, 38 figures

Peeking inside the Black Box: Interpreting Deep Learning Models for Exoplanet Atmospheric Retrievals

Authors: Kai Hou Yip, Quentin Changeat, Nikolaos Nikolaou, Mario Morvan, Billy Edwards, Ingo P. Waldmann, Giovanna Tinetti

Abstract: Deep learning algorithms are growing in popularity in the field of exoplanetary science due to their ability to model highly non-linear relations and solve interesting problems in a data-driven manner. Several works have attempted to perform fast retrievals of atmospheric parameters with the use of machine learning algorithms like deep neural networks (DNNs). Yet, despite their high predictive pow… ▽ More Deep learning algorithms are growing in popularity in the field of exoplanetary science due to their ability to model highly non-linear relations and solve interesting problems in a data-driven manner. Several works have attempted to perform fast retrievals of atmospheric parameters with the use of machine learning algorithms like deep neural networks (DNNs). Yet, despite their high predictive power, DNNs are also infamous for being 'black boxes'. It is their apparent lack of explainability that makes the astrophysics community reluctant to adopt them. What are their predictions based on? How confident should we be in them? When are they wrong and how wrong can they be? In this work, we present a number of general evaluation methodologies that can be applied to any trained model and answer questions like these. In particular, we train three different popular DNN architectures to retrieve atmospheric parameters from exoplanet spectra and show that all three achieve good predictive performance. We then present an extensive analysis of the predictions of DNNs, which can inform us - among other things - of the credibility limits for atmospheric parameters for a given instrument and model. Finally, we perform a perturbation-based sensitivity analysis to identify to which features of the spectrum the outcome of the retrieval is most sensitive. We conclude that for different molecules, the wavelength ranges to which the DNN's predictions are most sensitive, indeed coincide with their characteristic absorption regions. The methodologies presented in this work help to improve the evaluation of DNNs and to grant interpretability to their predictions. △ Less

Submitted 23 July, 2021; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: 32 pages, 22 figures. Submitted to AJ. Accepted July 20, 2021

PyLightcurve-torch: a transit modelling package for deep learning applications in PyTorch

Authors: Mario Morvan, Angelos Tsiaras, Nikolaos Nikolaou, Ingo P. Waldmann

Abstract: We present a new open source python package, based on PyLightcurve and PyTorch, tailored for efficient computation and automatic differentiation of exoplanetary transits. The classes and functions implemented are fully vectorised, natively GPU-compatible and differentiable with respect to the stellar and planetary parameters. This makes PyLightcurve-torch suitable for traditional forward computati… ▽ More We present a new open source python package, based on PyLightcurve and PyTorch, tailored for efficient computation and automatic differentiation of exoplanetary transits. The classes and functions implemented are fully vectorised, natively GPU-compatible and differentiable with respect to the stellar and planetary parameters. This makes PyLightcurve-torch suitable for traditional forward computation of transits, but also extends the range of possible applications with inference and optimisation algorithms requiring access to the gradients of the physical model. This endeavour is aimed at fostering the use of deep learning in exoplanets research, motivated by an ever increasing amount of stellar light curves data and various incentives for the improvement of detection and characterisation techniques. △ Less

Submitted 28 December, 2020; v1 submitted 3 November, 2020; originally announced November 2020.

Comments: 7 pages, 3 figures; submission status updated, fig 2 caption added

Lessons Learned from the 1st ARIEL Machine Learning Challenge: Correcting Transiting Exoplanet Light Curves for Stellar Spots

Authors: Nikolaos Nikolaou, Ingo P. Waldmann, Angelos Tsiaras, Mario Morvan, Billy Edwards, Kai Hou Yip, Giovanna Tinetti, Subhajit Sarkar, James M. Dawson, Vadim Borisov, Gjergji Kasneci, Matej Petkovic, Tomaz Stepisnik, Tarek Al-Ubaidi, Rachel Louise Bailey, Michael Granitzer, Sahib Julka, Roman Kern, Patrick Ofner, Stefan Wagner, Lukas Heppe, Mirko Bunse, Katharina Morik

Abstract: The last decade has witnessed a rapid growth of the field of exoplanet discovery and characterisation. However, several big challenges remain, many of which could be addressed using machine learning methodology. For instance, the most prolific method for detecting exoplanets and inferring several of their characteristics, transit photometry, is very sensitive to the presence of stellar spots. The… ▽ More The last decade has witnessed a rapid growth of the field of exoplanet discovery and characterisation. However, several big challenges remain, many of which could be addressed using machine learning methodology. For instance, the most prolific method for detecting exoplanets and inferring several of their characteristics, transit photometry, is very sensitive to the presence of stellar spots. The current practice in the literature is to identify the effects of spots visually and correct for them manually or discard the affected data. This paper explores a first step towards fully automating the efficient and precise derivation of transit depths from transit light curves in the presence of stellar spots. The methods and results we present were obtained in the context of the 1st Machine Learning Challenge organized for the European Space Agency's upcoming Ariel mission. We first present the problem, the simulated Ariel-like data and outline the Challenge while identifying best practices for organizing similar challenges in the future. Finally, we present the solutions obtained by the top-5 winning teams, provide their code and discuss their implications. Successful solutions either construct highly non-linear (w.r.t. the raw data) models with minimal preprocessing -deep neural networks and ensemble methods- or amount to obtaining meaningful statistics from the light curves, constructing linear models on which yields comparably good predictive performance. △ Less

Submitted 29 October, 2020; originally announced October 2020.

Comments: 20 pages, 7 figures, 2 tables, Submitted to The Astrophysics Journal (ApJ)

KELT-11 b: Abundances of water and constraints on carbon-bearing molecules from the Hubble transmission spectrum

Authors: Quentin Changeat, Billy Edwards, Ahmed F. Al-Refaie, Mario Morvan, Angelos Tsiaras, Ingo P. Waldmann, Giovanna Tinetti

Abstract: In the past decade, the analysis of exoplanet atmospheric spectra has revealed the presence of water vapour in almost all the planets observed, with the exception of a fraction of overcast planets. Indeed, water vapour presents a large absorption signature in the wavelength coverage of the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3), which is the main space-based observatory for atmo… ▽ More In the past decade, the analysis of exoplanet atmospheric spectra has revealed the presence of water vapour in almost all the planets observed, with the exception of a fraction of overcast planets. Indeed, water vapour presents a large absorption signature in the wavelength coverage of the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3), which is the main space-based observatory for atmospheric studies of exoplanets, making its detection very robust. However, while carbon-bearing species such as methane, carbon monoxide and carbon dioxide are also predicted from current chemical models, their direct detection and abundance characterisation has remained a challenge. Here we analyse the transmission spectrum of the puffy, clear hot-Jupiter KELT-11 b from the HST WFC3 camera. We find that the spectrum is consistent with the presence of water vapor and an additional absorption at longer wavelengths than 1.5um, which could well be explained by a mix of carbon bearing molecules. CO2, when included is systematically detected. One of the main difficulties to constrain the abundance of those molecules is their weak signatures across the HST WFC3 wavelength coverage, particularly when compared to those of water. Through a comprehensive retrieval analysis, we attempt to explain the main degeneracies present in this dataset and explore some of the recurrent challenges that are occurring in retrieval studies (e.g: the impact of model selection, the use of free vs self-consistent chemistry and the combination of instrument observations). Our results make this planet an exceptional example of chemical laboratory where to test current physical and chemical models of hot-Jupiters' atmospheres. △ Less

Submitted 3 October, 2020; originally announced October 2020.

Comments: 24 pages, 14 figures, Accepted in AJ

On The Compatibility of Ground-based and Space-based Data: WASP-96 b, An Example

Authors: Kai Hou Yip, Quentin Changeat, Billy Edwards, Mario Morvan, Katy L. Chubb, Angelos Tsiaras, Ingo P. Waldmann, Giovanna Tinetti

Abstract: The study of exoplanetary atmospheres relies on detecting minute changes in the transit depth at different wavelengths. To date, a number of ground and space based instruments have been used to obtain transmission spectra of exoplanets in different spectral band. One common practice is to combine observations from different instruments in order to achieve a broader wavelength coverage. We present… ▽ More The study of exoplanetary atmospheres relies on detecting minute changes in the transit depth at different wavelengths. To date, a number of ground and space based instruments have been used to obtain transmission spectra of exoplanets in different spectral band. One common practice is to combine observations from different instruments in order to achieve a broader wavelength coverage. We present here two inconsistent observations of WASP-96 b, one by Hubble Space Telescope (HST) and the other by the Very Large Telescope (VLT). We present two key findings in our investigation: 1.) a strong water signature is detected via the HST WFC3 observations. 2.) A notable offset in transit depth (> 1100 ppm) can be seen when the ground-based and space-based observations are combined together. The discrepancy raises the question of whether observations from different instruments could indeed be combined together. We attempt to align the observations by including an additional parameter in our retrieval studies but are unable to definitively ascertain that the aligned observations are indeed compatible. The case of WASP-96 b signals that compatibility of instruments should not be assumed. While wavelength overlaps between instruments can help, it should be noted that combining datasets remains a risky business. The difficulty in combining observations also strengthens the need for next generation instruments which will possess broader spectral coverage. △ Less

Submitted 19 October, 2020; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: 18 pages, 13 figures. Submitted to ApJ, currently under review

Alfnoor: A Retrieval Simulation of the Ariel Target List

Authors: Quentin Changeat, Ahmed Al-Refaie, Lorenzo V. Mugnai, Billy Edwards, Ingo P. Waldmann, Enzo Pascale, Giovanna Tinetti

Abstract: In this work, we present Alfnoor, a dedicated tool optimised for population studies of exoplanet atmospheres. Alfnoor combines the latest version of the retrieval algorithm TauREx 3, with the instrument noise simulator ArielRad and enables the simultaneous retrieval analysis of a large sample of exo-atmospheres. We applied this tool to the Ariel list of planetary candidates and focus on hydrogen d… ▽ More In this work, we present Alfnoor, a dedicated tool optimised for population studies of exoplanet atmospheres. Alfnoor combines the latest version of the retrieval algorithm TauREx 3, with the instrument noise simulator ArielRad and enables the simultaneous retrieval analysis of a large sample of exo-atmospheres. We applied this tool to the Ariel list of planetary candidates and focus on hydrogen dominated, cloudy atmospheres observed in transit with the Tier-2 mode (medium Ariel resolution). As a first experiment, we randomised the abundances - ranging from 10$^{-7}$ to 10$^{-2}$ - of the trace gases, which include H$_2$O, CH$_4$, CO, CO$_2$ and NH$_3$. This exercise allowed to estimate the detection limits for Ariel Tier-2 and Tier-3 modes when clouds are present. In a second experiment, we imposed an arbitrary trend between a chemical species and the effective temperature of the planet. A last experiment was run requiring molecular abundances being dictated by equilibrium chemistry at a certain temperature. Our results demonstrate the ability of Ariel Tier-2 and Tier-3 surveys to reveal trends between the chemistry and associated planetary parameters. Future work will focus on eclipse data, on atmospheres heavier than hydrogen and will be applied also to other observatories. △ Less

Submitted 9 October, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: 33 pages, 24 figures, Accepted in AJ

Disentangling Atmospheric Compositions of K2-18 b with Next Generation Facilities

Authors: Quentin Changeat, Billy Edwards, Ahmed F. Al-Refaie, Angelos Tsiaras, Ingo P. Waldmann, Giovanna Tinetti

Abstract: Recent analysis of the planet K2-18b has shown the presence of water vapour in its atmosphere. While the H2O detection is significant, the Hubble Space Telescope (HST) WFC3 spectrum suggests three possible solutions of very different nature which can equally match the data. The three solutions are a primary cloudy atmosphere with traces of water vapour (cloudy sub-Neptune), a secondary atmosphere… ▽ More Recent analysis of the planet K2-18b has shown the presence of water vapour in its atmosphere. While the H2O detection is significant, the Hubble Space Telescope (HST) WFC3 spectrum suggests three possible solutions of very different nature which can equally match the data. The three solutions are a primary cloudy atmosphere with traces of water vapour (cloudy sub-Neptune), a secondary atmosphere with a substantial amount (up to 50% Volume Mixing Ratio) of H2O (icy/water world) and/or an undetectable gas such as N2 (super-Earth). Additionally, the atmospheric pressure and the possible presence of a liquid/solid surface cannot be investigated with currently available observations. In this paper we used the best fit parameters from Tsiaras et al. (2019) to build James Webb Space Telescope (JWST) and Ariel simulations of the three scenarios. We have investigated 18 retrieval cases, which encompass the three scenarios and different observational strategies with the two observatories. Retrieval results show that twenty combined transits should be enough for the Ariel mission to disentangle the three scenarios, while JWST would require only two transits if combining NIRISS and NIRSpec data. This makes K2-18b an ideal target for atmospheric follow-ups by both facilities and highlights the capabilities of the next generation of space-based infrared observatories to provide a complete picture of low mass planets. △ Less

Submitted 21 August, 2021; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: Accepted in Experimental Astronomy, 18 pages, 13 figures

A comparison of exoplanet spectroscopic retrieval tools

Authors: Joanna K. Barstow, Quentin Changeat, Ryan Garland, Michael R. Line, Marco Rocchetto, Ingo P. Waldmann

Abstract: Over the last several years, spectroscopic observations of transiting exoplanets have begun to uncover information about their atmospheres, including atmospheric composition and indications of the presence of clouds and hazes. Spectral retrieval is the leading technique for interpretation of transmission spectra and is employed by several teams using a variety of forward models and parameter estim… ▽ More Over the last several years, spectroscopic observations of transiting exoplanets have begun to uncover information about their atmospheres, including atmospheric composition and indications of the presence of clouds and hazes. Spectral retrieval is the leading technique for interpretation of transmission spectra and is employed by several teams using a variety of forward models and parameter estimation algorithms. However, different model suites have mostly been used in isolation and so it is unknown whether the results from each are comparable. As we approach the launch of the James Webb Space Telescope we anticipate advances in wavelength coverage, precision, and resolution of transit spectroscopic data, so it is important that the tools that will be used to interpret these information rich spectra are validated. To this end, we present an inter-model comparison of three retrieval suites: TauREx, NEMESIS and CHIMERA. We demonstrate that the forward model spectra are in good agreement (residual deviations on the order of 20 - 40 ppm), and discuss the results of cross retrievals between the three tools. Generally, the constraints from the cross retrievals are consistent with each other and with input values to within 1 sigma However, for high precision scenarios with error envelopes of order 30 ppm, subtle differences in the simulated spectra result in discrepancies between the different retrieval suites, and inaccuracies in retrieved values of several sigma. This can be considered analogous to substantial systematic/astrophysical noise in a real observation, or errors/omissions in a forward model such as molecular linelist incompleteness or missing absorbers. △ Less

Submitted 9 March, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: 25 pages, 21 figures. Accepted in MNRAS

Global Chemistry and Thermal Structure Models for the Hot Jupiter WASP-43b and Predictions for JWST

Authors: Olivia Venot, Vivien Parmentier, Jasmina Blecic, Patricio E. Cubillos, Ingo P. Waldmann, Quentin Changeat, Julianne I. Moses, Pascal Tremblin, Nicolas Crouzet, Peter Gao, Diana Powell, Pierre-Olivier Lagage, Ian Dobbs-Dixon, Maria E. Steinrueck, Laura Kreidberg, Natalie Batalha, Jacob L. Bean, Kevin B. Stevenson, Sarah Casewell, Ludmila Carone

Abstract: The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, in… ▽ More The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $μ$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $μ$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition. △ Less

Submitted 14 January, 2020; originally announced January 2020.

Comments: 26 pages, 16 figures, accepted in The Astrophysical Journal

Detrending Exoplanetary Transit Light Curves with Long Short-Term Memory Networks

Authors: Mario Morvan, Nikolaos Nikolaou, Angelos Tsiaras, Ingo P. Waldmann

Abstract: The precise derivation of transit depths from transit light curves is a key component for measuring exoplanet transit spectra, and henceforth for the study of exoplanet atmospheres. However, it is still deeply affected by various kinds of systematic errors and noise. In this paper we propose a new detrending method by reconstructing the stellar flux baseline during transit time. We train a probabi… ▽ More The precise derivation of transit depths from transit light curves is a key component for measuring exoplanet transit spectra, and henceforth for the study of exoplanet atmospheres. However, it is still deeply affected by various kinds of systematic errors and noise. In this paper we propose a new detrending method by reconstructing the stellar flux baseline during transit time. We train a probabilistic Long Short-Term Memory (LSTM) network to predict the next data point of the light curve during the out-of-transit, and use this model to reconstruct a transit-free light curve - i.e. including only the systematics - during the in-transit. By making no assumption about the instrument, and using only the transit ephemeris, this provides a general way to correct the systematics and perform a subsequent transit fit. The name of the proposed model is TLCD-LSTM, standing for Transit Light Curve Detrending LSTM. Here we present the first results on data from six transit observations of HD 189733b with the IRAC camera on board the Spitzer Space Telescope, and discuss some of its possible further applications. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: 12 pages, 10 figures, 4 tables, accepted for publication in The Astronomical Journal

TauREx III: A fast, dynamic and extendable framework for retrievals

Authors: Ahmed F. Al-Refaie, Quentin Changeat, Ingo P. Waldmann, Giovanna Tinetti

Abstract: TauREx 3 is the next generation of the TauREx exoplanet atmospheric retrieval framework for Windows, Mac, and Linux. It is a complete rewrite with a full Python stack that makes it easy-to-use, high-performance, dynamic, and flexible. The new main TauREx program is built with modularity in mind, allowing the user to augment its functionalities with custom code and efficiently perform retrievals on… ▽ More TauREx 3 is the next generation of the TauREx exoplanet atmospheric retrieval framework for Windows, Mac, and Linux. It is a complete rewrite with a full Python stack that makes it easy-to-use, high-performance, dynamic, and flexible. The new main TauREx program is built with modularity in mind, allowing the user to augment its functionalities with custom code and efficiently perform retrievals on custom parameters. We achieve this result by dynamic determination of fitting parameters, whereby TauREx 3 can detect new parameters for retrieval from user code through a simple interface. TauREx 3 can act as a library with a simple 'import taurex' command, providing a rich set of classes and functions related to atmospheric modelling. A 10x speedup in forward model computations is achieved as compared to the previous version with a sixfold reduction in retrieval times while maintaining robust results. TauREx 3 is intended as a standalone, all-in-one package for retrievals while the TauREx 3 Python library can build or augment a user's custom data pipeline easily. △ Less

Submitted 18 August, 2021; v1 submitted 16 December, 2019; originally announced December 2019.

Journal ref: ApJ 917 (2021) 1-37

Identifiable Acetylene Features Predicted for Young Earth-like Exoplanets with Reducing Atmospheres undergoing Heavy Bombardment

Authors: P. B. Rimmer, M. Ferus, I. P. Waldmann, A Knížek, D. Kalvaitis, O. Ivanek, P. Kubelík, S. N. Yurchenko, T. Burian, J. Dostál, L. Juha, R. Dudžàk, M. Krůs, J. Tennyson, S. Civiš, A. T. Archibald, A. Granville-Willett

Abstract: The chemical environments of young planets are assumed to be largely influenced by impacts of bodies lingering on unstable trajectories after the dissolution of the protoplanetary disk. We explore the chemical consequences of impacts within the context of reducing planetary atmospheres dominated by carbon monoxide, methane and molecular nitrogen. A terawatt high-power laser was selected in order t… ▽ More The chemical environments of young planets are assumed to be largely influenced by impacts of bodies lingering on unstable trajectories after the dissolution of the protoplanetary disk. We explore the chemical consequences of impacts within the context of reducing planetary atmospheres dominated by carbon monoxide, methane and molecular nitrogen. A terawatt high-power laser was selected in order to simulate the airglow plasma and blast wave surrounding the impactor. The chemical results of these experiments are then applied to a theoretical atmospheric model. The impact simulation results in substantial volume mixing ratios within the reactor of 5% hydrogen cyanide (HCN), 8% acetylene (C$_2$H$_2$), 5% cyanoacetylene (HC$_3$N) and 1% ammonia (NH$_3$). These yields are combined with estimated impact rates for the Early Earth to predict surface boundary conditions for an atmospheric model. We show that impacts might have served as sources of energy that would have led to steady-state surface quantities of 0.4% C$_2$H$_2$, 400 ppm HCN and 40 ppm NH$_3$. We provide simulated transit spectra for an Earth-like exoplanet with this reducing atmosphere during and shortly after eras of intense impacts. We predict that acetylene is as observable as other molecular features on exoplanets with reducing atmospheres that have recently gone through their own `Heavy Bombardments', with prominent features at 3.05 $μ$m and 10.5 $μ$m. △ Less

Submitted 5 November, 2019; originally announced November 2019.

Comments: 26 pages, 12 figures, 2 tables, accepted for publication by The Astrophysical Journal

Water vapour in the atmosphere of the habitable-zone eight Earth-mass planet K2-18 b

Authors: Angelos Tsiaras, Ingo P. Waldmann, Giovanna Tinetti, Jonathan Tennyson, Sergey N. Yurchenko

Abstract: In the past decade, observations from space and ground have found H$_2$O to be the most abundant molecular species, after hydrogen, in the atmospheres of hot, gaseous, extrasolar planets. Being the main molecular carrier of oxygen, H$_2$O is a tracer of the origin and the evolution mechanisms of planets. For temperate, terrestrial planets, the presence of H$_2$O is of great significance as an indi… ▽ More In the past decade, observations from space and ground have found H$_2$O to be the most abundant molecular species, after hydrogen, in the atmospheres of hot, gaseous, extrasolar planets. Being the main molecular carrier of oxygen, H$_2$O is a tracer of the origin and the evolution mechanisms of planets. For temperate, terrestrial planets, the presence of H$_2$O is of great significance as an indicator of habitable conditions. Being small and relatively cold, these planets and their atmospheres are the most challenging to observe, and therefore no atmospheric spectral signatures have so far been detected. Super-Earths -- planets lighter than ten M$_\oplus$ -- around later-type stars may provide our first opportunity to study spectroscopically the characteristics of such planets, as they are best suited for transit observations. Here we report the detection of an H$_2$O spectroscopic signature in the atmosphere of \planet\ -- an eight M$_\oplus$ planet in the habitable-zone of an M-dwarf -- with high statistical confidence (ADI = 5.0, $\sim$3.6$σ$). In addition, the derived mean molecular weight suggests an atmosphere still containing some hydrogen. The observations were recorded with the Hubble Space Telescope/WFC3 camera, and analysed with our dedicated, publicly available, algorithms. While the suitability of M-dwarfs to host habitable worlds is still under discussion, \planet\ offers an unprecedented opportunity to get insight into the composition and climate of habitable-zone planets. △ Less

Submitted 11 September, 2019; originally announced September 2019.

Comments: Published in Nature Astronomy

Impact of planetary mass uncertainties on exoplanet atmospheric retrievals

Authors: Quentin Changeat, Luke Keyte, Ingo P Waldmann, Giovanna Tinetti

Abstract: In current models used to interpret exoplanet atmospheric observations, the planet mass is treated as a prior and is estimated independently with external methods, such as RV or TTV techniques. This approach is necessary as available spectroscopic data do not have sufficient wavelength coverage and/or SNR to infer the planetary mass. We examine here the impact of mass uncertainties on spectral ret… ▽ More In current models used to interpret exoplanet atmospheric observations, the planet mass is treated as a prior and is estimated independently with external methods, such as RV or TTV techniques. This approach is necessary as available spectroscopic data do not have sufficient wavelength coverage and/or SNR to infer the planetary mass. We examine here the impact of mass uncertainties on spectral retrieval analyses for a host of atmospheric scenarios. Our approach is both analytical and numerical: we first use simple approximations to extract analytically the influence of each parameter to the wavelength-dependent transit depth. We then adopt a fully Bayesian retrieval model to quantify the propagation of the mass uncertainty onto other atmospheric parameters. We found that for clear-sky, gaseous atmospheres the posterior distributions are the same when the mass is known or retrieved. The retrieved mass is very accurate, with a precision of more than 10%, provided the wavelength coverage and S/N are adequate. When opaque clouds are included in the simulations, the uncertainties in the retrieved mass increase, especially for high altitude clouds. However atmospheric parameters such as the temperature and trace-gas abundances are unaffected by the knowledge of the mass. Secondary atmospheres are more challenging due to the higher degree of freedom for the atmospheric main component, which is unknown. For broad wavelength range and adequate SNR, the mass can still be retrieved accurately and precisely if clouds are not present, and so are all the other atmospheric/planetary parameters. When clouds are added, we find that the mass uncertainties may impact substantially the retrieval of the mean molecular weight: an independent characterisation of the mass would therefore be helpful to capture/confirm the main atmospheric constituent. △ Less

Submitted 9 October, 2020; v1 submitted 17 August, 2019; originally announced August 2019.

Comments: 19 pages, 12 figures, Accepted in ApJ

Mapping Saturn using deep learning

Authors: Ingo P. Waldmann, Caitlin A. Griffith

Abstract: Clouds and aerosols provide unique insight into the chemical and physical processes of gas-giant planets. Mapping and characterizing the spectral features indicative of the cloud structure and composition enables an understand-ing of a planet's energy budget, chemistry and atmospheric dynamics. Current space missions to Solar System planets produce high-quality datasets, yet the sheer amount of da… ▽ More Clouds and aerosols provide unique insight into the chemical and physical processes of gas-giant planets. Mapping and characterizing the spectral features indicative of the cloud structure and composition enables an understand-ing of a planet's energy budget, chemistry and atmospheric dynamics. Current space missions to Solar System planets produce high-quality datasets, yet the sheer amount of data obtained often prohibits detailed 'by hand' analyses. Current techniques rely mainly on two approaches: identifying the existence of spectral features by dividing the fluxes of two or more spectral channels, and performing full radiative transfer calculations for individual spectra. The first method is not sufficiently accurate and the second is not easily scalable to the entire planetary surface. Here we have developed a deep learning algorithm, PlanetNet, that is able to quickly and accurately map spatial and spectral features across large, heterogeneous areas of a planet. We use PlanetNet to delineate the major components of the 2008 storm on Saturn, enhancing the scope of the area previously studied and indicating regions that can be probed more deeply with radiative transfer models. Our spectral-component maps indicate compositional and cloud variations of the vast region affected by the storm, showing regions of vertical upwelling, and diminished clouds at the centre of compact substorms △ Less

Submitted 29 April, 2019; originally announced April 2019.

Comments: Nature Astronomy, 28 pages 18 figures

Pushing the Limits of Exoplanet Discovery via Direct Imaging with Deep Learning

Authors: Kai Hou Yip, Nikolaos Nikolaou, Piero Coronica, Angelos Tsiaras, Billy Edwards, Quentin Changeat, Mario Morvan, Beth Biller, Sasha Hinkley, Jeffrey Salmond, Matthew Archer, Paul Sumption, Elodie Choquet, Remi Soummer, Laurent Pueyo, Ingo P. Waldmann

Abstract: Further advances in exoplanet detection and characterisation require sampling a diverse population of extrasolar planets. One technique to detect these distant worlds is through the direct detection of their thermal emission. The so-called direct imaging technique, is suitable for observing young planets far from their star. These are very low signal-to-noise-ratio (SNR) measurements and limited g… ▽ More Further advances in exoplanet detection and characterisation require sampling a diverse population of extrasolar planets. One technique to detect these distant worlds is through the direct detection of their thermal emission. The so-called direct imaging technique, is suitable for observing young planets far from their star. These are very low signal-to-noise-ratio (SNR) measurements and limited ground truth hinders the use of supervised learning approaches. In this paper, we combine deep generative and discriminative models to bypass the issues arising when directly training on real data. We use a Generative Adversarial Network to obtain a suitable dataset for training Convolutional Neural Network classifiers to detect and locate planets across a wide range of SNRs. Tested on artificial data, our detectors exhibit good predictive performance and robustness across SNRs. To demonstrate the limits of the detectors, we provide maps of the precision and recall of the model per pixel of the input image. On real data, the models can re-confirm bright source detections. △ Less

Submitted 28 January, 2020; v1 submitted 12 April, 2019; originally announced April 2019.

Comments: 16 Pages, 6 Figures, 3 Tables, Presented in ECML-PKDD 2019

Integrating light-curve and atmospheric modelling of transiting exoplanets

Authors: Kai Hou Yip, Ingo P. Waldmann, Angelos Tsiaras, Giovanna Tinetti

Abstract: Spectral retrieval techniques are currently our best tool to interpret the observed exoplanet atmospheric data. Said techniques retrieve the optimal atmospheric components and parameters by identifying the best fit to an observed transmission/emission spectrum. Over the past decade, our understanding of remote worlds in our galaxy has flourished thanks to the use of increasingly sophisticated spec… ▽ More Spectral retrieval techniques are currently our best tool to interpret the observed exoplanet atmospheric data. Said techniques retrieve the optimal atmospheric components and parameters by identifying the best fit to an observed transmission/emission spectrum. Over the past decade, our understanding of remote worlds in our galaxy has flourished thanks to the use of increasingly sophisticated spectral retrieval techniques and the collective effort of the community working on exoplanet atmospheric models. A new generation of instruments in space and from the ground is expected to deliver higher quality data in the next decade, it is therefore paramount to upgrade current models and improve their reliability, completeness and numerical speed with which they can be run. In this paper, we address the issue of reliability of the results provided by retrieval models in the presence of systematics of unknown origin. More specifically, we demonstrate that if we fit directly individual light-curves at different wavelengths (L-retrieval), instead of fitting transit or eclipse depths, as it is currently done (S-retrieval), the results obtained are more robust against astrophysical and instrumental noise. This new approach is tested, in particular, when discrepant simulated observations from HST/WFC3 and Spitzer/IRAC are combined. We find that while S-retrievals converge to an incorrect solution without any warning, L-retrievals are able to identify potential discrepancies between the data-sets. △ Less

Submitted 19 October, 2020; v1 submitted 12 November, 2018; originally announced November 2018.

Comments: 14 pages, 15 figures,Submitted to ApJ

Journal ref: The Astronomical Journal, Volume 160 (2020), Issue 4, id.171, 13 pp

ExoGAN: Retrieving Exoplanetary Atmospheres Using Deep Convolutional Generative Adversarial Networks

Authors: Tiziano Zingales, Ingo Peter Waldmann

Abstract: Atmospheric retrievals on exoplanets usually involve computationally intensive Bayesian sampling methods. Large parameter spaces and increasingly complex atmospheric models create a computational bottleneck forcing a trade-off between statistical sampling accuracy and model complexity. It is especially true for upcoming JWST and ARIEL observations. We introduce ExoGAN, the Exoplanet Generative Adv… ▽ More Atmospheric retrievals on exoplanets usually involve computationally intensive Bayesian sampling methods. Large parameter spaces and increasingly complex atmospheric models create a computational bottleneck forcing a trade-off between statistical sampling accuracy and model complexity. It is especially true for upcoming JWST and ARIEL observations. We introduce ExoGAN, the Exoplanet Generative Adversarial Network, a new deep learning algorithm able to recognise molecular features, atmospheric trace-gas abundances and planetary parameters using unsupervised learning. Once trained, ExoGAN is widely applicable to a large number of instruments and planetary types. The ExoGAN retrievals constitute a significant speed improvement over traditional retrievals and can be used either as a final atmospheric analysis or provide prior constraints to subsequent retrieval. △ Less

Submitted 9 October, 2018; v1 submitted 7 June, 2018; originally announced June 2018.

Comments: 19 pages, 17 figures, 7 tables

The Transiting Exoplanet Community Early Release Science Program for JWST

Authors: Jacob L. Bean, Kevin B. Stevenson, Natalie M. Batalha, Zachory Berta-Thompson, Laura Kreidberg, Nicolas Crouzet, Björn Benneke, Michael R. Line, David K. Sing, Hannah R. Wakeford, Heather A. Knutson, Eliza M. -R. Kempton, Jean-Michel Désert, Ian Crossfield, Natasha E. Batalha, Julien de Wit, Vivien Parmentier, Joseph Harrington, Julianne I. Moses, Mercedes Lopez-Morales, Munazza K. Alam, Jasmina Blecic, Giovanni Bruno, Aarynn L. Carter, John W. Chapman , et al. (77 additional authors not shown)

Abstract: The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other… ▽ More The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission. △ Less

Submitted 3 September, 2018; v1 submitted 13 March, 2018; originally announced March 2018.

Comments: PASP in press

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

Authors: Olivia Venot, Benjamin Drummond, Yamila Miguel, Ingo P. Waldmann, Enzo Pascale, Tiziano Zingales

Abstract: Since the discovery of the first extrasolar planet more than twenty years ago, we have discovered nearly four thousand planets orbiting stars other than the Sun. 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 planet… ▽ More Since the discovery of the first extrasolar planet more than twenty years ago, we have discovered nearly four thousand planets orbiting stars other than the Sun. 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. △ Less

Submitted 2 May, 2018; v1 submitted 22 November, 2017; originally announced November 2017.

Comments: Accepted for publication in Experimental Astronomy, ARIEL Special Issue

VUV-absorption cross section of carbon dioxide from 150 to 800 K and applications to warm exoplanetary atmospheres

Authors: Olivia Venot, Yves Bénilan, Nicolas Fray, Marie-Claire Gazeau, Franck Lefèvre, Et-touhami Es-sebbar, Eric Hébrard, Martin Schwell, Chiheb Bahrini, Franck Montmessin, Maxence Lefèvre, Ingo P. Waldmann

Abstract: Most exoplanets detected so far have atmospheric T significantly higher than 300K. Often close to their star, they receive an intense UV photons flux that triggers important photodissociation processes. The T dependency of VUV absorption cross sections are poorly known, leading to an undefined uncertainty in atmospheric models. Similarly, data measured at low T similar to that of the high atmosphe… ▽ More Most exoplanets detected so far have atmospheric T significantly higher than 300K. Often close to their star, they receive an intense UV photons flux that triggers important photodissociation processes. The T dependency of VUV absorption cross sections are poorly known, leading to an undefined uncertainty in atmospheric models. Similarly, data measured at low T similar to that of the high atmosphere of Mars, Venus, and Titan are often lacking. Our aim is to quantify the T dependency of the abs. cross section of important molecules in planetary atmospheres. We want to provide both high-resolution data at T prevailing in these media and a simple parameterization of the absorption in order to simplify its use in photochemical models. This study focuses on carbon dioxide. We performed experimental measurements of CO$_2$ absorption cross section with synchrotron radiation for the wavelength range (115--200nm). For longer wavelengths (195--230nm), we used a deuterium lamp and a 1.5m Jobin-Yvon spectrometer. We used these data in our 1D thermo-photochemical model in order to study their impact on the predicted atmospheric compositions. The cross section of CO$_2$ increases with T. It can be separated in two parts: a continuum and a fine structure superimposed on the continuum. The variation of the continuum of absorption can be represented by the sum of three gaussian functions. Using data at high T in thermo-photochemical models modifies significantly the abundance and the photodissociation rates of many species, in addition to CO$_2$, such as methane and ammonia. These deviations have an impact on synthetic transmission spectra, leading to variations of up to 5 ppm. We present a full set of HR ($Δλ$=0.03nm) absorption cross sections of CO$_2$ from 115 to 230nm for T ranging from 150 to 800K. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: 14 pages, 17 figures, accepted for pulication in A&A

Journal ref: A&A 609, A34 (2018)

A hybrid line list for CH$_4$ and hot methane continuum

Authors: Sergei N. Yurchenko, David S. Amundsen, Jonathan Tennyson, Ingo P. Waldmann

Abstract: Molecular line lists are important for modelling absorption and emission processes in atmospheres of different astronomical objects, such as cool stars and exoplanets. In order to be applicable for high temperatures, line lists for molecules like methane must contain billions of transitions, which makes their direct (line-by-line) application in radiative transfer calculations impracticable. Here… ▽ More Molecular line lists are important for modelling absorption and emission processes in atmospheres of different astronomical objects, such as cool stars and exoplanets. In order to be applicable for high temperatures, line lists for molecules like methane must contain billions of transitions, which makes their direct (line-by-line) application in radiative transfer calculations impracticable. Here we suggest a new, hybrid line list format to mitigate this problem, based on the idea of temperature-dependent absorption continuum. The line list is partitioned into a large set of relatively weak lines and a small set of important, stronger lines. The weaker lines are then used either to construct a temperature-dependent (but pressure-independent) set of intensity cross sections or are blended into a greatly reduced set of super-lines. The strong lines are kept in the form of temperature independent Einstein A coefficients. A line list for methane is constructed as a combination of 17 million strong absorption lines relative to the reference absorption spectra and a background methane continuum in two temperature-dependent forms, of cross sections and super-lines. This approach eases the use of large high temperature line lists significantly as the computationally expensive calculation of pressure dependent profiles only need to be performed for a relatively small number of lines. Both the line list and cross sections were generated using a new 34 billion methane line list (34to10), which extends the 10to10 line list to higher temperatures (up to 2000 K). The new hybrid scheme can be applied to any large line lists containing billions of transitions. We recommend to use super-lines generated on a high resolution grid based on resolving power (R = 1,000,000) to model the molecular continuum as a more flexible alternative to the temperature dependent cross sections. △ Less

Submitted 18 June, 2017; originally announced June 2017.

Journal ref: A&A 605, A95 (2017)

Community targets for JWST's early release science program: evaluation of WASP-63b

Authors: Brian M. Kilpatrick, Patricio E. Cubillos, Kevin B. Stevenson, Nikole K. Lewis, Hannah Wakeford, Ryan J. Macdonald, Nikku Madhusudhan, Jasmina Blecic, Giovanni Bruno, Adam Burrows, Drake Deming, Kevin Heng, Michael R. Line, Caroline V. Morley, Vivien Parmentier, Gregory S. Tucker, Jeff A. Valenti, Ingo P. Waldmann, Jacob L. Bean, Charles Beichman, Jonathan Fraine, J. E. Krick, Joshua D. Lothringer, Avi M. Mandell

Abstract: We present observations of WASP-63b by the Hubble Space Telescope (HST) as part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet for JWST Early Release Science". WASP-63b is one of the community targets under consideration for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. We present a spectrum derived from a single observation by HST Wide Fiel… ▽ More We present observations of WASP-63b by the Hubble Space Telescope (HST) as part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet for JWST Early Release Science". WASP-63b is one of the community targets under consideration for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. We present a spectrum derived from a single observation by HST Wide Field Camera 3 in the near infrared. We engaged groups across the transiting exoplanet community to participate in the analysis of the data and present results from each. There is general agreement amongst all results that we find an H2O absorption feature with 3.5-4.0 sigma significance. However, the feature is muted in comparison to a clear atmosphere at solar composition. Although the detection of the water feature is robust, the reasons for the muting of this feature are ambiguous due to a degeneracy between clouds and composition. The data does not yield robust detections of any molecular species other than H2O. The group was motivated to perform an additional set of retrieval exercises to investigate an apparent bump in the spectrum at ~ 1.55 um. We explore possible disequilibrium chemistry and find this feature is consistent with super-solar HCN abundance but it is questionable if the required mixing ratio of HCN is chemically and physically plausible. The ultimate goal of this study is to vet WASP-63b as a potential community target to best demonstrate the capabilities and systematics of JWST instruments for transiting exoplanet science. In the case of WASP-63b, the presence of a detectable water feature indicates that WASP-63b remains a plausible target for ERS observations. △ Less

Submitted 30 August, 2018; v1 submitted 24 April, 2017; originally announced April 2017.

Comments: repaired broken citations

A population study of gaseous exoplanets

Authors: A. Tsiaras, I. P. Waldmann, T. Zingales, M. Rocchetto, G. Morello, M. Damiano, K. Karpouzas, G. Tinetti, L. K. McKemmish, J. Tennyson, S. N. Yurchenko

Abstract: We present here the analysis of 30 gaseous extrasolar planets, with temperatures between 600 and 2400 K and radii between 0.35 and 1.9 $R_\mathrm{Jup}$. The quality of the HST/WFC3 spatially scanned data combined with our specialized analysis tools allow us to study the largest and most self-consistent sample of exoplanetary transmission spectra to date and examine the collective behavior of warm… ▽ More We present here the analysis of 30 gaseous extrasolar planets, with temperatures between 600 and 2400 K and radii between 0.35 and 1.9 $R_\mathrm{Jup}$. The quality of the HST/WFC3 spatially scanned data combined with our specialized analysis tools allow us to study the largest and most self-consistent sample of exoplanetary transmission spectra to date and examine the collective behavior of warm and hot gaseous planets rather than isolated case-studies. We define a new metric, the Atmospheric Detectability Index (ADI) to evaluate the statistical significance of an atmospheric detection and find statistically significant atmospheres around 16 planets out of the 30 analysed. For most of the Jupiters in our sample, we find the detectability of their atmospheres to be dependent on the planetary radius but not on the planetary mass. This indicates that planetary gravity plays a secondary role in the state of gaseous planetary atmospheres. We detect the presence of water vapour in all of the statistically detectable atmospheres, and we cannot rule out its presence in the atmospheres of the others. In addition, TiO and/or VO signatures are detected with 4$σ$ confidence in WASP-76 b, and they are most likely present in WASP-121 b. We find no correlation between expected signal-to-noise and atmospheric detectability for most targets. This has important implications for future large-scale surveys. △ Less

Submitted 17 March, 2018; v1 submitted 18 April, 2017; originally announced April 2017.

Comments: 14 pages, 12 figures, 3 tables, published in AJ

Exploring Biases of Atmospheric Retrievals in Simulated JWST Transmission Spectra of Hot Jupiters

Authors: M. Rocchetto, I. P. Waldmann, O. Venot, P. -O. Lagage, G. Tinetti

Abstract: With a scheduled launch in October 2018, the James Webb Space Telescope (JWST) is expected to revolutionise 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 cur… ▽ More With a scheduled launch in October 2018, the James Webb Space Telescope (JWST) is expected to revolutionise 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 parametrized 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 parametrized profile are usually within one sigma of the true state, and we found the retrieved uncertainties to be generally larger compared to the isothermal approximation. Interestingly, we found that using the parametrized temperature profile we could place tight constraints on the temperature structure. This opens the possibility to characterize 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. △ Less

Submitted 10 October, 2016; originally announced October 2016.

Comments: 16 pages, 10 figures. Accepted for publication in The Astrophysical Journal

Repeatability of Spitzer/IRAC exoplanetary eclipses with Independent Component Analysis

Authors: Giuseppe Morello, Ingo P. Waldmann, Giovanna Tinetti

Abstract: The research of effective and reliable detrending methods for Spitzer data is of paramount importance for the characterization of exoplanetary atmospheres. To date, the totality of exoplanetary observations in the mid- and far-infrared, at wavelengths $>$3 $μ$m, have been taken with Spitzer. In some cases, in the past years, repeated observations and multiple reanalyses of the same datasets led to… ▽ More The research of effective and reliable detrending methods for Spitzer data is of paramount importance for the characterization of exoplanetary atmospheres. To date, the totality of exoplanetary observations in the mid- and far-infrared, at wavelengths $>$3 $μ$m, have been taken with Spitzer. In some cases, in the past years, repeated observations and multiple reanalyses of the same datasets led to discrepant results, raising questions about the accuracy and reproducibility of such measurements. Morello et al. 2014, 2015 proposed a blind-source separation method based on the Independent Component Analysis of pixel time series (pixel-ICA) to analyze IRAC data, obtaining coherent results when applied to repeated transit observations previously debated in the literature. Here we introduce a variant to pixel-ICA through the use of wavelet transform, wavelet pixel-ICA, which extends its applicability to low-S/N cases. We describe the method and discuss the results obtained over twelve eclipses of the exoplanet XO3b observed during the "Warm Spitzer" era in the 4.5 $μ$m band. The final results will be reported also in Ingalls et al. (in prep.), together with results obtained with other detrending methods, and over ten synthetic eclipses that were analyzed for the "IRAC Data Challenge 2015". Our results are consistent within 1 $σ$ with the ones reported in Wong et al. 2014. The self-consistency of individual measurements of eclipse depth and phase curve slope over a span of more than three years proves the stability of Warm Spitzer/IRAC photometry within the error bars, at the level of 1 part in 10$^4$ in stellar flux. △ Less

Submitted 25 February, 2016; v1 submitted 15 January, 2016; originally announced January 2016.

Comments: accepted by ApJ (updated reference to Zellem et al. and acknowledgments)

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

Authors: A. Tsiaras, M. Rocchetto, I. P. Waldmann, O. Venot, R. Varley, G. Morello, M. Damiano, G. Tinetti, E. J. Barton, S. N. Yurchenko, J. Tennyson

Abstract: 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 lengt… ▽ More 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$σ$ 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 $μ$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. △ Less

Submitted 6 February, 2016; v1 submitted 28 November, 2015; originally announced November 2015.

Comments: 10 pages, 10 figures, 4 tables, Accepted for publication in ApJ

Dreaming of atmospheres

Authors: I. P. Waldmann

Abstract: Here we introduce the RobERt (Robotic Exoplanet Recognition) algorithm for the classification of exoplanetary emission spectra. Spectral retrievals of exoplanetary atmospheres frequently requires the preselection of molecular/atomic opacities to be defined by the user. In the era of open-source, automated and self-sufficient retrieval algorithms, manual input should be avoided. User dependent inpu… ▽ More Here we introduce the RobERt (Robotic Exoplanet Recognition) algorithm for the classification of exoplanetary emission spectra. Spectral retrievals of exoplanetary atmospheres frequently requires the preselection of molecular/atomic opacities to be defined by the user. In the era of open-source, automated and self-sufficient retrieval algorithms, manual input should be avoided. User dependent input could, in worst case scenarios, lead to incomplete models and biases in the retrieval. The RobERt algorithm is based on deep belief neural (DBN) networks trained to accurately recognise molecular signatures for a wide range of planets, atmospheric thermal profiles and compositions. Reconstructions of the learned features, also referred to as `dreams' of the network, indicate good convergence and an accurate representation of molecular features in the DBN. Using these deep neural networks, we work towards retrieval algorithms that themselves understand the nature of the observed spectra, are able to learn from current and past data and make sensible qualitative preselections of atmospheric opacities to be used for the quantitative stage of the retrieval process. △ Less

Submitted 11 February, 2016; v1 submitted 26 November, 2015; originally announced November 2015.

Comments: ApJ accepted

A new approach to analysing HST spatial scans: the transmission spectrum of HD 209458 b

Authors: A. Tsiaras, I. P. Waldmann, M. Rocchetto, R. Varley, G. Morello, M. Damiano, G. Tinetti

Abstract: The Wide Field Camera 3 (WFC3) on Hubble Space Telescope (HST) is currently one of the most widely used instruments for observing exoplanetary atmospheres, especially with the use of the spatial scanning technique. An increasing number of exoplanets have been studied using this technique as it enables the observation of bright targets without saturating the sensitive detectors. In this work we pre… ▽ More The Wide Field Camera 3 (WFC3) on Hubble Space Telescope (HST) is currently one of the most widely used instruments for observing exoplanetary atmospheres, especially with the use of the spatial scanning technique. An increasing number of exoplanets have been studied using this technique as it enables the observation of bright targets without saturating the sensitive detectors. In this work we present a new pipeline for analyzing the data obtained with the spatial scanning technique, starting from the raw data provided by the instrument. In addition to commonly used correction techniques, we take into account the geometric distortions of the instrument, whose impact may become important when combined to the scanning process. Our approach can improve the photometric precision for existing data and also push further the limits of the spatial scanning technique, as it allows the analysis of even longer spatial scans. As an application of our method and pipeline, we present the results from a reanalysis of the spatially scanned transit spectrum of HD 209458 b. We calculate the transit depth per wavelength channel with an average relative uncertainty of 40 ppm. We interpret the final spectrum with T-Rex, our fully Bayesian spectral retrieval code, which confirms the presence of water vapor and clouds in the atmosphere of HD 209458 b. The narrow wavelength range limits our ability to disentangle the degeneracies between the fitted atmospheric parameters. Additional data over a broader spectral range are needed to address this issue. △ Less

Submitted 15 October, 2016; v1 submitted 24 November, 2015; originally announced November 2015.

Comments: 13 pages, 15 figures, 7 tables, Accepted for publication in ApJ

Tau-REx II: Retrieval of emission spectra

Authors: Ingo P. Waldmann, Marco Rocchetto, Giovanna Tinetti, Emma J. Barton, Sergey N. Yurchenko, Jonathan Tennyson

Abstract: Tau-REx (Tau Retrieval of Exoplanets) is a novel, fully Bayesian atmospheric retrieval code custom built for extrasolar atmospheres. In Waldmann et al. (2015) the transmission spectroscopic case was introduced, here we present the emission spectroscopy spectral retrieval for the Tau-REx framework. Compared to transmission spectroscopy, the emission case is often significantly more degenerate due t… ▽ More Tau-REx (Tau Retrieval of Exoplanets) is a novel, fully Bayesian atmospheric retrieval code custom built for extrasolar atmospheres. In Waldmann et al. (2015) the transmission spectroscopic case was introduced, here we present the emission spectroscopy spectral retrieval for the Tau-REx framework. Compared to transmission spectroscopy, the emission case is often significantly more degenerate due to the need to retrieve the full atmospheric temperature-pressure (TP) profile. This is particularly true in the case of current measurements of exoplanetary atmospheres, which are either of low signal-to-noise, low spectral resolution or both. Here we present a new way of combining two existing approaches to the modelling of the said TP profile: 1) the parametric profile, where the atmospheric TP structure is analytically approximated by a few model parameters, 2) the Layer-by-Layer approach, where individual atmospheric layers are modelled. Both these approaches have distinct advantages and disadvantages in terms of convergence properties and potential model biases. The Tau-REx hybrid model presented here is a new two-stage TP profile retrieval, which combines the robustness of the analytic solution with the accuracy of the Layer-by-Layer approach. The retrieval process is demonstrated using simulations of the hot-Jupiter WASP-76b and the hot SuperEarth 55 Cnc e, as well as on the secondary eclipse measurements of HD189733b. △ Less

Submitted 9 September, 2015; v1 submitted 30 August, 2015; originally announced August 2015.

Comments: ApJ accepted

Revisiting Spitzer transit observations with Independent Component Analysis: new results for the GJ436 system

Authors: G. Morello, I. P. Waldmann, G. Tinetti, I. D. Howarth, G. Micela, F. Allard

Abstract: We analyzed four Spitzer/IRAC observations at 3.6 and 4.5 μm of the primary transit of the exoplanet GJ436b, by using blind source separation techniques. These observations are important to investigate the atmospheric composition of the planet GJ436b. Previous analyses claimed strong inter-epoch variations of the transit parameters due to stellar variability, casting doubts on the possibility to e… ▽ More We analyzed four Spitzer/IRAC observations at 3.6 and 4.5 μm of the primary transit of the exoplanet GJ436b, by using blind source separation techniques. These observations are important to investigate the atmospheric composition of the planet GJ436b. Previous analyses claimed strong inter-epoch variations of the transit parameters due to stellar variability, casting doubts on the possibility to extract conclusively an atmospheric signal; those analyses also reported discrepant results, hence the necessity of this reanalysis. The method we used has been proposed in Morello et al. (2014) to analyze 3.6 μm transit light-curves of the hot Jupiter HD189733b; it performes an Independent Component Analysis (ICA) on a set of pixel-light-curves, i.e. time series read by individual pixels, from the same photometric observation. Our method only assumes the independence of instrumental and astrophysical signals, and therefore guarantees a higher degree of objectivity compared to parametric detrending techniques published in the literature. The datasets we analyzed in this paper represent a more challenging test compared to the previous ones. Contrary to previous results reported in the literature, our results (1) do not support any detectable inter-epoch variations of orbital and stellar parameters, (2) are photometrically stable at the level 10e-4 in the IR, and (3) the transit depth measurements at the two wavelengths are consistent within 1σ. We also (4) detect a possible transit duration variation (TDV) of 80 s (2 σ significance level), that has not been pointed out in the literature, and (5) confirm no transit timing variations (TTVs) >30 s. △ Less

Submitted 23 January, 2015; originally announced January 2015.

Comments: accepted for publication in The Astrophysical Journal

The ground-based H, K, and L-band absolute emission spectra of HD 209458b

Authors: Robert T. Zellem, Caitlin A. Griffith, Pieter Deroo, Mark R. Swain, Ingo P. Waldmann

Abstract: Here we explore the capabilities of NASA's 3.0 meter Infrared Telescope Facility (IRTF) and SpeX spectrometer and the 5.08 meter Hale telescope with the TripleSpec spectrometer with near-infrared H, K, and L-band measurements of HD 209458b's secondary eclipse. Our IRTF/SpeX data are the first absolute L-band spectroscopic emission measurements of any exoplanet other than the hot Jupiter HD 189733b… ▽ More Here we explore the capabilities of NASA's 3.0 meter Infrared Telescope Facility (IRTF) and SpeX spectrometer and the 5.08 meter Hale telescope with the TripleSpec spectrometer with near-infrared H, K, and L-band measurements of HD 209458b's secondary eclipse. Our IRTF/SpeX data are the first absolute L-band spectroscopic emission measurements of any exoplanet other than the hot Jupiter HD 189733b. Previous measurements of HD 189733b's L-band indicate bright emission hypothesized to result from non-LTE CH$_{4}$ $ν_{3}$ fluorescence. We do not detect a similar bright 3.3 $μ$m feature to ~3$σ$, suggesting that fluorescence does not need to be invoked to explain HD 209458b's L-band measurements. The validity of our observation and reduction techniques, which decrease the flux variance by up to 2.8 orders of magnitude, is reinforced by 1$σ$ agreement with existent $Hubble$/NICMOS and $Spitzer$/IRAC1 observations which overlap the H, K, and L-bands, suggesting that both IRTF/SpeX and Palomar/TripleSpec can measure an exoplanet's emission with high precision. △ Less

Submitted 29 September, 2014; originally announced September 2014.

Comments: 10 pages, 12 figures, 4 tables, accepted to ApJ, Table 4 is published in its entirety in the electronic edition of ApJ

Tau-REx I: A next generation retrieval code for exoplanetary atmospheres

Authors: Ingo P. Waldmann, Giovanna Tinetti, Marco Rocchetto, Emma J. Barton, Sergey N. Yurchenko, Jonathan Tennyson

Abstract: Spectroscopy of exoplanetary atmospheres has become a well established method for the characterisation of extrasolar planets. We here present a novel inverse retrieval code for exoplanetary atmospheres. TauRex (Tau Retrieval for Exoplanets) is a line-by-line radiative transfer fully Bayesian retrieval framework. TauRex includes the following features: 1) the optimised use of molecular line-lists f… ▽ More Spectroscopy of exoplanetary atmospheres has become a well established method for the characterisation of extrasolar planets. We here present a novel inverse retrieval code for exoplanetary atmospheres. TauRex (Tau Retrieval for Exoplanets) is a line-by-line radiative transfer fully Bayesian retrieval framework. TauRex includes the following features: 1) the optimised use of molecular line-lists from the Exomol project; 2) an unbiased atmospheric composition prior selection, through custom built pattern recognition software; 3) the use of two independent algorithms to fully sample the Bayesian likelihood space: nested sampling as well as a more classical Markov Chain Monte Carlo approach; 4) iterative Bayesian parameter and model selection using the full Bayesian Evidence as well as the Savage-Dickey Ratio for nested models, and 5) the ability to fully map very large parameter spaces through optimal code parallelisation and scalability to cluster computing. In this publication we outline the TauRex framework and demonstrate, using a theoretical hot-Jupiter transmission spectrum, the parameter retrieval and model selection. We investigate the impact of Signal-to-Noise and spectral resolution on the retrievability of individual model parameters, both in terms of error bars on the temperature and molecular mixing ratios as well as its effect on the model's global Bayesian evidence. △ Less

Submitted 6 April, 2015; v1 submitted 8 September, 2014; originally announced September 2014.

Comments: ApJ published

Rossiter-McLaughlin Observations of 55 Cnc e

Authors: Mercedes Lopez-Morales, Amaury H. M. J. Triaud, Florian Rodler, Xavier Dumusque, Lars A. Buchhave, A. Harutyunyan, Sergio Hoyer, Roi Alonso, Michael Gillon, Nathan A. Kaib, David W. Latham, Christophe Lovis, Francesco Pepe, Didier Queloz, Sean N. Raymond, Damien Segransan, Ingo P. Waldmann, Stephane Udry

Abstract: We present Rossiter-McLaughlin observations of the transiting super-Earth 55 Cnc e collected during six transit events between January 2012 and November 2013 with HARPS and HARPS-N. We detect no radial-velocity signal above 35 cm/s (3-sigma) and confine the stellar v sin i to 0.2 +/- 0.5 km/s. The star appears to be a very slow rotator, producing a very low amplitude Rossiter-McLaughlin effect. Gi… ▽ More We present Rossiter-McLaughlin observations of the transiting super-Earth 55 Cnc e collected during six transit events between January 2012 and November 2013 with HARPS and HARPS-N. We detect no radial-velocity signal above 35 cm/s (3-sigma) and confine the stellar v sin i to 0.2 +/- 0.5 km/s. The star appears to be a very slow rotator, producing a very low amplitude Rossiter-McLaughlin effect. Given such a low amplitude, the Rossiter-McLaughlin effect of 55 Cnc e is undetected in our data, and any spin-orbit angle of the system remains possible. We also performed Doppler tomography and reach a similar conclusion. Our results offer a glimpse of the capacity of future instrumentation to study low amplitude Rossiter-McLaughlin effects produced by super-Earths. △ Less

Submitted 8 August, 2014; originally announced August 2014.

Comments: Accepted for publication in ApJ Letters

EChOSim: The Exoplanet Characterisation Observatory software simulator

Authors: E. Pascale, I. P. Waldmann, C. J. MacTavish, A. Papageorgiou, A. Amaral-Rogers, R. Varley, V. Coudé de Foresto, M. J. Griffin, M. Ollivier, S. Sarkar, L. Spencer, B. M. Swinyard, M. Tessenyi, G. Tinetti

Abstract: 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 EChO has to detect and characterize the atmospheres of transiting exoplanets, and through this revolutionize the knowledge we have of the Milky Way and of our place in the Galaxy. Here we discuss the details of the EChOSim implement… ▽ More 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 EChO has to detect and characterize the atmospheres of transiting exoplanets, and through this revolutionize the knowledge we have of the Milky Way and of our place in the Galaxy. 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 for exoplanet transits, including instruments on the Hubble Space Telescope and Spitzer, or 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. △ Less

Submitted 16 June, 2014; originally announced June 2014.

A new look at Spitzer primary transit observations of the exoplanet HD189733b

Authors: Giuseppe Morello, Ingo P. Waldmann, Giovanna Tinetti, Giovanni Peres, Giuseppina Micela, Ian D. Howarth

Abstract: Blind source separation techniques are used to reanalyse two exoplanetary transit lightcurves of the exoplanet HD189733b recorded with the IR camera IRAC on board the Spitzer Space Telescope at 3.6$μ$m during the "cold" era. These observations, together with observations at other IR wavelengths, are crucial to characterise the atmosphere of the planet HD189733b. Previous analyses of the same datas… ▽ More Blind source separation techniques are used to reanalyse two exoplanetary transit lightcurves of the exoplanet HD189733b recorded with the IR camera IRAC on board the Spitzer Space Telescope at 3.6$μ$m during the "cold" era. These observations, together with observations at other IR wavelengths, are crucial to characterise the atmosphere of the planet HD189733b. Previous analyses of the same datasets reported discrepant results, hence the necessity of the reanalyses. The method we used here is based on the Independent Component Analysis (ICA) statistical technique, which ensures a high degree of objectivity. The use of ICA to detrend single photometric observations in a self-consistent way is novel in the literature. The advantage of our reanalyses over previous work is that we do not have to make any assumptions on the structure of the unknown instrumental systematics. Such "admission of ignorance" may result in larger error bars than reported in the literature, up to a factor $1.6$. This is a worthwhile trade-off for much higher objectivity, necessary for trustworthy claims. Our main results are (1) improved and robust values of orbital and stellar parameters, (2) new measurements of the transit depths at 3.6$μ$m, (3) consistency between the parameters estimated from the two observations, (4) repeatability of the measurement within the photometric level of $\sim 2 \times 10^{-4}$ in the IR, (5) no evidence of stellar variability at the same photometric level within 1 year. △ Less

Submitted 12 March, 2014; originally announced March 2014.

Comments: 43 pages, 18 figures

Journal ref: Morello et al. 2014, ApJ, 786, 22

Data analysis Pipeline for EChO end-to-end simulations

Authors: Ingo P. Waldmann, Enzo Pascale

Abstract: Atmospheric spectroscopy of extrasolar planets is an intricate business. Atmospheric signatures typically require a photometric precision of $1 \times 10^{-4}$ in flux over several hours. Such precision demands high instrument stability as well as an understanding of stellar variability and an optimal data reduction and removal of systematic noise. In the context of the EChO mission concept, we he… ▽ More Atmospheric spectroscopy of extrasolar planets is an intricate business. Atmospheric signatures typically require a photometric precision of $1 \times 10^{-4}$ in flux over several hours. Such precision demands high instrument stability as well as an understanding of stellar variability and an optimal data reduction and removal of systematic noise. In the context of the EChO mission concept, we here discuss the data reduction and analysis pipeline developed for the EChO end-to-end simulator EChOSim. We present and discuss the step by step procedures required in order to obtain the final exoplanetary spectrum from the EChOSim`raw data' using a simulated observation of the secondary eclipse of the hot-Neptune 55 Cnc e. △ Less

Submitted 18 February, 2014; originally announced February 2014.

Comments: Submitted to Experimental Astronomy (EChO special issue)

On signals faint and sparse: the ACICA algorithm for blind de-trending of Exoplanetary transits with low signal-to-noise

Authors: I. P. Waldmann

Abstract: Independent Component Analysis (ICA) has recently been shown to be a promising new path in data analysis and de-trending of exoplanetary time series signals. Such approaches do not require or assume any prior or auxiliary knowledge on the data or instrument in order to de-convolve the astrophysical light curve signal from instrument or stellar systematic noise. These methods are often known as 'bl… ▽ More Independent Component Analysis (ICA) has recently been shown to be a promising new path in data analysis and de-trending of exoplanetary time series signals. Such approaches do not require or assume any prior or auxiliary knowledge on the data or instrument in order to de-convolve the astrophysical light curve signal from instrument or stellar systematic noise. These methods are often known as 'blind source separation' (BSS) algorithms. Unfortunately all BSS methods suffer from a amplitude and sign ambiguity of their de-convolved components which severely limits these methods in low signal-to-noise (S/N) observations where their scalings cannot be determined otherwise. Here we present a novel approach to calibrate ICA using sparse wavelet calibrators. The Amplitude Calibrated Independent Component Analysis (ACICA) allows for the direct retrieval of the independent components' scalings and the robust de-trending of low S/N data. Such an approach gives us an unique and unprecedented insight in the underlying morphology of a data set, making this method a powerful tool for exoplanetary data de-trending and signal diagnostics. △ Less

Submitted 27 February, 2013; originally announced February 2013.

Comments: submitted to ApJ

Photometric stability analysis of the Exoplanet Characterisation Observatory

Authors: I. P. Waldmann, E. Pascale, B. Swinyard, G. Tinetti, A. Amaral-Rogers, L. Spencer, M. Tessenyi, M. Ollivier, V. Coudé du Foresto

Abstract: Photometric stability is a key requirement for time-resolved spectroscopic observations of transiting extrasolar planets. In the context of the Exoplanet Characterisation Observatory (EChO) mission design, we here present and investigate means of translating spacecraft pointing instabilities as well as temperature fluctuation of its optical chain into an overall error budget of the exoplanetary sp… ▽ More Photometric stability is a key requirement for time-resolved spectroscopic observations of transiting extrasolar planets. In the context of the Exoplanet Characterisation Observatory (EChO) mission design, we here present and investigate means of translating spacecraft pointing instabilities as well as temperature fluctuation of its optical chain into an overall error budget of the exoplanetary spectrum to be retrieved. Given the instrument specifications as of date, we investigate the magnitudes of these photometric instabilities in the context of simulated observations of the exoplanet HD189733b secondary eclipse. △ Less

Submitted 26 February, 2013; originally announced February 2013.

Comments: submitted to MNRAS