Showing 1–15 of 15 results for author: Challener, R C

Multiple Clues for Dayside Aerosols and Temperature Gradients in WASP-69 b from a Panchromatic JWST Emission Spectrum

Authors: Everett Schlawin, Sagnick Mukherjee, Kazumasa Ohno, Taylor Bell, Thomas G. Beatty, Thomas P. Greene, Michael Line, Ryan C. Challener, Vivien Parmentier, Jonathan J. Fortney, Emily Rauscher, Lindsey Wiser, Luis Welbanks, Matthew Murphy, Isaac Edelman, Natasha Batalha, Sarah E. Moran, Nishil Mehta, Marcia Rieke

Abstract: WASP-69 b is a hot, inflated, Saturn-mass planet 0.26 Mjup with a zero-albedo equilibrium temperature of 963 K. Here, we report the JWST 2 to 12 um emission spectrum of the planet consisting of two eclipses observed with NIRCam grism time series and one eclipse observed with MIRI LRS. The emission spectrum shows absorption features of water vapor, carbon dioxide and carbon monoxide, but no strong… ▽ More WASP-69 b is a hot, inflated, Saturn-mass planet 0.26 Mjup with a zero-albedo equilibrium temperature of 963 K. Here, we report the JWST 2 to 12 um emission spectrum of the planet consisting of two eclipses observed with NIRCam grism time series and one eclipse observed with MIRI LRS. The emission spectrum shows absorption features of water vapor, carbon dioxide and carbon monoxide, but no strong evidence for methane. WASP-69 b's emission spectrum is poorly fit by cloud-free homogeneous models. We find three possible model scenarios for the planet: 1) a Scattering Model that raises the brightness at short wavelengths with a free Geometric Albedo parameter 2) a Cloud Layer model that includes high altitude silicate aerosols to moderate long wavelength emission and 3) a Two-Region model that includes significant dayside inhomogeneity and cloud opacity with two different temperature-pressure profiles. In all cases, aerosols are needed to fit the spectrum of the planet. The Scattering model requires an unexpectedly high Geometric Albedo of 0.64. Our atmospheric retrievals indicate inefficient redistribution of heat and an inhomogeneous dayside distribution, which is tentatively supported by MIRI LRS broadband eclipse maps that show a central concentration of brightness. Our more plausible models (2 and 3) retrieve chemical abundances enriched in heavy elements relative to solar composition by 6x to 14x solar and a C/O ratio of 0.65 to 0.94, whereas the less plausible highly reflective scenario (1) retrieves a slightly lower metallicity and lower C/O ratio. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 41 pages, 19 figures, accepted to the Astronomical Journal

Latitudinal Asymmetry in the Dayside Atmosphere of WASP-43b

Authors: Ryan C. Challener, Zafar Rustamkulov, Elspeth K. H. Lee, Nikole Lewis, David K. Sing, Stephan M. Birkmann, Nicolas Crouzet, Néstor Espinoza, Elena Manjavacas, Natalia Oliveros-Gomez, Jeff A. Valenti, Jingxuan Yang

Abstract: We present two-dimensional near-infrared temperature maps of the canonical hot Jupiter WASP-43b using a phase-curve observation with JWST NIRSpec/G395H. From the white-light planetary transit, we improve constraints on the planet's orbital parameters and measure a planet-to-star radius ratio of $0.15883^{+0.00056}_{-0.00053}$. Using the white-light phase curve, we measure a longitude of maximum br… ▽ More We present two-dimensional near-infrared temperature maps of the canonical hot Jupiter WASP-43b using a phase-curve observation with JWST NIRSpec/G395H. From the white-light planetary transit, we improve constraints on the planet's orbital parameters and measure a planet-to-star radius ratio of $0.15883^{+0.00056}_{-0.00053}$. Using the white-light phase curve, we measure a longitude of maximum brightness of $6.9^{+0^\circ.5}_{-0^\circ.5}$ east of the substellar point and a phase-curve offset of $10.0^{+0^\circ.8}_{-0^\circ.8}$. We also find an $\approx4σ$ detection of a latitudinal hotspot offset of $-13.4^{+3^\circ.2}_{-1^\circ.7}$, the first significant detection of a non-equatorial hotspot in an exoplanet atmosphere. We show that this detection is robust to variations within planetary parameter uncertainties, but only if the transit is used to improve constraints, showing the importance of transit observations to eclipse mapping. Maps retrieved from the NRS1 and NRS2 detectors are similar, with hotspot locations consistent between the two detectors at the $1σ$ level. Our JWST data show brighter (hotter) nightsides and a dimmer (colder) dayside at the shorter wavelengths relative to fits to \textit{Spitzer} 3.6 and 4.5 \microns\ phase curves. Through comparison between our phase curves and a set of general circulation models, we find evidence for clouds on the nightside and atmospheric drag or high metallicity reducing the eastward hotspot offset. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: 15 pages, 8 figures. In review at ApJ Letters

