Showing 1–44 of 44 results for author: Blecic, J

Sulphur dioxide in the mid-infrared transmission spectrum of WASP-39b

Authors: Diana Powell, Adina D. Feinstein, Elspeth K. H. Lee, Michael Zhang, Shang-Min Tsai, Jake Taylor, James Kirk, Taylor Bell, Joanna K. Barstow, Peter Gao, Jacob L. Bean, Jasmina Blecic, Katy L. Chubb, Ian J. M. Crossfield, Sean Jordan, Daniel Kitzmann, Sarah E. Moran, Giuseppe Morello, Julianne I. Moses, Luis Welbanks, Jeehyun Yang, Xi Zhang, Eva-Maria Ahrer, Aaron Bello-Arufe, Jonathan Brande , et al. (48 additional authors not shown)

Abstract: The recent inference of sulphur dioxide (SO$_2$) in the atmosphere of the hot ($\sim$1100 K), Saturn-mass exoplanet WASP-39b from near-infrared JWST observations suggests that photochemistry is a key process in high temperature exoplanet atmospheres. This is due to the low ($<$1 ppb) abundance of SO$_2$ under thermochemical equilibrium, compared to that produced from the photochemistry of H$_2$O a… ▽ More The recent inference of sulphur dioxide (SO$_2$) in the atmosphere of the hot ($\sim$1100 K), Saturn-mass exoplanet WASP-39b from near-infrared JWST observations suggests that photochemistry is a key process in high temperature exoplanet atmospheres. This is due to the low ($<$1 ppb) abundance of SO$_2$ under thermochemical equilibrium, compared to that produced from the photochemistry of H$_2$O and H$_2$S (1-10 ppm). However, the SO$_2$ inference was made from a single, small molecular feature in the transmission spectrum of WASP-39b at 4.05 $μ$m, and therefore the detection of other SO$_2$ absorption bands at different wavelengths is needed to better constrain the SO$_2$ abundance. Here we report the detection of SO$_2$ spectral features at 7.7 and 8.5 $μ$m in the 5-12 $μ$m transmission spectrum of WASP-39b measured by the JWST Mid-Infrared Instrument (MIRI) Low Resolution Spectrometer (LRS). Our observations suggest an abundance of SO$_2$ of 0.5-25 ppm (1$σ$ range), consistent with previous findings. In addition to SO$_2$, we find broad water vapour absorption features, as well as an unexplained decrease in the transit depth at wavelengths longer than 10 $μ$m. Fitting the spectrum with a grid of atmospheric forward models, we derive an atmospheric heavy element content (metallicity) for WASP-39b of $\sim$7.1-8.0 $\times$ solar and demonstrate that photochemistry shapes the spectra of WASP-39b across a broad wavelength range. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Published in Nature

Journal ref: Nature 626, 979-983 (2024)

Simultaneous retrieval of orbital phase resolved JWST/MIRI emission spectra of the hot Jupiter WASP-43b: evidence of water, ammonia and carbon monoxide

Authors: Jingxuan Yang, Mark Hammond, Anjali A. A. Piette, Jasmina Blecic, Taylor J. Bell, Patrick G. J. Irwin, Vivien Parmentier, Shang-Min Tsai, Joanna K. Barstow, Nicolas Crouzet, Laura Kreidberg, João M. Mendonça, Jake Taylor, Robin Baeyens, Kazumasa Ohno, Lucas Teinturier, Matthew C. Nixon

Abstract: Spectroscopic phase curves of hot Jupiters measure their emission spectra at multiple orbital phases, thus enabling detailed characterisation of their atmospheres. Precise constraints on the atmospheric composition of these exoplanets offer insights into their formation and evolution. We analyse four phase-resolved emission spectra of the hot Jupiter WASP-43b, generated from a phase curve observed… ▽ More Spectroscopic phase curves of hot Jupiters measure their emission spectra at multiple orbital phases, thus enabling detailed characterisation of their atmospheres. Precise constraints on the atmospheric composition of these exoplanets offer insights into their formation and evolution. We analyse four phase-resolved emission spectra of the hot Jupiter WASP-43b, generated from a phase curve observed with the MIRI/LRS onboard the JWST, to retrieve its atmospheric properties. Using a parametric 2D temperature model and assuming a chemically homogeneous atmosphere within the observed pressure region, we simultaneously fit the four spectra to constrain the abundances of atmospheric constituents, thereby yielding more precise constraints than previous work that analysed each spectrum independently. Our analysis reveals statistically significant evidence of NH3 (4$σ$) in a hot Jupiter's emission spectra for the first time, along with evidence of H2O (6.5$σ$), CO (3.1$σ$), and a non-detection of CH4. With our abundance constraints, we tentatively estimate the metallicity of WASP-43b at 0.6-6.5$\times$solar and its C/O ratio at 0.6-0.9. Our findings offer vital insights into the atmospheric conditions and formation history of WASP-43b by simultaneously constraining the abundances of carbon, oxygen, and nitrogen-bearing species. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 17 pages, 14 figures, accepted for publication in MNRAS. Comments welcome!

Debris Disks can Contaminate Mid-Infrared Exoplanet Spectra: Evidence for a Circumstellar Debris Disk around Exoplanet Host WASP-39

Authors: Laura Flagg, Alycia J. Weinberger, Taylor J. Bell, Luis Welbanks, Giuseppe Morello, Diana Powell, Jacob L. Bean, Jasmina Blecic, Nicolas Crouzet, Peter Gao, Julie Inglis, James Kirk, Mercedes Lopez-Morales, Karan Molaverdikhani, Nikolay Nikolov, Apurva V. Oza, Benjamin V. Rackham, Seth Redfield, Shang-Min Tsai, Ray Jayawardhana, Laura Kreidberg, Matthew C. Nixon, Kevin B. Stevenson, Jake D. Turner

Abstract: The signal from a transiting planet can be diluted by astrophysical contamination. In the case of circumstellar debris disks, this contamination could start in the mid-infrared and vary as a function of wavelength, which would then change the observed transmission spectrum for any planet in the system. The MIRI/LRS WASP-39b transmission spectrum shows an unexplained dip starting at $\sim$10 $μ$m t… ▽ More The signal from a transiting planet can be diluted by astrophysical contamination. In the case of circumstellar debris disks, this contamination could start in the mid-infrared and vary as a function of wavelength, which would then change the observed transmission spectrum for any planet in the system. The MIRI/LRS WASP-39b transmission spectrum shows an unexplained dip starting at $\sim$10 $μ$m that could be caused by astrophysical contamination. The spectral energy distribution displays excess flux at similar levels to that which are needed to create the dip in the transmission spectrum. In this article, we show that this dip is consistent with the presence of a bright circumstellar debris disk, at a distance of $>$2 au. We discuss how a circumstellar debris disk like that could affect the atmosphere of WASP-39b. We also show that even faint debris disks can be a source of contamination in MIRI exoplanet spectra. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: accepted to ApJL

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

Detection of Carbon Monoxide in the Atmosphere of WASP-39b Applying Standard Cross-Correlation Techniques to JWST NIRSpec G395H Data

Authors: Emma Esparza-Borges, Mercedes López-Morales, Jéa I. Adams Redai, Enric Pallé, James Kirk, Núria Casasayas-Barris, Natasha E. Batalha, Benjamin V. Rackham, Jacob L. Bean, S. L. Casewell, Leen Decin, Leonardo A. Dos Santos, Antonio García Muñoz, Joseph Harrington, Kevin Heng, Renyu Hu, Luigi Mancini, Karan Molaverdikhani, Giuseppe Morello, Nikolay K. Nikolov, Matthew C. Nixon, Seth Redfield, Kevin B. Stevenson, Hannah R. Wakeford, Munazza K. Alam , et al. (8 additional authors not shown)

Abstract: Carbon monoxide was recently reported in the atmosphere of the hot Jupiter WASP-39b using the NIRSpec PRISM transit observation of this planet, collected as part of the JWST Transiting Exoplanet Community Early Release Science (JTEC ERS) Program. This detection, however, could not be confidently confirmed in the initial analysis of the higher resolution observations with NIRSpec G395H disperser. H… ▽ More Carbon monoxide was recently reported in the atmosphere of the hot Jupiter WASP-39b using the NIRSpec PRISM transit observation of this planet, collected as part of the JWST Transiting Exoplanet Community Early Release Science (JTEC ERS) Program. This detection, however, could not be confidently confirmed in the initial analysis of the higher resolution observations with NIRSpec G395H disperser. Here we confirm the detection of CO in the atmosphere of WASP-39b using the NIRSpec G395H data and cross-correlation techniques. We do this by searching for the CO signal in the unbinned transmission spectrum of the planet between 4.6 and 5.0 $μ$m, where the contribution of CO is expected to be higher than that of other anticipated molecules in the planet's atmosphere. Our search results in a detection of CO with a cross-correlation function (CCF) significance of $6.6 σ$ when using a template with only ${\rm ^{12}C^{16}O}$ lines. The CCF significance of the CO signal increases to $7.5 σ$ when including in the template lines from additional CO isotopologues, with the largest contribution being from ${\rm ^{13}C^{16}O}$. Our results highlight how cross-correlation techniques can be a powerful tool for unveiling the chemical composition of exoplanetary atmospheres from medium-resolution transmission spectra, including the detection of isotopologues. △ Less

Submitted 31 August, 2023; originally announced September 2023.

Comments: Accepted for publication in The Astrophysical Journal Letters

