Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn Open-Source Bayesian Atmospheric Radiative Transfer (BART) Code: I. Design, Tests, and Application to Exoplanet HD 189733 b
371
0
0.0
(
0
)
Authors
Joseph Harrington
Ask ChatGPT about the research
No Arabic abstract
We present the open-source Bayesian Atmospheric Radiative Transfer (BART) retrieval package, which produces estimates and uncertainties for an atmospheres 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 papers scientific claims. BART and TEA produce the compendiums 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/.
rate research
Read 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.
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.
HYPERION is a new three-dimensional dust continuum Monte-Carlo radiative transfer code that is designed to be as generic as possible, allowing radiative transfer to be computed through a variety of three-dimensional grids. The main part of the code is problem-independent, and only requires an arbitrary three-dimensional density structure, dust properties, the position and properties of the illuminating sources, and parameters controlling the running and output of the code. HYPERION is parallelized, and is shown to scale well to thousands of processes. Two common benchmark models for protoplanetary disks were computed, and the results are found to be in excellent agreement with those from other codes. Finally, to demonstrate the capabilities of the code, dust temperatures, SEDs, and synthetic multi-wavelength images were computed for a dynamical simulation of a low-mass star formation region. HYPERION is being actively developed to include new features, and is publicly available (http://www.hyperion-rt.org).
Giant exoplanets orbiting very close to their parent star (hot Jupiters) are subject to tidal forces expected to synchronize their rotational and orbital periods on short timescales (tidal locking). However, spin rotation has never been measured directly for hot Jupiters. Furthermore, their atmospheres can show equatorial super-rotation via strong eastward jet streams, and/or high-altitude winds flowing from the day- to the night-side hemisphere. Planet rotation and atmospheric circulation broaden and distort the planet spectral lines to an extent that is detectable with measurements at high spectral resolution. We observed a transit of the hot Jupiter HD 189733 b around 2.3 {mu}m and at a spectral resolution of R~10$^5$ with CRIRES at the ESO Very Large Telescope. After correcting for the stellar absorption lines and their distortion during transit (the Rossiter-McLaughlin effect), we detect the absorption of carbon monoxide and water vapor in the planet transmission spectrum by cross-correlating with model spectra. The signal is maximized (7.6{sigma}) for a planet rotational velocity of $(3.4^{+1.3}_{-2.1})$ km/s, corresponding to a rotational period of $(1.7^{+2.9}_{-0.4})$ days. This is consistent with the planet orbital period of 2.2 days and therefore with tidal locking. We find that the rotation of HD 189733 b is longer than 1 day (3{sigma}). The data only marginally (1.5{sigma}) prefer models with rotation versus models without rotation. We measure a small day- to night-side wind speed of $(-1.7^{+1.1}_{-1.2})$ km/s. Compared to the recent detection of sodium blue-shifted by (8$pm$2) km/s, this likely implies a strong vertical wind shear between the pressures probed by near-infrared and optical transmission spectroscopy.
We present TransitFit, an open-source Python~3 package designed to fit exoplanetary transit light-curves for transmission spectroscopy studies (Available at https://github.com/joshjchayes/TransitFit and https://github.com/spearnet/TransitFit, with documentation at https://transitfit.readthedocs.io/). TransitFit employs nested sampling to offer efficient and robust multi-epoch, multi-wavelength fitting of transit data obtained from one or more telescopes. TransitFit allows per-telescope detrending to be performed simultaneously with parameter fitting, including the use of user-supplied detrending alogorithms. Host limb darkening can be fitted either independently (uncoupled) for each filter or combined (coupled) using prior conditioning from the PHOENIX stellar atmosphere models. For this TransitFit uses the Limb Darkening Toolkit (LDTk) together with filter profiles, including user-supplied filter profiles. We demonstrate the application of TransitFit in three different contexts. First, we model SPEARNET broadband optical data of the low-density hot-Neptune WASP-127~b. The data were obtained from a globally-distributed network of 0.5m--2.4m telescopes. We find clear improvement in our broadband results using the coupled mode over uncoupled mode, when compared against the higher spectral resolution GTC/OSIRIS transmission spectrum obtained by Chen et al. (2018). Using TransitFit, we fit 26 transit observations by TESS to recover improved ephemerides of the hot-Jupiter WASP-91~b and a transit depth determined to a precision of 170~ppm. Finally, we use TransitFit to conduct an investigation into the contested presence of TTV signatures in WASP-126~b using 126 transits observed by TESS, concluding that there is no statistically significant evidence for such signatures from observations spanning 31 TESS sectors.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