A warm Neptune's methane reveals core mass and vigorous atmospheric mixing

Authors: David K. Sing, Zafar Rustamkulov, Daniel P. Thorngren, Joanna K. Barstow, Pascal Tremblin, Catarina Alves de Oliveira, Tracy L. Beck, Stephan M. Birkmann, Ryan C. Challener, Nicolas Crouzet, Néstor Espinoza, Pierre Ferruit, Giovanna Giardino, Amélie Gressier, Elspeth K. H. Lee, Nikole K. Lewis, Roberto Maiolino, Elena Manjavacas, Bernard J. Rauscher, Marco Sirianni, Jeff A. Valenti

Abstract: Observations of transiting gas giant exoplanets have revealed a pervasive depletion of methane, which has only recently been identified atmospherically. The depletion is thought to be maintained by disequilibrium processes such as photochemistry or mixing from a hotter interior. However, the interiors are largely unconstrained along with the vertical mixing strength and only upper limits on the CH… ▽ More Observations of transiting gas giant exoplanets have revealed a pervasive depletion of methane, which has only recently been identified atmospherically. The depletion is thought to be maintained by disequilibrium processes such as photochemistry or mixing from a hotter interior. However, the interiors are largely unconstrained along with the vertical mixing strength and only upper limits on the CH$_4$ depletion have been available. The warm Neptune WASP-107 b stands out among exoplanets with an unusually low density, reported low core mass, and temperatures amenable to CH$_4$ though previous observations have yet to find the molecule. Here we present a JWST NIRSpec transmission spectrum of WASP-107 b which shows features from both SO$_2$ and CH$_4$ along with H$_2$O, CO$_2$, and CO. We detect methane with 4.2$σ$ significance at an abundance of 1.0$\pm$0.5 ppm, which is depleted by 3 orders of magnitude relative to equilibrium expectations. Our results are highly constraining for the atmosphere and interior, which indicate the envelope has a super-solar metallicity of 43$\pm$8$\times$ solar, a hot interior with an intrinsic temperature of T$_{\rm int}$=460$\pm$40 K, and vigorous vertical mixing which depletes CH4 with a diffusion coefficient of Kzz = 10$^{11.6\pm0.1}$ cm$^2$/s. Photochemistry has a negligible effect on the CH$_4$ abundance, but is needed to account for the SO$_2$. We infer a core mass of 11.5$_{-3.6}^{+3.0}$ M$_{\odot}$, which is much higher than previous upper limits, releasing a tension with core-accretion models. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: Published in Nature at this URL, https://www.nature.com/articles/s41586-024-07395-z . This is the authors version of the manuscript, 20 pages including Methods. Data from Figs. 1, 2, and 3 are available at this URL, https://doi.org/10.5281/zenodo.10891400

Two-Dimensional Eclipse Mapping of the Hot Jupiter WASP-43b with JWST MIRI/LRS

Authors: Mark Hammond, Taylor J. Bell, Ryan C. Challener, Neil T. Lewis, Megan Weiner Mansfield, Isaac Malsky, Emily Rauscher, Jacob L. Bean, Ludmila Carone, João M. Mendonça, Lucas Teinturier, Xianyu Tan, Nicolas Crouzet, Laura Kreidberg, Giuseppe Morello, Vivien Parmentier, Jasmina Blecic, Jean-Michel Désert, Christiane Helling, Pierre-Olivier Lagage, Karan Molaverdikhani, Matthew C. Nixon, Benjamin V. Rackham, Jingxuan Yang