Awesome SOSS: Transmission Spectroscopy of WASP-96b with NIRISS/SOSS

Authors: Michael Radica, Luis Welbanks, Néstor Espinoza, Jake Taylor, Louis-Philippe Coulombe, Adina D. Feinstein, Jayesh Goyal, Nicholas Scarsdale, Loic Albert, Priyanka Baghel, Jacob L. Bean, Jasmina Blecic, David Lafrenière, Ryan J. MacDonald, Maria Zamyatina, Romain Allart, Étienne Artigau, Natasha E. Batalha, Neil James Cook, Nicolas B. Cowan, Lisa Dang, René Doyon, Marylou Fournier-Tondreau, Doug Johnstone, Michael R. Line , et al. (8 additional authors not shown)

Abstract: The future is now - after its long-awaited launch in December 2021, JWST began science operations in July 2022 and is already revolutionizing exoplanet astronomy. The Early Release Observations (ERO) program was designed to provide the first images and spectra from JWST, covering a multitude of science cases and using multiple modes of each on-board instrument. Here, we present transmission spectr… ▽ More The future is now - after its long-awaited launch in December 2021, JWST began science operations in July 2022 and is already revolutionizing exoplanet astronomy. The Early Release Observations (ERO) program was designed to provide the first images and spectra from JWST, covering a multitude of science cases and using multiple modes of each on-board instrument. Here, we present transmission spectroscopy observations of the hot-Saturn WASP-96b with the Single Object Slitless Spectroscopy (SOSS) mode of the Near Infrared Imager and Slitless Spectrograph, observed as part of the ERO program. As the SOSS mode presents some unique data reduction challenges, we provide an in-depth walk-through of the major steps necessary for the reduction of SOSS data: including background subtraction, correction of 1/f noise, and treatment of the trace order overlap. We furthermore offer potential routes to correct for field star contamination, which can occur due to the SOSS mode's slitless nature. By comparing our extracted transmission spectrum with grids of atmosphere models, we find an atmosphere metallicity between 1x and 5x solar, and a solar carbon-to-oxygen ratio. Moreover, our models indicate that no grey cloud deck is required to fit WASP-96b's transmission spectrum, but find evidence for a slope shortward of 0.9$μ$m, which could either be caused by enhanced Rayleigh scattering or the red wing of a pressure-broadened Na feature. Our work demonstrates the unique capabilities of the SOSS mode for exoplanet transmission spectroscopy and presents a step-by-step reduction guide for this new and exciting instrument. △ Less

Submitted 20 June, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: MNRAS, in press. Updated to reflect published version

Awesome SOSS: Atmospheric Characterisation of WASP-96 b using the JWST Early Release Observations

Authors: Jake Taylor, Michael Radica, Luis Welbanks, Ryan J. MacDonald, Jasmina Blecic, Maria Zamyatina, Alexander Roth, Jacob L. Bean, Vivien Parmentier, Louis-Philippe Coulombe, Adina D. Feinstein, Néstor Espinoza, Björn Benneke, David Lafrenière, René Doyon, Eva-Maria Ahrer

Abstract: The newly operational JWST offers the potential to study the atmospheres of distant worlds with precision that has not been achieved before. One of the first exoplanets observed by JWST in the summer of 2022 was WASP-96 b, a hot-Saturn orbiting a G8 star. As part of the Early Release Observations program, one transit of WASP-96 b was observed with NIRISS/SOSS to capture its transmission spectrum f… ▽ More The newly operational JWST offers the potential to study the atmospheres of distant worlds with precision that has not been achieved before. One of the first exoplanets observed by JWST in the summer of 2022 was WASP-96 b, a hot-Saturn orbiting a G8 star. As part of the Early Release Observations program, one transit of WASP-96 b was observed with NIRISS/SOSS to capture its transmission spectrum from 0.6-2.85 microns. In this work, we utilise four retrieval frameworks to report precise and robust measurements of WASP-96 b's atmospheric composition. We constrain the logarithmic volume mixing ratios of multiple chemical species in its atmosphere, including: H$_2$O = $-3.59 ^{+ 0.35 }_{- 0.35 }$, CO$_2$ = $-4.38 ^{+ 0.47 }_{- 0.57 }$ and K = $-8.04 ^{+ 1.22 }_{- 1.71 }$. Notably, our results offer a first abundance constraint on potassium in WASP-96 b's atmosphere, and important inferences on carbon-bearing species such as CO$_2$ and CO. Our short wavelength NIRISS/SOSS data are best explained by the presence of an enhanced Rayleigh scattering slope, despite previous inferences of a clear atmosphere - although we find no evidence for a grey cloud deck. Finally, we explore the data resolution required to appropriately interpret observations using NIRISS/SOSS. We find that our inferences are robust against different binning schemes. That is, from low $R = 125$ to the native resolution of the instrument, the bulk atmospheric properties of the planet are consistent. Our systematic analysis of these exquisite observations demonstrates the power of NIRISS/SOSS to detect and constrain multiple molecular and atomic species in the atmospheres of hot giant planets. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 12 pages, 5 Figures. Accepted for publication in MNRAS. Companion paper to Radica et al., 2023

A broadband thermal emission spectrum of the ultra-hot Jupiter WASP-18b

Authors: Louis-Philippe Coulombe, Björn Benneke, Ryan Challener, Anjali A. A. Piette, Lindsey S. Wiser, Megan Mansfield, Ryan J. MacDonald, Hayley Beltz, Adina D. Feinstein, Michael Radica, Arjun B. Savel, Leonardo A. Dos Santos, Jacob L. Bean, Vivien Parmentier, Ian Wong, Emily Rauscher, Thaddeus D. Komacek, Eliza M. -R. Kempton, Xianyu Tan, Mark Hammond, Neil T. Lewis, Michael R. Line, Elspeth K. H. Lee, Hinna Shivkumar, Ian J. M. Crossfield , et al. (51 additional authors not shown)

Abstract: Close-in giant exoplanets with temperatures greater than 2,000 K (''ultra-hot Jupiters'') have been the subject of extensive efforts to determine their atmospheric properties using thermal emission measurements from the Hubble and Spitzer Space Telescopes. However, previous studies have yielded inconsistent results because the small sizes of the spectral features and the limited information conten… ▽ More Close-in giant exoplanets with temperatures greater than 2,000 K (''ultra-hot Jupiters'') have been the subject of extensive efforts to determine their atmospheric properties using thermal emission measurements from the Hubble and Spitzer Space Telescopes. However, previous studies have yielded inconsistent results because the small sizes of the spectral features and the limited information content of the data resulted in high sensitivity to the varying assumptions made in the treatment of instrument systematics and the atmospheric retrieval analysis. Here we present a dayside thermal emission spectrum of the ultra-hot Jupiter WASP-18b obtained with the NIRISS instrument on JWST. The data span 0.85 to 2.85 $μ$m in wavelength at an average resolving power of 400 and exhibit minimal systematics. The spectrum shows three water emission features (at $>$6$σ$ confidence) and evidence for optical opacity, possibly due to H$^-$, TiO, and VO (combined significance of 3.8$σ$). Models that fit the data require a thermal inversion, molecular dissociation as predicted by chemical equilibrium, a solar heavy element abundance (''metallicity'', M/H = 1.03$_{-0.51}^{+1.11}$ $\times$ solar), and a carbon-to-oxygen (C/O) ratio less than unity. The data also yield a dayside brightness temperature map, which shows a peak in temperature near the sub-stellar point that decreases steeply and symmetrically with longitude toward the terminators. △ Less

Submitted 20 January, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: JWST ERS bright star observations. Uploaded to inform JWST Cycle 2 proposals. Manuscript under review. 50 pages, 14 figures, 2 tables

Early Release Science of the exoplanet WASP-39b with JWST NIRISS

Authors: Adina D. Feinstein, Michael Radica, Luis Welbanks, Catriona Anne Murray, Kazumasa Ohno, Louis-Philippe Coulombe, Néstor Espinoza, Jacob L. Bean, Johanna K. Teske, Björn Benneke, Michael R. Line, Zafar Rustamkulov, Arianna Saba, Angelos Tsiaras, Joanna K. Barstow, Jonathan J. Fortney, Peter Gao, Heather A. Knutson, Ryan J. MacDonald, Thomas Mikal-Evans, Benjamin V. Rackham, Jake Taylor, Vivien Parmentier, Natalie M. Batalha, Zachory K. Berta-Thompson , et al. (64 additional authors not shown)

Abstract: Transmission spectroscopy provides insight into the atmospheric properties and consequently the formation history, physics, and chemistry of transiting exoplanets. However, obtaining precise inferences of atmospheric properties from transmission spectra requires simultaneously measuring the strength and shape of multiple spectral absorption features from a wide range of chemical species. This has… ▽ More Transmission spectroscopy provides insight into the atmospheric properties and consequently the formation history, physics, and chemistry of transiting exoplanets. However, obtaining precise inferences of atmospheric properties from transmission spectra requires simultaneously measuring the strength and shape of multiple spectral absorption features from a wide range of chemical species. This has been challenging given the precision and wavelength coverage of previous observatories. Here, we present the transmission spectrum of the Saturn-mass exoplanet WASP-39b obtained using the SOSS mode of the NIRISS instrument on the JWST. This spectrum spans $0.6 - 2.8 μ$m in wavelength and reveals multiple water absorption bands, the potassium resonance doublet, as well as signatures of clouds. The precision and broad wavelength coverage of NIRISS-SOSS allows us to break model degeneracies between cloud properties and the atmospheric composition of WASP-39b, favoring a heavy element enhancement ("metallicity") of $\sim 10 - 30 \times$ the solar value, a sub-solar carbon-to-oxygen (C/O) ratio, and a solar-to-super-solar potassium-to-oxygen (K/O) ratio. The observations are best explained by wavelength-dependent, non-gray clouds with inhomogeneous coverage of the planet's terminator. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 48 pages, 12 figures, 2 tables. Under review at Nature

