-
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.
-
Large Interferometer For Exoplanets (LIFE): I. Improved exoplanet detection yield estimates for a large mid-infrared space-interferometer mission
Authors:
S. P. Quanz,
M. Ottiger,
E. Fontanet,
J. Kammerer,
F. Menti,
F. Dannert,
A. Gheorghe,
O. Absil,
V. S. Airapetian,
E. Alei,
R. Allart,
D. Angerhausen,
S. Blumenthal,
L. A. Buchhave,
J. Cabrera,
Ó. Carrión-González,
G. Chauvin,
W. C. Danchi,
C. Dandumont,
D. Defrère,
C. Dorn,
D. Ehrenreich,
S. Ertel,
M. Fridlund,
A. García Muñoz
, et al. (46 additional authors not shown)
Abstract:
One of the long-term goals of exoplanet science is the atmospheric characterization of dozens of small exoplanets in order to understand their diversity and search for habitable worlds and potential biosignatures. Achieving this goal requires a space mission of sufficient scale. We seek to quantify the exoplanet detection performance of a space-based mid-infrared nulling interferometer that measur…
▽ More
One of the long-term goals of exoplanet science is the atmospheric characterization of dozens of small exoplanets in order to understand their diversity and search for habitable worlds and potential biosignatures. Achieving this goal requires a space mission of sufficient scale. We seek to quantify the exoplanet detection performance of a space-based mid-infrared nulling interferometer that measures the thermal emission of exoplanets. For this, we have developed an instrument simulator that considers all major astrophysical noise sources and coupled it with Monte Carlo simulations of a synthetic exoplanet population around main-sequence stars within 20 pc. This allows us to quantify the number (and types) of exoplanets that our mission concept could detect over a certain time period. Two different scenarios to distribute the observing time among the stellar targets are discussed and different apertures sizes and wavelength ranges are considered. Within a 2.5-year initial search phase, an interferometer consisting of four 2 m apertures with a total instrument throughput of 5% covering a wavelength range between 4 and 18.5 $μ$m could detect up to ~550 exoplanets with radii between 0.5 and 6 R$_\oplus$ with an integrated SNR$\ge$7. At least ~160 of the detected exoplanets have radii $\le$1.5 R$_\oplus$. Depending on the observing scenario, ~25-45 rocky exoplanets (objects with radii between 0.5 and 1.5 $_{\oplus}$) orbiting within the empirical habitable zone (eHZ) of their host stars are among the detections. With an aperture size of 3.5 m, the total number of detections can increase to up to ~770, including ~60-80 rocky, eHZ planets. With 1 m aperture size, the maximum detection yield is ~315 exoplanets, including $\le$20 rocky, eHZ planets. In terms of predicted detection yield, such a mission can compete with large single-aperture reflected light missions. (abridged)
△ Less
Submitted 20 April, 2022; v1 submitted 19 January, 2021;
originally announced January 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.
-
Helium in the eroding atmosphere of an exoplanet
Authors:
J. J. Spake,
D. K. Sing,
T. M. Evans,
A. Oklopčić,
V. Bourrier,
L. Kreidberg,
B. V. Rackham,
J. Irwin,
D. Ehrenreich,
A. Wyttenbach,
H. R. Wakeford,
Y. Zhou,
K. L. Chubb,
N. Nikolov,
J. M. Goyal,
G. W. Henry,
M. H. Williamson,
S. Blumenthal,
D. R. Anderson,
C. Hellier,
D. Charbonneau,
S. Udry,
N. Madhusudhan
Abstract:
Helium is the second-most abundant element in the Universe after hydrogen and is one of the main constituents of gas-giant planets in our Solar System. Early theoretical models predicted helium to be among the most readily detectable species in the atmospheres of exoplanets, especially in extended and escaping atmospheres. Searches for helium, however, have hitherto been unsuccessful. Here we repo…
▽ More
Helium is the second-most abundant element in the Universe after hydrogen and is one of the main constituents of gas-giant planets in our Solar System. Early theoretical models predicted helium to be among the most readily detectable species in the atmospheres of exoplanets, especially in extended and escaping atmospheres. Searches for helium, however, have hitherto been unsuccessful. Here we report observations of helium on an exoplanet, at a confidence level of 4.5 standard deviations. We measured the near- infrared transmission spectrum of the warm gas giant WASP-107b and identified the narrow absorption feature of excited metastable helium at 10,833 angstroms. The amplitude of the feature, in transit depth, is 0.049 +/- 0.011 per cent in a bandpass of 98 angstroms, which is more than five times greater than what could be caused by nominal stellar chromospheric activity. This large absorption signal suggests that WASP-107b has an extended atmosphere that is eroding at a total rate of 10^10 to 3 x 10^11 grams per second (0.1-4 per cent of its total mass per billion years), and may have a comet-like tail of gas shaped by radiation pressure.
△ Less
Submitted 3 May, 2018;
originally announced May 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.