Abstract: We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order… ▽ More We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order $\ell_{\rm max}=2$, alongside the planetary, orbital, stellar, and systematic parameters. This yields a map with a meridionally-averaged eastward hot-spot shift of $(7.75 \pm 0.36)^{\circ}$, with no significant degeneracy between the map and the additional parameters. We show the latitudinal and longitudinal contributions of the day-side emission structure to the eclipse shape, finding a latitudinal signal of $\sim$200 ppm and a longitudinal signal of $\sim$250 ppm. To investigate the sensitivity of the map to the method, we fix the non-mapping parameters and derive an "eigenmap" fitted with an optimised number of orthogonal phase curves, which yields a similar map to the $\ell_{\rm max}=2$ map. We also fit a map up to $\ell_{\rm max}=3$, which shows a smaller hot-spot shift, with a larger uncertainty. These maps are similar to those produced by atmospheric simulations. We conclude that there is a significant mapping signal which constrains the spherical harmonic components of our model up to $\ell_{\rm max}=2$. Alternative mapping models may derive different structures with smaller-scale features; we suggest that further observations of WASP-43b and other planets will drive the development of more robust methods and more accurate maps. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Accepted for publication in The Astronomical Journal

Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b

Authors: Taylor J. Bell, Nicolas Crouzet, Patricio E. Cubillos, Laura Kreidberg, Anjali A. A. Piette, Michael T. Roman, Joanna K. Barstow, Jasmina Blecic, Ludmila Carone, Louis-Philippe Coulombe, Elsa Ducrot, Mark Hammond, João M. Mendonça, Julianne I. Moses, Vivien Parmentier, Kevin B. Stevenson, Lucas Teinturier, Michael Zhang, Natalie M. Batalha, Jacob L. Bean, Björn Benneke, Benjamin Charnay, Katy L. Chubb, Brice-Olivier Demory, Peter Gao , et al. (58 additional authors not shown)

Abstract: Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5… ▽ More Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $μ$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$σ$ upper limit of 1-6 parts per million, depending on model assumptions). △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 61 pages, 13 figures, 4 tables. This preprint has been submitted to and accepted in principle for publication in Nature Astronomy without significant changes

Bringing 2D Eclipse Mapping out of the Shadows with Leave-one-out Cross-validation

Authors: Ryan C. Challener, Luis Welbanks, Peter McGill

Abstract: Eclipse mapping is a technique for inferring 2D brightness maps of transiting exoplanets from the shape of an eclipse light curve. With JWST's unmatched precision, eclipse mapping is now possible for a large number of exoplanets. However, eclipse mapping has only been applied to two planets and the nuances of fitting eclipse maps are not yet fully understood. Here, we use Leave-one-out Cross- Vali… ▽ More Eclipse mapping is a technique for inferring 2D brightness maps of transiting exoplanets from the shape of an eclipse light curve. With JWST's unmatched precision, eclipse mapping is now possible for a large number of exoplanets. However, eclipse mapping has only been applied to two planets and the nuances of fitting eclipse maps are not yet fully understood. Here, we use Leave-one-out Cross- Validation (LOO-CV) to investigate eclipse mapping, with application to a JWST NIRISS/SOSS observation of the ultra-hot Jupiter WASP-18b. LOO-CV is a technique that provides insight into the out-of-sample predictive power of models on a data-point-by-data-point basis. We show that constraints on planetary brightness patterns behave as expected, with large-scale variations driven by the phase-curve variation in the light curve and smaller-scale structures constrained by the eclipse ingress and egress. For WASP-18b we show that the need for higher model complexity (smaller-scale features) is driven exclusively by the shape of the eclipse ingress and egress. We use LOO-CV to investigate the relationship between planetary brightness map components when mapping under a positive-flux constraint to better understand the need for complex models. Finally, we use LOO-CV to understand the degeneracy between the competing ``hotspot'' and ``plateau'' brightness map models of WASP-18b, showing that the plateau model is driven by the ingress shape and the hotspot model is driven by the egress shape, but preference for neither model is due to outliers or unmodeled signals. Based on this analysis, we make recommendations for the use of LOO-CV in future eclipse-mapping studies. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 12 pages, 8 figures. Submitted to AJ. Note that browsers can have difficulty displaying the data-heavy figures