Photochemically-produced SO$_2$ in the atmosphere of WASP-39b

Authors: Shang-Min Tsai, Elspeth K. H. Lee, Diana Powell, Peter Gao, Xi Zhang, Julianne Moses, Eric Hébrard, Olivia Venot, Vivien Parmentier, Sean Jordan, Renyu Hu, Munazza K. Alam, Lili Alderson, Natalie M. Batalha, Jacob L. Bean, Björn Benneke, Carver J. Bierson, Ryan P. Brady, Ludmila Carone, Aarynn L. Carter, Katy L. Chubb, Julie Inglis, Jérémy Leconte, Mercedes Lopez-Morales, Yamila Miguel , et al. (60 additional authors not shown)

Abstract: Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $μ$m arising from SO$_2$ in the atmosphere of WA… ▽ More Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $μ$m arising from SO$_2$ in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 M$_J$) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of $\sim$1100 K. The most plausible way of generating SO$_2$ in such an atmosphere is through photochemical processes. Here we show that the SO$_2$ distribution computed by a suite of photochemical models robustly explains the 4.05 $μ$m spectral feature identified by JWST transmission observations with NIRSpec PRISM (2.7$σ$) and G395H (4.5$σ$). SO$_2$ is produced by successive oxidation of sulphur radicals freed when hydrogen sulphide (H$_2$S) is destroyed. The sensitivity of the SO$_2$ feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of $\sim$10$\times$ solar. We further point out that SO$_2$ also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations. △ Less

Submitted 24 March, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 39 pages, 14 figures, accepted to be published in Nature

Early Release Science of the exoplanet WASP-39b with JWST NIRCam

Authors: Eva-Maria Ahrer, Kevin B. Stevenson, Megan Mansfield, Sarah E. Moran, Jonathan Brande, Giuseppe Morello, Catriona A. Murray, Nikolay K. Nikolov, Dominique J. M. Petit dit de la Roche, Everett Schlawin, Peter J. Wheatley, Sebastian Zieba, Natasha E. Batalha, Mario Damiano, Jayesh M Goyal, Monika Lendl, Joshua D. Lothringer, Sagnick Mukherjee, Kazumasa Ohno, Natalie M. Batalha, Matthew P. Battley, Jacob L. Bean, Thomas G. Beatty, Björn Benneke, Zachory K. Berta-Thompson , et al. (74 additional authors not shown)

Abstract: Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength covera… ▽ More Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution, and high precision that, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWST's Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0 - 4.0 $μ$m, exhibit minimal systematics, and reveal well-defined molecular absorption features in the planet's spectrum. Specifically, we detect gaseous H$_2$O in the atmosphere and place an upper limit on the abundance of CH$_4$. The otherwise prominent CO$_2$ feature at 2.8 $μ$m is largely masked by H$_2$O. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100$\times$ solar (i.e., an enrichment of elements heavier than helium relative to the Sun) and a sub-stellar carbon-to-oxygen (C/O) ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation or disequilibrium processes in the upper atmosphere. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 35 pages, 13 figures, 3 tables, Nature, accepted

Early Release Science of the Exoplanet WASP-39b with JWST NIRSpec G395H

Authors: Lili Alderson, Hannah R. Wakeford, Munazza K. Alam, Natasha E. Batalha, Joshua D. Lothringer, Jea Adams Redai, Saugata Barat, Jonathan Brande, Mario Damiano, Tansu Daylan, Néstor Espinoza, Laura Flagg, Jayesh M. Goyal, David Grant, Renyu Hu, Julie Inglis, Elspeth K. H. Lee, Thomas Mikal-Evans, Lakeisha Ramos-Rosado, Pierre-Alexis Roy, Nicole L. Wallack, Natalie M. Batalha, Jacob L. Bean, Björn Benneke, Zachory K. Berta-Thompson , et al. (67 additional authors not shown)

Abstract: Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems. Access to an exoplanet's chemical inventory requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based and high-resolution ground-based facilities. Here we report the m… ▽ More Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems. Access to an exoplanet's chemical inventory requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based and high-resolution ground-based facilities. Here we report the medium-resolution (R$\sim$600) transmission spectrum of an exoplanet atmosphere between 3-5 $μ$m covering multiple absorption features for the Saturn-mass exoplanet WASP-39b, obtained with JWST NIRSpec G395H. Our observations achieve 1.46x photon precision, providing an average transit depth uncertainty of 221 ppm per spectroscopic bin, and present minimal impacts from systematic effects. We detect significant absorption from CO$_2$ (28.5$σ$) and H$_2$O (21.5$σ$), and identify SO$_2$ as the source of absorption at 4.1 $μ$m (4.8$σ$). Best-fit atmospheric models range between 3 and 10x solar metallicity, with sub-solar to solar C/O ratios. These results, including the detection of SO$_2$, underscore the importance of characterising the chemistry in exoplanet atmospheres, and showcase NIRSpec G395H as an excellent mode for time series observations over this critical wavelength range. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 44 pages, 11 figures, 3 tables. Resubmitted after revision to Nature

Early Release Science of the exoplanet WASP-39b with JWST NIRSpec PRISM

Authors: Z. Rustamkulov, D. K. Sing, S. Mukherjee, E. M. May, J. Kirk, E. Schlawin, M. R. Line, C. Piaulet, A. L. Carter, N. E. Batalha, J. M. Goyal, M. López-Morales, J. D. Lothringer, R. J. MacDonald, S. E. Moran, K. B. Stevenson, H. R. Wakeford, N. Espinoza, J. L. Bean, N. M. Batalha, B. Benneke, Z. K. Berta-Thompson, I. J. M. Crossfield, P. Gao, L. Kreidberg , et al. (69 additional authors not shown)

Abstract: Transmission spectroscopy of exoplanets has revealed signatures of water vapor, aerosols, and alkali metals in a few dozen exoplanet atmospheres. However, these previous inferences with the Hubble and Spitzer Space Telescopes were hindered by the observations' relatively narrow wavelength range and spectral resolving power, which precluded the unambiguous identification of other chemical species… ▽ More Transmission spectroscopy of exoplanets has revealed signatures of water vapor, aerosols, and alkali metals in a few dozen exoplanet atmospheres. However, these previous inferences with the Hubble and Spitzer Space Telescopes were hindered by the observations' relatively narrow wavelength range and spectral resolving power, which precluded the unambiguous identification of other chemical species$-$in particular the primary carbon-bearing molecules. Here we report a broad-wavelength 0.5-5.5 $μ$m atmospheric transmission spectrum of WASP-39 b, a 1200 K, roughly Saturn-mass, Jupiter-radius exoplanet, measured with JWST NIRSpec's PRISM mode as part of the JWST Transiting Exoplanet Community Early Release Science Team program. We robustly detect multiple chemical species at high significance, including Na (19$σ$), H$_2$O (33$σ$), CO$_2$ (28$σ$), and CO (7$σ$). The non-detection of CH$_4$, combined with a strong CO$_2$ feature, favours atmospheric models with a super-solar atmospheric metallicity. An unanticipated absorption feature at 4$μ$m is best explained by SO$_2$ (2.7$σ$), which could be a tracer of atmospheric photochemistry. These observations demonstrate JWST's sensitivity to a rich diversity of exoplanet compositions and chemical processes. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 41 pages, 4 main figures, 10 extended data figures, 4 tables. Under review in Nature

Identification of carbon dioxide in an exoplanet atmosphere

Authors: The JWST Transiting Exoplanet Community Early Release Science Team, Eva-Maria Ahrer, Lili Alderson, Natalie M. Batalha, Natasha E. Batalha, Jacob L. Bean, Thomas G. Beatty, Taylor J. Bell, Björn Benneke, Zachory K. Berta-Thompson, Aarynn L. Carter, Ian J. M. Crossfield, Néstor Espinoza, Adina D. Feinstein, Jonathan J. Fortney, Neale P. Gibson, Jayesh M. Goyal, Eliza M. -R. Kempton, James Kirk, Laura Kreidberg, Mercedes López-Morales, Michael R. Line, Joshua D. Lothringer, Sarah E. Moran, Sagnick Mukherjee , et al. (107 additional authors not shown)

Abstract: Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres… ▽ More Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2 but have not yielded definitive detections due to the lack of unambiguous spectroscopic identification. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science Program (ERS). The data used in this study span 3.0 to 5.5 μm in wavelength and show a prominent CO2 absorption feature at 4.3 μm (26σ significance). The overall spectrum is well matched by one-dimensional, 10x solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide, and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 μm that is not reproduced by these models. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

GCM-Motivated Multidimensional Temperature Parametrization Scheme for Phasecurve Retrieval

Authors: Ian Dobbs-Dixon, Jasmina Blecic

