How do we optimally sample model grids of exoplanet spectra?

Abstract

The construction and implementation of atmospheric model grids is a popular tool in exoplanet characterisation. These typically vary a number of parameters linearly, containing one model for every combination of parameter values. Here we investigate alternative methods of sampling parameters, including random sampling and Latin hypercube (LH) sampling, and how these compare to linearly sampled grids. We use a random forest to analyse the performance of these grids for two different models, as well as investigate the information content of the particular model grid from Goyal et al. 2019. We also use nested-sampling to implement mock atmospheric retrievals on simulated JWST transmission spectra by interpolating on linearly sampled model grids. Our results show that random or LH sampling out-performs linear sampling in parameter predictability for our higher dimensional models, requiring fewer models in the grid, and thus allowing for more computationally intensive forward models to be used. We also find that using a traditional retrieval with interpolation on a linear grid can produce biased posterior distributions, especially for parameters with non-linear effects on the spectrum. In particular, we advise caution when performing linear interpolation on the C/O ratio, cloud properties, and metallicity. Finally, we find that the information content analysis of the grid from Goyal et al. 2019 is able to highlight key areas of the spectra where the presence or absence of certain molecules can be detected, providing good indicators for parameters such as temperature and C/O ratio.

14 pages, 8 figures. Accepted for publication in ApJ

References in corpus (17)

• Rayleigh scattering in the transit spectrum of HD 189733b
• Fingering convection and cloudless models for cool brown dwarf atmospheres
• Retrieving scattering clouds and disequilibrium chemistry in the atmosphere of HR 8799e
• Observing transiting planets with JWST -- Prime targets and their synthetic spectral observations
• The Effects of Consistent Chemical Kinetics Calculations on the Pressure-Temperature Profiles and Emission Spectra of Hot Jupiters
• HELIOS-K 2.0 Opacity Calculator and Open-source Opacity Database for Exoplanetary Atmospheres
• The Atmospheric Circulation of Ultra-hot Jupiters
• The Impact of Non-Uniform Thermal Structure on the Interpretation of Exoplanet Emission Spectra
• Exo-Transmit: An Open-Source Code for Calculating Transmission Spectra for Exoplanet Atmospheres of Varied Composition
• A Significant Increase in Detection of High-Resolution Emission Spectra Using a Three-Dimensional Atmospheric Model of a Hot Jupiter
• An Ensemble of Bayesian Neural Networks for Exoplanetary Atmospheric Retrieval
• A Library of Self-Consistent Simulated Exoplanet Atmospheres
• An analytical formalism accounting for clouds and other "surfaces" for exoplanet transmission spectroscopy
• Assessment of Supervised Machine Learning for Atmospheric Retrieval of Exoplanets
• Supervised Machine Learning for Intercomparison of Model Grids of Brown Dwarfs: Application to GJ 570D and the Epsilon Indi B Binary System
• On physical interpretations of the reference transit radius of gas-giant exoplanets
• The Sonora Brown Dwarf Atmosphere and Evolution Models I. Model Description and Application to Cloudless Atmospheres in Rainout Chemical Equilibrium

Cited by in corpus (2)

• Grid-Based Atmospheric Retrievals for Reflected-Light Spectra of Exoplanets using PSGnest
• Lessons from Hubble and Spitzer: 1D Self-Consistent Model Grids for 19 Hot Jupiter Emission Spectra