The Eclipse Mapping Null Space: Comparing Theoretical Predictions with Observed Maps

Authors: Ryan C. Challener, Emily Rauscher

Abstract: High-precision exoplanet eclipse light curves, like those possible with JWST, enable flux and temperature mapping of exoplanet atmospheres. These eclipse maps will have unprecedented precision, providing an opportunity to constrain current theoretical predictions of exoplanet atmospheres. However, eclipse mapping has unavoidable mathematical limitations because many map patterns are unobservable.… ▽ More High-precision exoplanet eclipse light curves, like those possible with JWST, enable flux and temperature mapping of exoplanet atmospheres. These eclipse maps will have unprecedented precision, providing an opportunity to constrain current theoretical predictions of exoplanet atmospheres. However, eclipse mapping has unavoidable mathematical limitations because many map patterns are unobservable. This ``null space'' has implications for making comparisons between predictions from general circulation models (GCMs) and the observed planet maps, and, thus, affects our understanding of the physical processes driving the observed maps. We describe the eclipse-mapping null space and show how GCM forward models can be transformed to their observable modes for more appropriate comparison with retrieved eclipse maps, demonstrated with applications to synthetic data of an ultra-hot Jupiter and a cloudy warm Jupiter under JWST-best-case- and extreme-precision observing scenarios. We show that the effects of the null space can be mitigated and manipulated through observational design, and JWST exposure times are short enough to not increase the size of the null space. Furthermore, we show the mathematical connection between the null space and the ``eigenmapping'' method, demonstrating how eigenmaps can be used to understand the null space in a model-independent way. We leverage this connection to incorporate null-space uncertainties in retrieved maps, which increases the uncertainties to now encompass the ground truth for synthetic data. The comparisons between observed maps and forward models that are enabled by this work, and the improved eclipse-mapping uncertainties, will be critical to our interpretation of multidimensional aspects of exoplanets in the JWST era. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 20 pages, 12 figures. Accepted for publication in The Astronomical Journal. Note that PDF readers may blur figures 1 and 3, which can be fixed by zooming in

ThERESA: Three-Dimensional Eclipse Mapping with Application to Synthetic JWST Data

Authors: Ryan C. Challener, Emily Rauscher

Abstract: Spectroscopic eclipse observations, like those possible with the James Webb Space Telescope, should enable 3D mapping of exoplanet daysides. However, fully-flexible 3D planet models are overly complex for the data and computationally infeasible for data-fitting purposes. Here, we present ThERESA, a method to retrieve the 3D thermal structure of an exoplanet from eclipse observations by first retri… ▽ More Spectroscopic eclipse observations, like those possible with the James Webb Space Telescope, should enable 3D mapping of exoplanet daysides. However, fully-flexible 3D planet models are overly complex for the data and computationally infeasible for data-fitting purposes. Here, we present ThERESA, a method to retrieve the 3D thermal structure of an exoplanet from eclipse observations by first retrieving 2D thermal maps at each wavelength and then placing them vertically in the atmosphere. This approach allows the 3D model to include complex thermal structures with a manageable number of parameters, hastening fit convergence and limiting overfitting. An analysis runs in a matter of days. We enforce consistency of the 3D model by comparing vertical placement of the 2D maps with their corresponding contribution functions. To test this approach, we generated a synthetic JWST NIRISS-like observation of a single hot-Jupiter eclipse using a global circulation model of WASP-76b and retrieved its 3D thermal structure. We find that a model which places the 2D maps at different depths depending on latitude and longitude is preferred over a model with a single pressure for each 2D map, indicating that ThERESA is able to retrieve 3D atmospheric structure from JWST observations. We successfully recover the temperatures of the planet's dayside, the eastward shift of its hotspot, and the thermal inversion. ThERESA is open-source and publicly available as a tool for the community. △ Less

Submitted 4 January, 2022; v1 submitted 15 October, 2021; originally announced October 2021.

Comments: 21 pages, 15 figures. Accepted for publication in The Astronomical Journal

Spitzer Dayside Emission of WASP-34b

Authors: Ryan C. Challener, Joseph Harrington, Patricio E. Cubillos, Jasmina Blecic, Barry Smalley