Abstract: We present a novel physically motivated, parametrized temperature model for phase-curve retrieval, able to self-consistently assess the variation in thermal structure in multidimensions. To develop this approach, we drew motivation from both full three-dimensional general circulation models and analytic formulations, accounting for the dominant dynamical feature of tidally locked planets, the plan… ▽ More We present a novel physically motivated, parametrized temperature model for phase-curve retrieval, able to self-consistently assess the variation in thermal structure in multidimensions. To develop this approach, we drew motivation from both full three-dimensional general circulation models and analytic formulations, accounting for the dominant dynamical feature of tidally locked planets, the planetary jet. Our formulation shows notable flexibility. It can generate planetary jets of various characteristics and redistribution efficiencies seen in the literature, including both standard eastward and unusual westward offset hotspots, as well as more exotic configurations for potential future observations. In our modeling scheme we utilize a tractable set of parameters efficient enough to enable future Bayesian analysis and, in addition to the resolved temperature structure, we return physical insights not yet derived from retrievals: the amplitude and the phase offset, and the location and the extent of the equatorial jet. △ Less

Submitted 8 September, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 10 pages, 12 figures

Journal ref: The Astrophysical Journal, 929:46 (10pp), 2022 April 10

Probing the atmosphere of WASP-69 b with low- and high-resolution transmission spectroscopy

Authors: S. Khalafinejad, K. Molaverdikhani, J. Blecic, M. Mallonn, L. Nortmann, J. A. Caballero, H. Rahmati, A. Kaminski, S. Sadegi, E. Nagel, L. Carone, P. J. Amado, M. Azzaro, F. F. Bauer, N. Casasayas-Barris, S. Czesla, C. von Essen, L. Fossati, M. Güdel, Th. Henning, M. López-Puertas, M. Lendl, T. Lüftinger, D. Montes, M. Oshagh , et al. (9 additional authors not shown)

Abstract: Consideration of both low- and high-resolution transmission spectroscopy is key for obtaining a comprehensive picture of exoplanet atmospheres. In studies of transmission spectra, the continuum information is well established with low-resolution spectra, while the shapes of individual lines are best constrained with high-resolution observations. In this work, we aim to merge high- with low-resolut… ▽ More Consideration of both low- and high-resolution transmission spectroscopy is key for obtaining a comprehensive picture of exoplanet atmospheres. In studies of transmission spectra, the continuum information is well established with low-resolution spectra, while the shapes of individual lines are best constrained with high-resolution observations. In this work, we aim to merge high- with low-resolution transmission spectroscopy. We present the analysis of three primary transits of WASP-69b in the VIS channel of the CARMENES instrument and perform a combined low- and high-resolution analysis using additional data from HARPS-N, OSIRIS/GTC, and WFC3/HST already available in the literature. During the first CARMENES observing night, we detected the planet Na D$_{2}$ and D$_{1}$ lines at $\sim 7σ$ and $\sim 3σ$ significance levels, respectively. We measured a D$_{2}$/D$_{1}$ intensity ratio of 2.5$\pm$0.7, which is in agreement with previous HARPS-N observations. Our modelling of WFC3 and OSIRIS data suggests strong Rayleigh scattering, solar to super-solar water abundance, and a highly muted Na feature in the atmosphere of this planet, in agreement with previous investigations of this target. We use the continuum information retrieved from the low-resolution spectroscopy as a prior to break the degeneracy between the Na abundance, reference pressure, and thermosphere temperature for the high-resolution spectroscopic analysis. We fit the Na D$_{1}$ and D$_{2}$ lines individually and find that the posterior distributions of the model parameters agree with each other within 1$σ$. Our results suggest that local thermodynamic equilibrium processes can explain the observed D$_{2}$/D$_{1}$ ratio because the presence of haze opacity mutes the absorption features. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 20 pages, 18 figures, accepted for publication in A&A

Journal ref: A&A 656, A142 (2021)

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

The Pyrat Bay Framework for Exoplanet Atmospheric Modeling: A Population Study of Hubble/WFC3 Transmission Spectra

Authors: Patricio E. Cubillos, Jasmina Blecic

