-
An open-source Bayesian atmospheric radiative transfer (BART) code: III. Initialization, atmospheric profile generator, post-processing routines, and application to exoplanet WASP-43b
Authors:
Jasmina Blecic,
Joseph Harrington,
Patricio E. Cubillos,
M. Oliver Bowman,
Patricio Rojo,
Madison Stemm,
Ryan C. Challener,
Michael D. Himes,
Austin J. Foster,
Ian Dobbs-Dixon,
Andrew S. D. Foster,
Nathaniel B. Lust,
Sarah D. Blumenthal,
Dylan Bruce,
Thomas J. Loredo
Abstract:
This and companion papers by Harrington et al. 2021, submitted and Cubillos et al. 2021, submitted describe an open-source retrieval framework, Bayesian Atmospheric Radiative Transfer (BART), available to the community under the reproducible-research license via https://github.com/exosports/BART . BART is a radiative-transfer code (transit, https://github.com/exosports/transit , Rojo 2009, 2009ASP…
▽ More
This and companion papers by Harrington et al. 2021, submitted and Cubillos et al. 2021, submitted describe an open-source retrieval framework, Bayesian Atmospheric Radiative Transfer (BART), available to the community under the reproducible-research license via https://github.com/exosports/BART . BART is a radiative-transfer code (transit, https://github.com/exosports/transit , Rojo 2009, 2009ASPC..420..321R), initialized by the Thermochemical Equilibrium Abundances (TEA, https://github.com/dzesmin/TEA , Blecic et al. 2016, arXiv:1505.06392) code, and driven through the parameter phase space by a differential-evolution Markov-chain Monte Carlo (MC3, https://github.com/pcubillos/mc3 , Cubillos et al. 2017, arXiv:1610.01336) sampler. In this paper we give a brief description of the framework, and its modules that can be used separately for other scientific purposes; outline the retrieval analysis flow; present the initialization routines, describing in detail the atmospheric profile generator and the temperature and species parameterizations; and specify the post-processing routines and outputs, concentrating on the spectrum band integrator, the best-fit model selection, and the contribution functions. We also present an atmospheric analysis of WASP-43b secondary eclipse data obtained from space- and ground-based observations. We compare our results with the results from the literature, and investigate how the inclusion of additional opacity sources influence the best-fit model.
△ Less
Submitted 26 April, 2021;
originally announced April 2021.
-
An Open-source Bayesian Atmospheric Radiative Transfer (BART) Code: II. The Transit Radiative-transfer Module and Retrieval of HAT-P-11b
Authors:
Patricio E. Cubillos,
Joseph Harrington,
Jasmina Blecic,
Michael D. Himes,
Patricio Rojo,
Thomas J. Loredo,
Nate B. Lust,
Ryan C. Challener,
Austin J. Foster,
Madison M. Stemm,
Andrew S. D. Foster,
Sarah D. Blumenthal
Abstract:
This and companion papers by Harrington et al. and Blecic et al. present the Bayesian Atmospheric Radiative Transfer (BART) code, an open-source, open-development package to characterize extrasolar-planet atmospheres. BART combines a thermochemical equilibrium abundances (TEA), a radiative-transfer (Transit), and a Bayesian statistical (MC3) module to constrain atmospheric temperatures and molecul…
▽ More
This and companion papers by Harrington et al. and Blecic et al. present the Bayesian Atmospheric Radiative Transfer (BART) code, an open-source, open-development package to characterize extrasolar-planet atmospheres. BART combines a thermochemical equilibrium abundances (TEA), a radiative-transfer (Transit), and a Bayesian statistical (MC3) module to constrain atmospheric temperatures and molecular abundances for given spectroscopic observations. Here, we describe the Transit radiative-transfer package, an efficient line-by-line radiative-transfer C code for one-dimensional atmospheres, developed by P. Rojo and further modified by the UCF exoplanet group. This code produces transmission and hemisphere-integrated emission spectra. Transit handles line-by-line opacities from HITRAN, Partridge \& Schwenke ({\water}), Schwenke (TiO), and Plez (VO); and collision-induced absorption from Borysow, HITRAN, and ExoMol. Transit emission-spectra models agree with models from C. Morley (priv. comm.) within a few percent. We applied BART to the {\Spitzer} and {\Hubble} transit observations of the Neptune-sized planet HAT-P-11b. Our results generally agree with those from previous studies, constraining the {\water} abundance and finding an atmosphere enhanced in heavy elements. Different conclusions start to emerge when we make different assumptions from other studies. The BART source code and documentation are available at https://github.com/exosports/BART.
△ Less
Submitted 2 December, 2021; v1 submitted 26 April, 2021;
originally announced April 2021.
-
An Open-Source Bayesian Atmospheric Radiative Transfer (BART) Code: I. Design, Tests, and Application to Exoplanet HD 189733 b
Authors:
Joseph Harrington,
Michael D. Himes,
Patricio E. Cubillos,
Jasmina Blecic,
Patricio M. Rojo,
Ryan C. Challener,
Nate B. Lust,
M. Oliver Bowman,
Sarah D. Blumenthal,
Ian Dobbs-Dixon,
Andrew S. D. Foster,
Austin J. Foster,
M. R. Green,
Thomas J. Loredo,
Kathleen J. McIntyre,
Madison M. Stemm,
David C. Wright
Abstract:
We present the open-source Bayesian Atmospheric Radiative Transfer (BART) retrieval package, which produces estimates and uncertainties for an atmosphere's thermal profile and chemical abundances from observations. Several BART components are also stand-alone packages, including the parallel Multi-Core Markov chain Monte Carlo (MC3), which implements several Bayesian samplers; a line-by-line radia…
▽ More
We present the open-source Bayesian Atmospheric Radiative Transfer (BART) retrieval package, which produces estimates and uncertainties for an atmosphere's thermal profile and chemical abundances from observations. Several BART components are also stand-alone packages, including the parallel Multi-Core Markov chain Monte Carlo (MC3), which implements several Bayesian samplers; a line-by-line radiative-transfer model, transit; a code that calculates Thermochemical Equilibrium Abundances, TEA; and a test suite for verifying radiative-transfer and retrieval codes, BARTTest. The codes are in Python and C. BART and TEA are under a Reproducible Research (RR) license, which requires reviewed-paper authors to publish a compendium of all inputs, codes, and outputs supporting the paper's scientific claims. BART and TEA produce the compendium's content. Otherwise, these codes are under permissive open-source terms, as are MC3 and BARTTest, for any purpose. This paper presents an overview of the code, BARTTest, and an application to eclipse data for exoplanet HD 189733 b. Appendices address RR methodology for accelerating science, a reporting checklist for retrieval papers, the spectral resolution required for synthetic tests, and a derivation of the effective sample size required to estimate any Bayesian posterior distribution to a given precision, which determines how many iterations to run. Paper II, by Cubillos et al., presents the underlying radiative-transfer scheme and an application to transit data for exoplanet HAT-P-11b. Paper III, by Blecic et al., discusses the initialization and post-processing routines, with an application to eclipse data for exoplanet WASP-43b. We invite the community to use and improve BART and its components at http://GitHub.com/ExOSPORTS/BART/.
△ Less
Submitted 26 April, 2021;
originally announced April 2021.
-
An optical transmission spectrum for the ultra-hot Jupiter WASP-121b measured with the Hubble Space Telescope
Authors:
Thomas M. Evans,
David K. Sing,
Jayesh Goyal,
Nikolay Nikolov,
Mark S. Marley,
Kevin Zahnle,
Gregory W. Henry,
Joanna K. Barstow,
Munazza K. Alam,
Jorge Sanz-Forcada,
Tiffany Kataria,
Nikole K. Lewis,
Panayotis Lavvas,
Gilda E. Ballester,
Lotfi Ben-Jaffel,
Sarah D. Blumenthal,
Vincent Bourrier,
Benjamin Drummond,
Antonio Garcia Munoz,
Mercedes Lopez-Morales,
Pascal Tremblin,
David Ehrenreich,
Hannah R. Wakeford,
Lars A. Buchhave,
Alain Lecavelier des Etangs
, et al. (2 additional authors not shown)
Abstract:
We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 m…
▽ More
We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 micron). Wavelength-dependent variations in the opacity rule out a gray cloud deck at a confidence level of 3.8-sigma and may instead be explained by VO spectral bands. We find a cloud-free model assuming chemical equilibrium for a temperature of 1500K and metal enrichment of 10-30x solar matches these data well. Using a free-chemistry retrieval analysis, we estimate a VO abundance of -6.6(-0.3,+0.2) dex. We find no evidence for TiO and place a 3-sigma upper limit of -7.9 dex on its abundance, suggesting TiO may have condensed from the gas phase at the day-night limb. The opacity rises steeply at the shortest wavelengths, increasing by approximately five pressure scale heights from 0.47 to 0.3 micron in wavelength. If this feature is caused by Rayleigh scattering due to uniformly-distributed aerosols, it would imply an unphysically high temperature of 6810+/-1530K. One alternative explanation for the short-wavelength rise is absorption due to SH (mercapto radical), which has been predicted as an important product of non-equilibrium chemistry in hot Jupiter atmospheres. Irrespective of the identity of the NUV absorber, it likely captures a significant amount of incident stellar radiation at low pressures, thus playing a significant role in the overall energy budget, thermal structure, and circulation of the atmosphere.
△ Less
Submitted 25 October, 2018;
originally announced October 2018.
-
A Comparison of Simulated JWST Observations Derived from Equilibrium and Non-Equilibrium Chemistry Models of Giant Exoplanets
Authors:
Sarah D. Blumenthal,
Avi M. Mandell,
Eric Hébrard,
Natasha E. Batalha,
Patricio E. Cubillos,
Sarah Rugheimer,
Hannah R. Wakeford
Abstract:
We aim to see if the difference between equilibrium and disequilibrium chemistry is observable in the atmospheres of transiting planets by the James Webb Space Telescope (JWST). We perform a case study comparing the dayside emission spectra of three planets like HD 189733b, WASP-80b, and GJ436b, in and out of chemical equilibrium at two metallicities each. These three planets were chosen because t…
▽ More
We aim to see if the difference between equilibrium and disequilibrium chemistry is observable in the atmospheres of transiting planets by the James Webb Space Telescope (JWST). We perform a case study comparing the dayside emission spectra of three planets like HD 189733b, WASP-80b, and GJ436b, in and out of chemical equilibrium at two metallicities each. These three planets were chosen because they span a large range of planetary masses and equilibrium temperatures, from hot and Jupiter-sized to warm and Neptune-sized. We link the one-dimensional disequilibrium chemistry model from Venot et al. (2012) in which thermochemical kinetics, vertical transport, and photochemistry are taken into account, to the one-dimensional, pseudo line-by-line radiative transfer model, Pyrat Bay, developed especially for hot Jupiters, and then simulate JWST spectra using PandExo for comparing the effects of temperature, metallicity, and radius. We find the most significant differences from 4 to 5 $μ$m due to disequilibrium from CO and CO$_{2}$ abundances, and also H$_{2}$O for select cases. Our case study shows a certain "sweet spot" of planetary mass, temperature, and metallicity where the difference between equilibrium and disequilibrium is observable. For a planet similar to WASP-80b, JWST's NIRSpec G395M can detect differences due to disequilibrium chemistry with one eclipse event. For a planet similar to GJ 436b, the observability of differences due to disequilibrium chemistry is possible at low metallicity given five eclipse events, but not possible at the higher metallicity.
△ Less
Submitted 4 December, 2017;
originally announced December 2017.