Abstract: We analyzed two eclipse observations of the low-density transiting, likely grazing, exoplanet WASP-34b with the Spitzer Space Telescope's InfraRed Array Camera (IRAC) using two techniques to correct for intrapixel sensitivity variation: Pixel-Level Decorrelation (PLD) and BiLinearly Interpolated Subpixel Sensitivity (BLISS). When jointly fitting both light curves, timing results are consistent wit… ▽ More We analyzed two eclipse observations of the low-density transiting, likely grazing, exoplanet WASP-34b with the Spitzer Space Telescope's InfraRed Array Camera (IRAC) using two techniques to correct for intrapixel sensitivity variation: Pixel-Level Decorrelation (PLD) and BiLinearly Interpolated Subpixel Sensitivity (BLISS). When jointly fitting both light curves, timing results are consistent within 0.7$σ$ between the two models and eclipse depths are consistent within 1.1$σ$, where the difference is due to photometry methods, not the models themselves. By combining published radial velocity data, amateur and professional transit observations, and our eclipse timings, we improved upon measurements of orbital parameters and found an eccentricity consistent with zero (0.0). Atmospheric retrieval, using our Bayesian Atmospheric Radiative Transfer code (BART), shows that the planetary spectrum most resembles a blackbody, with no constraint on molecular abundances or vertical temperature variation. WASP-34b is redder than other warm Jupiters with a similar temperature, hinting at unique chemistry, although further observations are necessary to confirm this. △ Less

Submitted 9 August, 2021; originally announced August 2021.

Comments: 12 pages, 6 figures. Accepted for publication at the Planetary Science Journal

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.

Comments: 16 pages, 14 figures

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.

Comments: Accepted for publication at PSJ

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.

Comments: 38 pages, 17 figures, submitted to The Planetary Science Journal

Identification and Mitigation of a Vibrational Telescope Systematic with Application to Spitzer

Authors: Ryan C. Challener, Joseph Harrington, James Jenkins, Nicolás T. Kurtovic, Ricardo Ramirez, Kathleen J. McIntyre, Michael D. Himes, Eloy Rodríguez, Guillem Anglada-Escudé, Stefan Dreizler, Aviv Ofir, Pablo A. Peña Rojas, Ignasi Ribas, Patricio Rojo, David Kipping, R. Paul Butler, Pedro J. Amado, Cristina Rodríguez-López, Enric Palle, Felipe Murgas

Abstract: We observed Proxima Centauri with the Spitzer Space Telescope InfraRed Array Camera (IRAC) five times in 2016 and 2017 to search for transits of Proxima Centauri b. Following standard analysis procedures, we found three asymmetric, transit-like events that are now understood to be vibrational systematics. This systematic is correlated with the width of the point-response function (PRF), which we m… ▽ More We observed Proxima Centauri with the Spitzer Space Telescope InfraRed Array Camera (IRAC) five times in 2016 and 2017 to search for transits of Proxima Centauri b. Following standard analysis procedures, we found three asymmetric, transit-like events that are now understood to be vibrational systematics. This systematic is correlated with the width of the point-response function (PRF), which we measure with rotated and non-rotated Gaussian fits with respect to the detecor array. We show that the systematic can be removed with a novel application of an adaptive elliptical-aperture photometry technique, and compare the performance of this technique with fixed and variable circular-aperture photometry, using both BiLinearly Interpolated Subpixel Sensitivity (BLISS) maps and non-binned Pixel-Level Decorrelation (PLD). With BLISS maps, elliptical photometry results in a lower standard deviation of normalized residuals, and reduced or similar correlated noise when compared to circular apertures. PLD prefers variable, circular apertures, but generally results in more correlated noise than BLISS. This vibrational effect is likely present in other telescopes and Spitzer observations, where correction could improve results. Our elliptical apertures can be applied to any photometry observations, and may be even more effective when applied to more circular PRFs than Spitzer's. △ Less

Submitted 10 November, 2020; originally announced November 2020.

Comments: 18 pages, 12 figures, accepted for publication in PSJ

Proxima Centauri b is not a transiting exoplanet