Abstract: We present the open-source pyratbay framework for exoplanet atmospheric modeling, spectral synthesis, and Bayesian retrieval. The modular design of the code allows the users to generate atmospheric 1D parametric models of the temperature, abundances (in thermochemical equilibrium or constant-with-altitude), and altitude profiles in hydrostatic equilibrium; sample ExoMol and HITRAN line-by-line cro… ▽ More We present the open-source pyratbay framework for exoplanet atmospheric modeling, spectral synthesis, and Bayesian retrieval. The modular design of the code allows the users to generate atmospheric 1D parametric models of the temperature, abundances (in thermochemical equilibrium or constant-with-altitude), and altitude profiles in hydrostatic equilibrium; sample ExoMol and HITRAN line-by-line cross sections with custom resolving power and line-wing cutoff values; compute emission or transmission spectra considering cross sections from molecular line transitions, collision-induced absorption, Rayleigh scattering, gray clouds, and alkali resonance lines; and perform Markov chain Monte Carlo atmospheric retrievals for a given transit or eclipse dataset. We benchmarked the pyratbay framework by reproducing line-by-line cross-section sampling of ExoMol cross sections, producing transmission and emission spectra consistent with petitRADTRANS models, accurately retrieving the atmospheric properties of simulated transmission and emission observations generated with TauREx models, and closely reproducing Aura retrieval analyses of the space-based transmission spectrum of HD 209458b. Finally, we present a retrieval analysis of a population of transiting exoplanets, focusing on those observed in transmission with the HST WFC3/G141 grism. We found that this instrument alone can confidently identify when a dataset shows H2O-absorption features; however, it cannot distinguish whether a muted H2O feature is caused by clouds, high atmospheric metallicity, or low H2O abundance. Our results are consistent with previous retrieval analyses. The pyratbay code is available at PyPI (pip install pyratbay) and conda. The code is heavily documented (https://pyratbay.readthedocs.io) and tested to provide maximum accessibility to the community and long-term development stability. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: Accepted for publication at MNRAS

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

Global Chemistry and Thermal Structure Models for the Hot Jupiter WASP-43b and Predictions for JWST

Authors: Olivia Venot, Vivien Parmentier, Jasmina Blecic, Patricio E. Cubillos, Ingo P. Waldmann, Quentin Changeat, Julianne I. Moses, Pascal Tremblin, Nicolas Crouzet, Peter Gao, Diana Powell, Pierre-Olivier Lagage, Ian Dobbs-Dixon, Maria E. Steinrueck, Laura Kreidberg, Natalie Batalha, Jacob L. Bean, Kevin B. Stevenson, Sarah Casewell, Ludmila Carone

Abstract: The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, in… ▽ More The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $μ$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $μ$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition. △ Less

Submitted 14 January, 2020; originally announced January 2020.

Comments: 26 pages, 16 figures, accepted in The Astrophysical Journal

The Need for Laboratory Measurements and Ab Initio Studies to Aid Understanding of Exoplanetary Atmospheres

Authors: Jonathan J. Fortney, Tyler D. Robinson, Shawn Domagal-Goldman, Anthony D. Del Genio, Iouli E. Gordon, Ehsan Gharib-Nezhad, Nikole Lewis, Clara Sousa-Silva, Vladimir Airapetian, Brian Drouin, Robert J. Hargreaves, Xinchuan Huang, Tijs Karman, Ramses M. Ramirez, Gregory B. Rieker, Jonathan Tennyson, Robin Wordsworth, Sergei N Yurchenko, Alexandria V Johnson, Timothy J. Lee, Chuanfei Dong, Stephen Kane, Mercedes Lopez-Morales, Thomas Fauchez, Timothy Lee , et al. (63 additional authors not shown)

Abstract: We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as c… ▽ More We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: Submitted as an Astro2020 Science White Paper

Imaging Cool Giant Planets in Reflected Light: Science Investigations and Synergy with Habitable Planets

Authors: Mark Marley, Nikole Lewis, Giada Arney, Vanessa Bailey, Natasha Batalha, Charles Beichman, Björn Benneke, Jasmina Blecic, Kerri Cahoy, Jeffrey Chilcote, Shawn Domagal-Goldman, Courtney Dressing, Michael Fitzgerald, Jonathan Fortney, Richard Freedman, Dawn Gelino, John Gizis, Olivier Guyon, Thomas Greene, Heidi Hammel, Yasuhiro Hasegawa, Nemanja Jovanovic, Quinn Konopacky, Ravi Kopparapu, Michael Liu , et al. (16 additional authors not shown)

Abstract: Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the… ▽ More Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the decades to follow. △ Less

Submitted 26 March, 2019; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro 2020 Decadal Survey on Astronomy and Astrophysics. Replace version to fix typo in co-signer name and add figure credits

Toward More Reliable Analytic Thermochemical-equilibrium Abundances

Authors: Patricio Cubillos, Jasmina Blecic, Ian Dobbs-Dixon

Abstract: Heng & Tsai (2016) developed an analytic framework to obtain thermochemical-equilibrium abundances for H$_{2}$O, CO, CO$_2$, CH$_4$, C$_{2}$H$_{2}$, C$_{2}$H$_{4}$, HCN, NH$_3$, and N$_2$ for a system with known temperature, pressure, and elemental abundances (hydrogen, carbon, nitrogen, and oxygen). However, the implementation of their approach can become numerically unstable under certain circum… ▽ More Heng & Tsai (2016) developed an analytic framework to obtain thermochemical-equilibrium abundances for H$_{2}$O, CO, CO$_2$, CH$_4$, C$_{2}$H$_{2}$, C$_{2}$H$_{4}$, HCN, NH$_3$, and N$_2$ for a system with known temperature, pressure, and elemental abundances (hydrogen, carbon, nitrogen, and oxygen). However, the implementation of their approach can become numerically unstable under certain circumstances, leading to inaccurate solutions (e.g., ${\rm C/O} \ge 1$ atmospheres at low pressures). Building up on their approach, we identified the conditions that prompt inaccurate solutions, and developed a new framework to avoid them, providing a reliable implementation for arbitrary values of temperature (200 to $\sim$2000 K), pressure ($10^{-8}$ to $10^{3}$ bar), and CNO abundances ($10^{-3}$ to $\sim 10^{2}\times$ solar elemental abundances), for hydrogen-dominated atmospheres. The accuracy our analytic framework is better than 10% for the more abundant species that have mixing fractions larger than $10^{-10}$, whereas the accuracy is better than 50% for the less abundant species. Additionally, we added the equilibrium-abundance calculation of atomic and molecular hydrogen into the system, and explored the physical limitations of this approach. Efficient and reliable tools, such as this one, are highly valuable for atmospheric Bayesian studies, which need to evaluate a large number of models. We implemented our analytic framework into the \textsc{rate} Python open-source package, available at https://github.com/pcubillos/rate . △ Less

Submitted 31 January, 2019; v1 submitted 11 January, 2019; originally announced January 2019.

Comments: Updated after ApJ proof corrections. The Reproducible Research Compendium is available at https://doi.org/10.5281/zenodo.2529532

The Implications of 3D Thermal Structure on 1D Atmospheric Retrieval

Authors: Jasmina Blecic, Ian Dobbs-Dixon, Thomas Greene

Abstract: Using the atmospheric structure from a 3D global radiation-hydrodynamic simulation of HD 189733b and the open-source BART code, we investigate the difference between the secondary-eclipse temperature structure produced with a 3D simulation and the best-fit 1D retrieved model. Synthetic data are generated by integrating the 3D models over the Spitzer, HST, and JWST bandpasses, covering the waveleng… ▽ More Using the atmospheric structure from a 3D global radiation-hydrodynamic simulation of HD 189733b and the open-source BART code, we investigate the difference between the secondary-eclipse temperature structure produced with a 3D simulation and the best-fit 1D retrieved model. Synthetic data are generated by integrating the 3D models over the Spitzer, HST, and JWST bandpasses, covering the wavelength range between 1 and 11 um. Using the data from different observing instruments, we present detailed comparisons between the temperature-pressure profiles recovered by BART and those from the 3D simulations. We calculate several averages of the 3D thermal structure and implement two temperature parameterizations to investigate different thermal profile shapes. To assess which part of the thermal structure is best constrained by the data, we generate contribution functions for both our theoretical model and each of our retrieved models. Our conclusions are strongly affected by the spectral resolution of the instruments included, their wavelength coverage, and the number of data points combined. We also see some limitations in each of the temperature parametrizations. The results show that our 1D retrieval is recovering a temperature and pressure profile that most closely matches the arithmetic average of the 3D thermal structure. When we use a higher resolution, more data points, and a parametrized temperature profile that allows more flexibility in the middle part of the atmosphere, we find a better match between the retrieved temperature and pressure profile and the arithmetic average. △ Less

Submitted 18 March, 2018; originally announced March 2018.

Comments: 23 pages

Journal ref: The Astrophysical Journal, 848:127 (24pp), 2017 October 20

The Transiting Exoplanet Community Early Release Science Program for JWST

Authors: Jacob L. Bean, Kevin B. Stevenson, Natalie M. Batalha, Zachory Berta-Thompson, Laura Kreidberg, Nicolas Crouzet, Björn Benneke, Michael R. Line, David K. Sing, Hannah R. Wakeford, Heather A. Knutson, Eliza M. -R. Kempton, Jean-Michel Désert, Ian Crossfield, Natasha E. Batalha, Julien de Wit, Vivien Parmentier, Joseph Harrington, Julianne I. Moses, Mercedes Lopez-Morales, Munazza K. Alam, Jasmina Blecic, Giovanni Bruno, Aarynn L. Carter, John W. Chapman , et al. (77 additional authors not shown)

Abstract: The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other… ▽ More The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission. △ Less

Submitted 3 September, 2018; v1 submitted 13 March, 2018; originally announced March 2018.

Comments: PASP in press

Dust in brown dwarfs and extra-solar planets VI. Assessing seed formation across the brown dwarf and exoplanet regimes

Authors: Elspeth K. H. Lee, Jasmina Blecic, Christiane Helling

Abstract: The cloud formation process starts with the formation of seed particles, after which, surface chemical reactions grow or erode the cloud particles. We investigate which materials may form cloud condensation seeds in the gas temperature and pressure regimes (T$_{\rm gas}$ = 100-2000 K, p$_{\rm gas}$ = 10$^{-8}$-100 bar) expected to occur in planetary and brown dwarf atmospheres. We apply modified… ▽ More The cloud formation process starts with the formation of seed particles, after which, surface chemical reactions grow or erode the cloud particles. We investigate which materials may form cloud condensation seeds in the gas temperature and pressure regimes (T$_{\rm gas}$ = 100-2000 K, p$_{\rm gas}$ = 10$^{-8}$-100 bar) expected to occur in planetary and brown dwarf atmospheres. We apply modified classical nucleation theory which requires surface tensions and vapour pressure data for each solid species, which are taken from the literature. We calculate the seed formation rates of TiO$_{2}$[s] and SiO[s] and find that they efficiently nucleate at high temperatures of T$_{\rm gas}$ = 1000-1750 K. Cr[s], KCl[s] and NaCl[s] are found to efficiently nucleate across an intermediate temperature range of T$_{\rm gas}$ = 500-1000 K. We find CsCl[s] may serve as the seed particle for the water cloud layers in cool sub-stellar atmospheres. Four low temperature ice species, H$_{2}$O[s/l], NH$_{3}$[s], H$_{2}$S[s/l] and CH$_{4}$[s], nucleation rates (T$_{\rm gas}$ = 100-250 K) are also investigated for the coolest sub-stellar/planetary atmospheres. Our results suggest a possibly, (T$_{\rm gas}$, p$_{\rm gas}$) distributed hierarchy of seed particle formation regimes throughout the sub-stellar and planetary atmospheric temperature-pressure space. In order to improve the accuracy of the nucleation rate calculation, further research into the small cluster thermochemical data for each cloud species is warranted. The validity of these seed particle scenarios will be tested by applying it to more complete cloud models in the future. △ Less

Submitted 25 January, 2018; originally announced January 2018.

Comments: 23 pages, 16 figures

Journal ref: A&A 614, A126 (2018)

Equilibrium chemistry down to 100 K - Impact of silicates and phyllosilicates on carbon/oxygen ratio

Authors: P. Woitke, Ch. Helling, G. H. Hunter, J. D. Millard, G. E. Turner, M. Worters, J. Blecic, J. W. Stock

Abstract: We introduce a fast and versatile computer code, GGchem, to determine the chemical composition of gases in thermo-chemical equilibrium down to 100 K, with or without equilibrium condensation. We review the data for molecular equilibrium constants, kp(T), from several sources and discuss which functional fits are most suitable for low temperatures. We benchmark our results against another chemical… ▽ More We introduce a fast and versatile computer code, GGchem, to determine the chemical composition of gases in thermo-chemical equilibrium down to 100 K, with or without equilibrium condensation. We review the data for molecular equilibrium constants, kp(T), from several sources and discuss which functional fits are most suitable for low temperatures. We benchmark our results against another chemical equilibrium code. We collect Gibbs free energies, dG(T), for about 200 solid and liquid species from the NIST-JANAF database and the geophysical database SUPCRTBL. We discuss the condensation sequence of the elements with solar abundances in phase equilibrium down to 100 K. Once the major magnesium silicates Mg2SiO4[s] and MgSiO3[s] have formed, the dust/gas mass ratio jumps to a value of about 0.0045 which is significantly lower than the often assumed value of 0.01. Silicate condensation is found to increase the carbon/oxygen ratio (C/O) in the gas from its solar value of ~0.55 up to ~0.71, and, by the additional intake of water and hydroxyl into the solid matrix, the formation of phyllosilicates at temperatures below ~400 K increases the gaseous C/O further to about 0.83. Metallic tungsten (W) is the first condensate found to become thermodynamically stable around 1600 - 2200 K (depending on pressure), several hundreds of Kelvin before subsequent materials like zirconium dioxide (ZrO2) or corundum (Al2O3) can condense. We briefly discuss whether tungsten, despite its low abundance of ~2.E-7 times the silicon abundance, could provide the first seed particles for astrophysical dust formation. The GGchem code is publicly available at https://github.com/pw31/GGchem. △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: accepted by A&A, 18 pages plus 16 pages appendix, 16 figures, and 5 long tables

Journal ref: A&A 614, A1 (2018)

Community targets for JWST's early release science program: evaluation of WASP-63b

Authors: Brian M. Kilpatrick, Patricio E. Cubillos, Kevin B. Stevenson, Nikole K. Lewis, Hannah Wakeford, Ryan J. Macdonald, Nikku Madhusudhan, Jasmina Blecic, Giovanni Bruno, Adam Burrows, Drake Deming, Kevin Heng, Michael R. Line, Caroline V. Morley, Vivien Parmentier, Gregory S. Tucker, Jeff A. Valenti, Ingo P. Waldmann, Jacob L. Bean, Charles Beichman, Jonathan Fraine, J. E. Krick, Joshua D. Lothringer, Avi M. Mandell

Abstract: We present observations of WASP-63b by the Hubble Space Telescope (HST) as part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet for JWST Early Release Science". WASP-63b is one of the community targets under consideration for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. We present a spectrum derived from a single observation by HST Wide Fiel… ▽ More We present observations of WASP-63b by the Hubble Space Telescope (HST) as part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet for JWST Early Release Science". WASP-63b is one of the community targets under consideration for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. We present a spectrum derived from a single observation by HST Wide Field Camera 3 in the near infrared. We engaged groups across the transiting exoplanet community to participate in the analysis of the data and present results from each. There is general agreement amongst all results that we find an H2O absorption feature with 3.5-4.0 sigma significance. However, the feature is muted in comparison to a clear atmosphere at solar composition. Although the detection of the water feature is robust, the reasons for the muting of this feature are ambiguous due to a degeneracy between clouds and composition. The data does not yield robust detections of any molecular species other than H2O. The group was motivated to perform an additional set of retrieval exercises to investigate an apparent bump in the spectrum at ~ 1.55 um. We explore possible disequilibrium chemistry and find this feature is consistent with super-solar HCN abundance but it is questionable if the required mixing ratio of HCN is chemically and physically plausible. The ultimate goal of this study is to vet WASP-63b as a potential community target to best demonstrate the capabilities and systematics of JWST instruments for transiting exoplanet science. In the case of WASP-63b, the presence of a detectable water feature indicates that WASP-63b remains a plausible target for ERS observations. △ Less

Submitted 30 August, 2018; v1 submitted 24 April, 2017; originally announced April 2017.

Comments: repaired broken citations

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

On Correlated-noise Analyses Applied To Exoplanet Light Curves

Authors: Patricio Cubillos, Joseph Harrington, Thomas J. Loredo, Nate B. Lust, Jasmina Blecic, Madison Stemm

Abstract: Time-correlated noise is a significant source of uncertainty when modeling exoplanet light-curve data. A correct assessment of correlated noise is fundamental to determine the true statistical significance of our findings. Here we review three of the most widely used correlated-noise estimators in the exoplanet field, the time-averaging, residual-permutation, and wavelet-likelihood methods. We arg… ▽ More Time-correlated noise is a significant source of uncertainty when modeling exoplanet light-curve data. A correct assessment of correlated noise is fundamental to determine the true statistical significance of our findings. Here we review three of the most widely used correlated-noise estimators in the exoplanet field, the time-averaging, residual-permutation, and wavelet-likelihood methods. We argue that the residual-permutation method is unsound in estimating the uncertainty of parameter estimates. We thus recommend to refrain from this method altogether. We characterize the behavior of the time averaging's rms-vs.-bin-size curves at bin sizes similar to the total observation duration, which may lead to underestimated uncertainties. For the wavelet-likelihood method, we note errors in the published equations and provide a list of corrections. We further assess the performance of these techniques by injecting and retrieving eclipse signals into synthetic and real Spitzer light curves, analyzing the results in terms of the relative-accuracy and coverage-fraction statistics. Both the time-averaging and wavelet-likelihood methods significantly improve the estimate of the eclipse depth over a white-noise analysis (a Markov-chain Monte Carlo exploration assuming uncorrelated noise). However, the corrections are not perfect, when retrieving the eclipse depth from Spitzer datasets, these methods covered the true (injected) depth within the 68\% credible region in only $\sim$45--65\% of the trials. Lastly, we present our open-source model-fitting tool, Multi-Core Markov-Chain Monte Carlo ({MC$^3$}). This package uses Bayesian statistics to estimate the best-fitting values and the credible regions for the parameters for a (user-provided) model. {MC$^3$} is a Python/C code, available at https://github.com/pcubillos/MCcubed. △ Less

Submitted 5 October, 2016; originally announced October 2016.

Comments: Accepted for publication in ApJ

Observations, Thermochemical Calculations, and Modeling of Exoplanetary Atmospheres

Authors: Jasmina Blecic

Abstract: This dissertation as a whole aims to provide means to better understand hot-Jupiter planets through observing, performing thermochemical calculations, and modeling their atmospheres. We used Spitzer multi-wavelength secondary-eclipse observations and targets with high signal-to-noise ratios, as their deep eclipses allow us to detect signatures of spectral features and assess planetary atmospheric… ▽ More This dissertation as a whole aims to provide means to better understand hot-Jupiter planets through observing, performing thermochemical calculations, and modeling their atmospheres. We used Spitzer multi-wavelength secondary-eclipse observations and targets with high signal-to-noise ratios, as their deep eclipses allow us to detect signatures of spectral features and assess planetary atmospheric structure and composition with greater certainty. Chapter 1 gives a short introduction. Chapter 2 presents the Spitzer secondary-eclipse analysis and atmospheric characterization of WASP-14b. WASP-14b is a highly irradiated, transiting hot Jupiter. By applying a Bayesian approach in the atmospheric analysis, we found an absence of thermal inversion contrary to theoretical predictions. Chapter 3 describes the infrared observations of WASP-43b Spitzer secondary eclipses, data analysis, and atmospheric characterization. WASP-43b is one of the closest-orbiting hot Jupiters, orbiting one of the coolest stars with a hot Jupiter. The atmospheric analysis ruled out a strong thermal inversion in its dayside atmosphere. Chapter 4 presents an open-source Thermochemical Equilibrium Abundances (TEA) code and its application to several hot Jupiters. TEA calculates the abundances of gaseous species using the Gibbs free-energy minimization method within an iterative Lagrangian optimization scheme. The code is written in Python and available to the community via http://github.com/dzesmin/TEA. Chapter 5 presents my contributions to an open-source Bayesian Atmospheric Radiative Transfer (BART) code, and its application to WASP-43b. BART characterizes planetary atmospheres based on the observed spectroscopic information. It initializes a planetary atmospheric model, performs radiative-transfer calculations to produce models of planetary spectra, and using a statistical module compares models with observations. △ Less

Submitted 10 April, 2016; originally announced April 2016.

Comments: Ph.D. Dissertation

TEA: A Code for Calculating Thermochemical Equilibrium Abundances

Authors: Jasmina Blecic, Joseph Harrington, M. Oliver Bowman

Abstract: We present an open-source Thermochemical Equilibrium Abundances (TEA) code that calculates the abundances of gaseous molecular species. The code is based on the methodology of White et al. (1958) and Eriksson (1971). It applies Gibbs free-energy minimization using an iterative, Lagrangian optimization scheme. Given elemental abundances, TEA calculates molecular abundances for a particular temperat… ▽ More We present an open-source Thermochemical Equilibrium Abundances (TEA) code that calculates the abundances of gaseous molecular species. The code is based on the methodology of White et al. (1958) and Eriksson (1971). It applies Gibbs free-energy minimization using an iterative, Lagrangian optimization scheme. Given elemental abundances, TEA calculates molecular abundances for a particular temperature and pressure or a list of temperature-pressure pairs. We tested the code against the method of Burrows & Sharp (1999), the free thermochemical equilibrium code CEA (Chemical Equilibrium with Applications), and the example given by White et al. (1958). Using their thermodynamic data, TEA reproduces their final abundances, but with higher precision. We also applied the TEA abundance calculations to models of several hot-Jupiter exoplanets, producing expected results. TEA is written in Python in a modular format. There is a start guide, a user manual, and a code document in addition to this theory paper. TEA is available under a reproducible-research, open-source license via https://github.com/dzesmin/TEA. △ Less

Submitted 23 May, 2015; originally announced May 2015.

Comments: 14 pages, 8 figures, article is submitted to ApJS, posted on arXiv for public commentary, please send comments to the lead author

WASP-8b: Characterization of a Cool and Eccentric Exoplanet with Spitzer

Authors: Patricio Cubillos, Joseph Harrington, Nikku Madhusudhan, Kevin B. Stevenson, Ryan A. Hardy, Jasmina Blecic, David R. Anderson, Matthew Hardin, Christopher J. Campo

Abstract: WASP-8b has 2.18 times Jupiter's mass and is on an eccentric ($e=0.31$) 8.16-day orbit. With a time-averaged equilibrium temperature of 948 K, it is one of the least-irradiated hot Jupiters observed with the Spitzer Space Telescope. We have analyzed six photometric light curves of WASP-8b during secondary eclipse observed in the 3.6, 4.5, and 8.0 {\microns} Infrared Array Camera bands. The eclipse… ▽ More WASP-8b has 2.18 times Jupiter's mass and is on an eccentric ($e=0.31$) 8.16-day orbit. With a time-averaged equilibrium temperature of 948 K, it is one of the least-irradiated hot Jupiters observed with the Spitzer Space Telescope. We have analyzed six photometric light curves of WASP-8b during secondary eclipse observed in the 3.6, 4.5, and 8.0 {\microns} Infrared Array Camera bands. The eclipse depths are $0.113\pm 0.018$%, $0.069\pm 0.007$%, and $0.093\pm 0.023$%, respectively, giving respective brightness temperatures of 1552, 1131, and 938 K. We characterized the atmospheric thermal profile and composition of the planet using a line-by-line radiative transfer code and a Markov-chain Monte Carlo sampler. The data indicated no thermal inversion, independently of any assumption about chemical composition. We noted an anomalously high 3.6-{\microns} brightness temperature (1552 K); by modeling the eccentricity-caused thermal variation, we found that this temperature is plausible for radiative time scales less than $\sim 10^{2}$ hours. However, as no model spectra fit all three data points well, the temperature discrepancy remains as an open question. △ Less

Submitted 21 March, 2013; originally announced March 2013.

Comments: 12 pages, 9 figures, accepted for publication in ApJ

Spitzer observations of the thermal emission from WASP-43b

Authors: Jasmina Blecic, Joseph Harrington, Nikku Madhusudhan, Kevin B. Stevenson, Ryan A. Hardy, Patricio E. Cubillos, Matthew Hardin, Sarah Nymeyer, David R. Anderson, Coel Hellier, Alexis M. S. Smith, Andrew Collier Cameron

Abstract: WASP-43b is one of the closest-orbiting hot Jupiters, with a semimajor axis of a = 0.01526 +/- 0.00018 AU and a period of only 0.81 days. However, it orbits one of the coolest stars with a hot Jupiter (Tstar = 4520 +/- 120 K), giving the planet a modest equilibrium temperature of Teq = 1440 +/- 40 K, assuming zero Bond albedo and uniform planetary energy redistribution. The eclipse depths and brig… ▽ More WASP-43b is one of the closest-orbiting hot Jupiters, with a semimajor axis of a = 0.01526 +/- 0.00018 AU and a period of only 0.81 days. However, it orbits one of the coolest stars with a hot Jupiter (Tstar = 4520 +/- 120 K), giving the planet a modest equilibrium temperature of Teq = 1440 +/- 40 K, assuming zero Bond albedo and uniform planetary energy redistribution. The eclipse depths and brightness temperatures from our jointly fit model are 0.347% +/- 0.013% and 1670 +/- 23 K at 3.6 μm and 0.382% +/- 0.015% and 1514 +/- 25 K at 4.5 μm. The eclipse timings improved the estimate of the orbital period, P, by a factor of three (P = 0.81347436 +/- 1.4*10-7 days) and put an upper limit on the eccentricity (e = 0.010+0.010 -0.007). We use our Spitzer eclipse depths along with four previously reported ground-based photometric observations in the near-infrared to constrain the atmospheric properties of WASP-43b. The data rule out a strong thermal inversion in the dayside atmosphere of WASP-43b. Model atmospheres with no thermal inversions and fiducial oxygen-rich compositions are able to explain all the available data. However, a wide range of metallicities and C/O ratios can explain the data. The data suggest low day-night energy redistribution in the planet, consistent with previous studies, with a nominal upper limit of about 35% for the fraction of energy incident on the dayside that is redistributed to the nightside. △ Less

Submitted 18 February, 2014; v1 submitted 27 February, 2013; originally announced February 2013.

Comments: 11 pages, 9 figures

Journal ref: 2014ApJ...781..116B

Two nearby sub-Earth-sized exoplanet candidates in the GJ 436 system

Authors: Kevin B. Stevenson, Joseph Harrington, Nate B. Lust, Nikole K. Lewis, Guillaume Montagnier, Julianne I. Moses, Channon Visscher, Jasmina Blecic, Ryan A. Hardy, Patricio Cubillos, Christopher J. Campo

Abstract: We report the detection of UCF-1.01, a strong exoplanet candidate with a radius 0.66 +/- 0.04 times that of Earth (R_{\oplus}). This sub-Earth-sized planet transits the nearby M-dwarf star GJ 436 with a period of 1.365862 +/- 8x10^{-6} days. We also report evidence of a 0.65 +/- 0.06 R_{\oplus} exoplanet candidate (labeled UCF-1.02) orbiting the same star with an undetermined period. Using the Spi… ▽ More We report the detection of UCF-1.01, a strong exoplanet candidate with a radius 0.66 +/- 0.04 times that of Earth (R_{\oplus}). This sub-Earth-sized planet transits the nearby M-dwarf star GJ 436 with a period of 1.365862 +/- 8x10^{-6} days. We also report evidence of a 0.65 +/- 0.06 R_{\oplus} exoplanet candidate (labeled UCF-1.02) orbiting the same star with an undetermined period. Using the Spitzer Space Telescope, we measure the dimming of light as the planets pass in front of their parent star to assess their sizes and orbital parameters. If confirmed, UCF-1.01 and UCF-1.02 would be called GJ 436c and GJ 436d, respectively, and would be part of the first multiple-transiting-planet system outside of the Kepler field. Assuming Earth-like densities of 5.515 g/cm^3, we predict both candidates to have similar masses (~0.28 Earth-masses, M_{\oplus}, 2.6 Mars-masses) and surface gravities of ~0.65 g (where g is the gravity on Earth). UCF-1.01's equilibrium temperature (T_{eq}, where emitted and absorbed radiation balance for an equivalent blackbody) is 860 K, making the planet unlikely to harbor life as on Earth. Its weak gravitational field and close proximity to its host star imply that UCF-1.01 is unlikely to have retained its original atmosphere; however, a transient atmosphere is possible if recent impacts or tidal heating were to supply volatiles to the surface. We also present additional observations of GJ 436b during secondary eclipse. The 3.6-micron light curve shows indications of stellar activity, making a reliable secondary eclipse measurement impossible. A second non-detection at 4.5 microns supports our previous work in which we find a methane-deficient and carbon monoxide-rich dayside atmosphere. △ Less

Submitted 17 July, 2012; originally announced July 2012.

Comments: Accepted for publication with ApJ

Infrared Eclipses of the Strongly Irradiated Planet WASP-33b, and Oscillations of its Host Star

Authors: Drake Deming, Jonathan D. Fraine, Pedro V. Sada, Nikku Madhusudhan, Heather A. Knutson, Joseph Harrington, Jasmina Blecic, Sarah Nymeyer, Alexis M. S. Smith, Brian Jackson

Abstract: We observe two secondary eclipses of the strongly irradiated transiting planet WASP-33b in the Ks band, and one secondary eclipse each at 3.6- and 4.5 microns using Warm Spitzer. This planet orbits an A5V delta-Scuti star that is known to exhibit low amplitude non-radial p-mode oscillations at about 0.1-percent semi-amplitude. We detect stellar oscillations in all of our infrared eclipse data, and… ▽ More We observe two secondary eclipses of the strongly irradiated transiting planet WASP-33b in the Ks band, and one secondary eclipse each at 3.6- and 4.5 microns using Warm Spitzer. This planet orbits an A5V delta-Scuti star that is known to exhibit low amplitude non-radial p-mode oscillations at about 0.1-percent semi-amplitude. We detect stellar oscillations in all of our infrared eclipse data, and also in one night of observations at J-band out of eclipse. The oscillation amplitude, in all infrared bands except Ks, is about the same as in the optical. However, the stellar oscillations in Ks band have about twice the amplitude as seen in the optical, possibly because the Brackett-gamma line falls in this bandpass. We use our best-fit values for the eclipse depth, as well as the 0.9 micron eclipse observed by Smith et al., to explore possible states of the exoplanetary atmosphere, based on the method of Madhusudhan and Seager. On this basis we find two possible states for the atmospheric structure of WASP-33b. One possibility is a non-inverted temperature structure in spite of the strong irradiance, but this model requires an enhanced carbon abundance (C/O>1). The alternative model has solar composition, but an inverted temperature structure. Spectroscopy of the planet at secondary eclipse, using a spectral resolution that can resolve the water vapor band structure, should be able to break the degeneracy between these very different possible states of the exoplanetary atmosphere. However, both of those model atmospheres absorb nearly all of the stellar irradiance with minimal longitudinal re-distribution of energy, strengthening the hypothesis of Cowan et al. that the most strongly irradiated planets circulate energy poorly. Our measurement of the central phase of the eclipse yields e*cos(omega)=0.0003 +/-0.00013, which we regard as being consistent with a circular orbit. △ Less

Submitted 4 June, 2012; originally announced June 2012.

Comments: 23 pages, 9 figures, 3 tables, accepted for the Astrophysical Journal

Thermal emission from WASP-24b at 3.6 and 4.5 μm

Authors: A. M. S. Smith, D. R. Anderson, N. Madhusudhan, J. Southworth, A. Collier Cameron, J. Blecic, J. Harrington, C. Hellier, P. F. L. Maxted, D. Pollacco, D. Queloz, B. Smalley, A. H. M . J. Triaud, P. J. Wheatley

Abstract: Aims. We observe occultations of WASP-24b to measure brightness temperatures and to determine whether or not its atmosphere exhibits a thermal inversion (stratosphere). Methods. We observed occultations of WASP-24b at 3.6 and 4.5 μm using the Spitzer Space Telescope. It has been suggested that there is a correlation between stellar activity and the presence of inversions, so we analysed existing H… ▽ More Aims. We observe occultations of WASP-24b to measure brightness temperatures and to determine whether or not its atmosphere exhibits a thermal inversion (stratosphere). Methods. We observed occultations of WASP-24b at 3.6 and 4.5 μm using the Spitzer Space Telescope. It has been suggested that there is a correlation between stellar activity and the presence of inversions, so we analysed existing HARPS spectra in order to calculate log R'HK for WASP-24 and thus determine whether or not the star is chromospherically active. We also observed a transit of WASP-24b in the Strömgren u and y bands, with the CAHA 2.2-m telescope. Results. We measure occultation depths of 0.159 \pm 0.013 per cent at 3.6 μm and 0.202 \pm 0.018 per cent at 4.5 μm. The corresponding planetary brightness temperatures are 1974 \pm 71 K and 1944 \pm 85 K respectively. Atmosphere models with and without a thermal inversion fit the data equally well; we are unable to constrain the presence of an inversion without additional occultation measurements in the near-IR. We find log R'HK = -4.98 \pm 0.12, indicating that WASP-24 is not a chromospherically active star. Our global analysis of new and previously-published data has refined the system parameters, and we find no evidence that the orbit of WASP-24b is non-circular. Conclusions. These results emphasise the importance of complementing Spitzer measurements with observations at shorter wavelengths to gain a full understanding of hot Jupiter atmospheres. △ Less

Submitted 20 August, 2012; v1 submitted 27 March, 2012; originally announced March 2012.

Comments: 7 pages, 4 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 545, A93 (2012)

Thermal Emission of WASP-14b Revealed with Three Spitzer Eclipses

Authors: Jasmina Blecic, Joseph Harrington, Nikku Madhusudhan, Kevin B. Stevenson, Ryan A. Hardy, Patricio Cubillos, Matthew Hardin, Christopher J. Campo, William C. Bowman, Sarah Nymeyer, Thomas J. Loredo, David R. Anderson, Pierre F. L. Maxted

Abstract: Exoplanet WASP-14b is a highly irradiated, transiting hot Jupiter. Joshi et al. calculate an equilibrium temperature Teq of 1866 K for zero albedo and reemission from the entire planet, a mass of 7.3 +/- 0.5 Jupiter masses and a radius of 1.28 +/- 0.08 Jupiter radii. Its mean density of 4.6 g/cm3 is one of the highest known for planets with periods less than 3 days. We obtained three secondary ecl… ▽ More Exoplanet WASP-14b is a highly irradiated, transiting hot Jupiter. Joshi et al. calculate an equilibrium temperature Teq of 1866 K for zero albedo and reemission from the entire planet, a mass of 7.3 +/- 0.5 Jupiter masses and a radius of 1.28 +/- 0.08 Jupiter radii. Its mean density of 4.6 g/cm3 is one of the highest known for planets with periods less than 3 days. We obtained three secondary eclipse light curves with the Spitzer Space Telescope. The eclipse depths from the best jointly fit model are $0.224\%$ +/- $0.018\%$ at 4.5 μm and $0.181\%$ +/- $0.022\%$ at 8.0 μm. The corresponding brightness temperatures are 2212 +/- 94 K and 1590 +/- 116 K. A slight ambiguity between systematic models suggests a conservative 3.6 μm eclipse depth of $0.19\%$ +/- $0.01\%$ and brightness temperature of 2242 +/- 55 K. Although extremely irradiated, WASP-14b does not show any distinct evidence of a thermal inversion. In addition, the present data nominally favor models with day night energy redistribution less than $~30\%$. The current data are generally consistent with oxygen-rich as well as carbon-rich compositions, although an oxygen-rich composition provides a marginally better fit. We confirm a significant eccentricity of e = 0.087 +/- 0.002 and refine other orbital parameters. △ Less

Submitted 19 November, 2013; v1 submitted 9 November, 2011; originally announced November 2011.

Comments: 16 pages, 16 figures

Journal ref: 2013 ApJ 779 5

High C/O Ratio and Weak Thermal Inversion in the Very Hot Atmosphere of Exoplanet WASP-12b

Authors: Nikku Madhusudhan, Joseph Harrington, Kevin B. Stevenson, Sarah Nymeyer, Christopher J. Campo, Peter J. Wheatley, Drake Deming, Jasmina Blecic, Ryan A. Hardy, Nate B. Lust, David R. Anderson, Andrew Collier-Cameron, Christopher B. T. Britt, William C. Bowman, Leslie Hebb, Coel Hellier, Pierre F. L. Maxted, Don Pollacco, Richard G. West

Abstract: The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated interior as opposed to a silicate-dominated composition as found on Earth, and the atmospheres can also differ from those in the Solar System. The solar C/O is 0.54. Here we report an analysis of dayside multi… ▽ More The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated interior as opposed to a silicate-dominated composition as found on Earth, and the atmospheres can also differ from those in the Solar System. The solar C/O is 0.54. Here we report an analysis of dayside multi-wavelength photometry of the transiting hot-Jupiter WASP-12b that reveals C/O >= 1 in its atmosphere. The atmosphere is abundant in CO. It is depleted in water vapor and enhanced in methane by over two orders of magnitude each compared to a solar-abundance chemical-equilibrium model at the expected temperatures. We also find that the extremely irradiated atmosphere (T > 2,500 K) of WASP-12b lacks a prominent thermal inversion, or a stratosphere, and has very efficient day- night energy circulation. The absence of a strong thermal inversion is in stark contrast to theoretical predictions for the most highly irradiated hot-Jupiter atmospheres. △ Less

Submitted 7 December, 2010; originally announced December 2010.

Comments: Discovery of the first carbon-rich planet. To be published in Nature online on 9 Dec 2010. Accepted version is posted here. Please see published paper for final version with minor proof edits

Spitzer Secondary Eclipses of WASP-18b

Authors: Sarah Nymeyer, Joseph Harrington, Ryan A. Hardy, Kevin B. Stevenson, Christopher J. Campo, Nikku Madhusudhan, Andrew Collier-Cameron, Thomas J. Loredo, Jasmina Blecic, William C. Bowman, Christopher B. T. Britt, Patricio Cubillos, Coel Hellier, Michael Gillon, Pierre F. L. Maxted, Leslie Hebb, Peter J. Wheatley, Don Pollacco, David R. Anderson

Abstract: The transiting exoplanet WASP-18b was discovered in 2008 by the Wide Angle Search for Planets (WASP) project. The Spitzer Exoplanet Target of Opportunity Program observed secondary eclipses of WASP-18b using Spitzer's Infrared Array Camera (IRAC) in the 3.6 micron and 5.8 micron bands on 2008 December 20, and in the 4.5 micron and 8.0 micron bands on 2008 December 24. We report eclipse depths of 0… ▽ More The transiting exoplanet WASP-18b was discovered in 2008 by the Wide Angle Search for Planets (WASP) project. The Spitzer Exoplanet Target of Opportunity Program observed secondary eclipses of WASP-18b using Spitzer's Infrared Array Camera (IRAC) in the 3.6 micron and 5.8 micron bands on 2008 December 20, and in the 4.5 micron and 8.0 micron bands on 2008 December 24. We report eclipse depths of 0.30 +/- 0.02%, 0.39 +/- 0.02%, 0.37 +/- 0.03%, 0.41 +/- 0.02%, and brightness temperatures of 3100 +/- 90, 3310 +/- 130, 3080 +/- 140 and 3120 +/- 110 K in order of increasing wavelength. WASP-18b is one of the hottest planets yet discovered - as hot as an M-class star. The planet's pressure-temperature profile most likely features a thermal inversion. The observations also require WASP-18b to have near-zero albedo and almost no redistribution of energy from the day-side to the night side of the planet. △ Less

Submitted 16 August, 2011; v1 submitted 6 May, 2010; originally announced May 2010.

Comments: 10 pages, 13 figures. Electronic supplements to appear with published paper

Journal ref: 2011 ApJ 742 35

On the Orbit of Exoplanet WASP-12b

Authors: Christopher J. Campo, Joseph Harrington, Ryan A. Hardy, Kevin B. Stevenson, Sarah Nymeyer, Darin Ragozzine, Nate B. Lust, David R. Anderson, Andrew Collier-Cameron, Jasmina Blecic, Christopher B. T. Britt, William C. Bowman, Peter J. Wheatley, Thomas J. Loredo, Drake Deming, Leslie Hebb, Coel Hellier, Pierre F. L. Maxted, Don Pollaco, Richard G. West

Abstract: We observed two secondary eclipses of the exoplanet WASP-12b using the Infrared Array Camera on the Spitzer Space Telescope. The close proximity of WASP-12b to its G-type star results in extreme tidal forces capable of inducing apsidal precession with a period as short as a few decades. This precession would be measurable if the orbit had a significant eccentricity, leading to an estimate of the t… ▽ More We observed two secondary eclipses of the exoplanet WASP-12b using the Infrared Array Camera on the Spitzer Space Telescope. The close proximity of WASP-12b to its G-type star results in extreme tidal forces capable of inducing apsidal precession with a period as short as a few decades. This precession would be measurable if the orbit had a significant eccentricity, leading to an estimate of the tidal Love number and an assessment of the degree of central concentration in the planetary interior. An initial ground-based secondary eclipse phase reported by Lopez-Morales et al. (0.510 +/- 0.002) implied eccentricity at the 4.5 sigma level. The spectroscopic orbit of Hebb et al. has eccentricity 0.049 +/- 0.015, a 3 sigma result, implying an eclipse phase of 0.509 +/- 0.007. However, there is a well documented tendency of spectroscopic data to overestimate small eccentricities. Our eclipse phases are 0.5010 +/- 0.0006 (3.6 and 5.8 microns) and 0.5006 +/- 0.0007 (4.5 and 8.0 microns). An unlikely orbital precession scenario invoking an alignment of the orbit during the Spitzer observations could have explained this apparent discrepancy, but the final eclipse phase of Lopez-Morales et al. (0.510 -0.006 / +0.007) is consistent with a circular orbit at better than 2 sigma. An orbit fit to all the available transit, eclipse, and radial-velocity data indicates precession at <1 sigma; a non-precessing solution fits better. We also comment on analysis and reporting for Spitzer exoplanet data in light of recent re-analyses. △ Less

Submitted 9 December, 2010; v1 submitted 14 March, 2010; originally announced March 2010.

Comments: 12 pages, 8 figures, electronic supplements to appear with published paper