Authors: James S. Jenkins, Joseph Harrington, Ryan C. Challener, Nicolás T. Kurtovic, Ricardo Ramirez, Jose Peña, Kathleen J. McIntyre, Michael D. Himes, Eloy Rodríguez, Guillem Anglada-Escudé, Stefan Dreizler, Aviv Ofir, Pablo A. Peña Rojas, Ignasi Ribas, Patricio Rojo, David Kipping, R. Paul Butler, Pedro J. Amado, Cristina Rodríguez-López, Eliza M. -R. Kempton, Enric Palle, Felipe Murgas

Abstract: We report Spitzer Space Telescope observations during predicted transits of the exoplanet Proxima Centauri b. As the nearest terrestrial habitable-zone planet we will ever discover, any potential transit of Proxima b would place strong constraints on its radius, bulk density, and atmosphere. Subsequent transmission spectroscopy and secondary-eclipse measurements could then probe the atmospheric ch… ▽ More We report Spitzer Space Telescope observations during predicted transits of the exoplanet Proxima Centauri b. As the nearest terrestrial habitable-zone planet we will ever discover, any potential transit of Proxima b would place strong constraints on its radius, bulk density, and atmosphere. Subsequent transmission spectroscopy and secondary-eclipse measurements could then probe the atmospheric chemistry, physical processes, and orbit, including a search for biosignatures. However, our photometric results rule out planetary transits at the 200~ppm level at 4.5$~μm$, yielding a 3$σ$ upper radius limit of 0.4~$R_\rm{\oplus}$ (Earth radii). Previous claims of possible transits from optical ground- and space-based photometry were likely correlated noise in the data from Proxima Centauri's frequent flaring. Follow-up observations should focus on planetary radio emission, phase curves, and direct imaging. Our study indicates dramatically reduced stellar activity at near-to-mid infrared wavelengths, compared to the optical. Proxima b is an ideal target for space-based infrared telescopes, if their instruments can be configured to handle Proxima's brightness. △ Less

Submitted 3 May, 2019; originally announced May 2019.

Comments: 8 pages, 3 figures, 2 tables, accepted for publication in MNRAS

The Atmosphere and Interior Structure of HAT-P-13b from Spitzer Secondary Eclipses

Authors: Ryan A. Hardy, Joseph Harrington, Matthew R. Hardin, Nikku Madhusudhan, Thomas J. Loredo, Ryan C. Challener, Andrew S. D. Foster, Patricio E. Cubillos, Jasmina Blecic

Abstract: We present {\em Spitzer} secondary-eclipse observations of the hot Jupiter HAT-P-13 b in the 3.6 {\micron} and 4.5 {\micron} bands. HAT-P-13 b inhabits a two-planet system with a configuration that enables constraints on the planet's second Love number, \math{k\sb{2}}, from precise eccentricity measurements, which in turn constrains models of the planet's interior structure. We exploit the direct… ▽ More We present {\em Spitzer} secondary-eclipse observations of the hot Jupiter HAT-P-13 b in the 3.6 {\micron} and 4.5 {\micron} bands. HAT-P-13 b inhabits a two-planet system with a configuration that enables constraints on the planet's second Love number, \math{k\sb{2}}, from precise eccentricity measurements, which in turn constrains models of the planet's interior structure. We exploit the direct measurements of \math{e \cos ω} from our secondary-eclipse data and combine them with previously published radial velocity data to generate a refined model of the planet's orbit and thus an improved estimate on the possible interval for \math{k\sb{2}}. We report eclipse phases of \math{0.49154 \pm 0.00080} and \math{0.49711 \pm 0.00083} and corresponding \math{e \cos ω} estimates of \math{-0.0136 \pm 0.0013} and \math{-0.0048 \pm 0.0013}. Under the assumptions of previous work, our estimate of \math{k\sb{2}} of 0.81 {\pm} 0.10 is consistent with the lower extremes of possible core masses found by previous models, including models with no solid core. This anomalous result challenges both interior models and the dynamical assumptions that enable them, including the essential assumption of apsidal alignment. We also report eclipse depths of 0.081\% {\pm} 0.008\% in the 3.6 {\micron} channel and 0.088 \% {\pm} 0.028 \% in the 4.5 {\micron} channel. These photometric results are non-uniquely consistent with solar-abundance composition without any thermal inversion. △ Less

Submitted 3 January, 2017; originally announced January 2017.

Comments: Accepted for publication by The Astrophysical Journal December 12, 2016