Showing 1–50 of 170 results for author: Lewis, N

Latitudinal Asymmetry in the Dayside Atmosphere of WASP-43b

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

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

Submitted 14 June, 2024; originally announced June 2024.

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

Identifying and Fitting Eclipse Maps of Exoplanets with Cross-Validation

Authors: Mark Hammond, Neil T. Lewis, Sasha Boone, Xueqing Chen, João M. Mendonça, Vivien Parmentier, Jake Taylor, Taylor Bell, Leonardo dos Santos, Nicolas Crouzet, Laura Kreidberg, Michael Radica, Michael Zhang

Abstract: Eclipse mapping uses the shape of the eclipse of an exoplanet to measure its two-dimensional structure. Light curves are mostly composed of longitudinal information, with the latitudinal information only contained in the brief ingress and egress of the eclipse. This imbalance can lead to a spuriously confident map, where the longitudinal structure is constrained by out-of-eclipse data and the lati… ▽ More Eclipse mapping uses the shape of the eclipse of an exoplanet to measure its two-dimensional structure. Light curves are mostly composed of longitudinal information, with the latitudinal information only contained in the brief ingress and egress of the eclipse. This imbalance can lead to a spuriously confident map, where the longitudinal structure is constrained by out-of-eclipse data and the latitudinal structure is wrongly determined by the priors on the map. We present a new method to address this issue. The method tests for the presence of an eclipse mapping signal by using k-fold cross-validation to compare the performance of a simple mapping model to the null hypothesis of a uniform disk. If a signal is found, the method fits a map with more degrees of freedom, optimising its information content. The information content is varied by penalising the model likelihood by a factor proportional to the spatial entropy of the map, optimised by cross-validation. We demonstrate this method for simulated datasets then apply it to three observational datasets. The method identifies an eclipse mapping signal for JWST MIRI/LRS observations of WASP-43b but does not identify a signal for JWST NIRISS/SOSS observations of WASP-18b or Spitzer Space Telescope observations of HD 189733b. It is possible to fit eclipse maps to these datasets, but we suggest that these maps are overfitting the eclipse shape. We fit a new map with more spatial freedom to the WASP-43b dataset and show a flatter east-west structure than previously derived. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

JWST-TST High Contrast: JWST/NIRCam observations of the young giant planet $β$ Pic b

Authors: Jens Kammerer, Kellen Lawson, Marshall D. Perrin, Isabel Rebollido, Christopher C. Stark, Tomas Stolker, Julien H. Girard, Laurent Pueyo, William O. Balmer, Kadin Worthen, Christine Chen, Roeland P. van der Marel, Nikole K. Lewis, Kimberly Ward-Duong, Jeff A. Valenti, Mark Clampin, C. Matt Mountain

Abstract: We present the first JWST/NIRCam observations of the directly-imaged gas giant exoplanet $β$ Pic b. Observations in six filters using NIRCam's round coronagraphic masks provide a high signal-to-noise detection of $β$ Pic b and the archetypal debris disk around $β$ Pic over a wavelength range of $\sim$1.7-5 $μ$m. This paper focuses on the detection of $β$ Pic b and other potential point sources in… ▽ More We present the first JWST/NIRCam observations of the directly-imaged gas giant exoplanet $β$ Pic b. Observations in six filters using NIRCam's round coronagraphic masks provide a high signal-to-noise detection of $β$ Pic b and the archetypal debris disk around $β$ Pic over a wavelength range of $\sim$1.7-5 $μ$m. This paper focuses on the detection of $β$ Pic b and other potential point sources in the NIRCam data, following a paper by Rebollido et al. which presented the NIRCam and MIRI view of the debris disk around $β$ Pic. We develop and validate approaches for obtaining accurate photometry of planets in the presence of bright, complex circumstellar backgrounds. By simultaneously fitting the planet's PSF and a geometric model for the disk, we obtain planet photometry that is in good agreement with previous measurements from the ground. The NIRCam data supports the cloudy nature of $β$ Pic b's atmosphere and the discrepancy between its mass as inferred from evolutionary models and the dynamical mass reported in the literature. We further identify five additional localized sources in the data, but all of them are found to be background stars or galaxies based on their color or spatial extent. We can rule out additional planets in the disk midplane above 1 Jupiter mass outward of 2 arcsec ($\sim$40 au) and away from the disk midplane above 0.05 Jupiter masses outward of 4 arcsec ($\sim$80 au). The inner giant planet $β$ Pic c remains undetected behind the coronagraphic masks of NIRCam in our observations. △ Less

Submitted 3 July, 2024; v1 submitted 28 May, 2024; originally announced May 2024.

Comments: 30 pages, 14 figures, published in AJ

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

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

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

Submitted 17 May, 2024; originally announced May 2024.

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

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

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

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

Submitted 25 April, 2024; originally announced April 2024.

Comments: Accepted for publication in The Astronomical Journal

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

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

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

Submitted 23 January, 2024; originally announced January 2024.

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

JWST-TST High Contrast: Asymmetries, dust populations and hints of a collision in the $β$ Pictoris disk with NIRCam and MIRI

Authors: Isabel Rebollido, Christopher C. Stark, Jens Kammerer, Marshall D. Perrin, Kellen Lawson, Laurent Pueyo, Christine Chen, Dean Hines, Julien H. Girard, Kadin Worthen, Carl Ingerbretsen, Sarah Betti, Mark Clampin, David Golimowski, Kielan Hoch, Nikole K. Lewis, Cicero X. Lu, Roeland P. van der Marel, Emily Rickman, Sara Seager, Remi Soummer, Jeff A. Valenti, Kimberly Ward-Duong, C. Matt Mountain

Abstract: We present the first JWST MIRI and NIRCam observations of the prominent debris disk around Beta Pictoris. Coronagraphic observations in 8 filters spanning from 1.8 to 23~$μ$m provide an unprecedentedly clear view of the disk at these wavelengths. The objectives of the observing program were to investigate the dust composition and distribution, and to investigate the presence of planets in the syst… ▽ More We present the first JWST MIRI and NIRCam observations of the prominent debris disk around Beta Pictoris. Coronagraphic observations in 8 filters spanning from 1.8 to 23~$μ$m provide an unprecedentedly clear view of the disk at these wavelengths. The objectives of the observing program were to investigate the dust composition and distribution, and to investigate the presence of planets in the system. In this paper, we focus on the disk components, providing surface brightness measurements for all images and a detailed investigation of the asymmetries observed. A companion paper by Kammerer et al. will focus on the planets in this system using the same data. We report for the first time the presence of an extended secondary disk in thermal emission, with a curved extension bent away from the plane of the disk. This feature, which we refer to as the ``cat's tail", seems to be connected with the previously reported CO clump, mid-infrared asymmetry detected in the southwest side, and the warp observed in scattered light. We present a model of this secondary disk sporadically producing dust that broadly reproduces the morphology, flux, and color of the cat's tail, as well as other features observed in the disk, and suggests the secondary disk is composed largely of porous, organic refractory dust grains. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: Accepted to be published in AJ. 37 pages, 21 figures

The impact of spectral line wing cut-off: Recommended standard method with application to MAESTRO opacity database

Authors: Ehsan, Gharib-Nezhad, Natasha E. Batalha, Katy Chubb, Richard Freedman, Iouli E. Gordon, Robert R. Gamache, Robert J. Hargreaves, Nikole K. Lewis, Jonathan Tennyson, Sergei N. Yurchenko

Abstract: When computing cross-sections from a line list, the result depends not only on the line strength, but also the line shape, pressure-broadening parameters, and line wing cut-off (i.e., the maximum distance calculated from each line centre). Pressure-broadening can be described using the Lorentz lineshape, but it is known to not represent the true absorption in the far wings. Both theory and experim… ▽ More When computing cross-sections from a line list, the result depends not only on the line strength, but also the line shape, pressure-broadening parameters, and line wing cut-off (i.e., the maximum distance calculated from each line centre). Pressure-broadening can be described using the Lorentz lineshape, but it is known to not represent the true absorption in the far wings. Both theory and experiment have shown that far from the line centre, non-Lorentzian behaviour controls the shape of the wings and the Lorentz lineshape fails to accurately characterize the absorption, leading to an underestimation or overestimation of the opacity continuum depending on the molecular species involved. The line wing cut-off is an often overlooked parameter when calculating absorption cross sections, but can have a significant effect on the appearance of the spectrum since it dictates the extent of the line wing that contributes to the calculation either side of every line centre. Therefore, when used to analyse exoplanet and brown dwarf spectra, an inaccurate choice for the line wing cut-off can result in errors in the opacity continuum, which propagate into the modeled transit spectra, and ultimately impact/bias the interpretation of observational spectra, and the derived composition and thermal structure. Here, we examine the different methods commonly utilized to calculate the wing cut-off and propose a standard practice procedure (i.e., absolute value of 25~cm$^{-1}$ for $P\leqslant$~200~bar and 100~cm$^{-1}$ for $P >$ ~200~bar) to generate molecular opacities which will be used by the open-access {\tt MAESTRO} (Molecules and Atoms in Exoplanet Science: Tools and Resources for Opacities) database. The pressing need for new measurements and theoretical studies of the far-wings is highlighted. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: 13 pages, 6 figures

Journal ref: RAS Techniques and Instruments, 2023, rzad058

Near-Infrared Transmission Spectroscopy of HAT-P-18$\,$b with NIRISS: Disentangling Planetary and Stellar Features in the Era of JWST

Authors: Marylou Fournier-Tondreau, Ryan J. MacDonald, Michael Radica, David Lafrenière, Luis Welbanks, Caroline Piaulet, Louis-Philippe Coulombe, Romain Allart, Kim Morel, Étienne Artigau, Loïc Albert, Olivia Lim, René Doyon, Björn Benneke, Jason F. Rowe, Antoine Darveau-Bernier, Nicolas B. Cowan, Nikole K. Lewis, Neil James Cook, Laura Flagg, Frédéric Genest, Stefan Pelletier, Doug Johnstone, Lisa Dang, Lisa Kaltenegger , et al. (2 additional authors not shown)

Abstract: The JWST Early Release Observations (ERO) included a NIRISS/SOSS (0.6-2.8$\,μ$m) transit of the $\sim\,$850$\,$K Saturn-mass exoplanet HAT-P-18$\,$b. Initial analysis of these data reported detections of water, escaping helium, and haze. However, active K dwarfs like HAT-P-18 possess surface heterogeneities $-$ starspots and faculae $-$ that can complicate the interpretation of transmission spectr… ▽ More The JWST Early Release Observations (ERO) included a NIRISS/SOSS (0.6-2.8$\,μ$m) transit of the $\sim\,$850$\,$K Saturn-mass exoplanet HAT-P-18$\,$b. Initial analysis of these data reported detections of water, escaping helium, and haze. However, active K dwarfs like HAT-P-18 possess surface heterogeneities $-$ starspots and faculae $-$ that can complicate the interpretation of transmission spectra, and indeed, a spot-crossing event is present in HAT-P-18$\,$b's NIRISS/SOSS light curves. Here, we present an extensive reanalysis and interpretation of the JWST ERO transmission spectrum of HAT-P-18$\,$b, as well as HST/WFC3 and $\textit{Spitzer}$/IRAC transit observations. We detect H$_2$O (12.5$\,σ$), CO$_2$ (7.3$\,σ$), a cloud deck (7.4$\,σ$), and unocculted starspots (5.8$\,σ$), alongside hints of Na (2.7$\,σ$). We do not detect the previously reported CH$_4$ ($\log$ CH$_4$ $<$ -6 to 2$\,σ$). We obtain excellent agreement between three independent retrieval codes, which find a sub-solar H$_2$O abundance ($\log$ H$_2$O $\approx -4.4 \pm 0.3$). However, the inferred CO$_2$ abundance ($\log$ CO$_2$ $\approx -4.8 \pm 0.4$) is significantly super-solar and requires further investigation into its origin. We also introduce new stellar heterogeneity considerations by fitting for the active regions' surface gravities $-$ a proxy for the effects of magnetic pressure. Finally, we compare our JWST inferences to those from HST/WFC3 and $\textit{Spitzer}$/IRAC. Our results highlight the exceptional promise of simultaneous planetary atmosphere and stellar heterogeneity constraints in the era of JWST and demonstrate that JWST transmission spectra may warrant more complex treatments of the transit light source effect. △ Less

Submitted 18 December, 2023; v1 submitted 23 October, 2023; originally announced October 2023.

JWST-TST High Contrast: Achieving direct spectroscopy of faint substellar companions next to bright stars with the NIRSpec IFU

Authors: Jean-Baptiste Ruffio, Marshall D. Perrin, Kielan K. W. Hoch, Jens Kammerer, Quinn M. Konopacky, Laurent Pueyo, Alex Madurowicz, Emily Rickman, Christopher A. Theissen, Shubh Agrawal, Alexandra Z. Greenbaum, Brittany E. Miles, Travis S. Barman, William O. Balmer, Jorge Llop-Sayson, Julien H. Girard, Isabel Rebollido, Rémi Soummer, Natalie H. Allen, Jay Anderson, Charles A. Beichman, Andrea Bellini, Geoffrey Bryden, Néstor Espinoza, Ana Glidden , et al. (11 additional authors not shown)

Abstract: The JWST NIRSpec integral field unit (IFU) presents a unique opportunity to observe directly imaged exoplanets from 3-5 um at moderate spectral resolution (R~2,700) and thereby better constrain the composition, disequilibrium chemistry, and cloud properties of their atmospheres. In this work, we present the first NIRSpec IFU high-contrast observations of a substellar companion that requires starli… ▽ More The JWST NIRSpec integral field unit (IFU) presents a unique opportunity to observe directly imaged exoplanets from 3-5 um at moderate spectral resolution (R~2,700) and thereby better constrain the composition, disequilibrium chemistry, and cloud properties of their atmospheres. In this work, we present the first NIRSpec IFU high-contrast observations of a substellar companion that requires starlight suppression techniques. We develop specific data reduction strategies to study faint companions around bright stars, and assess the performance of NIRSpec at high contrast. First, we demonstrate an approach to forward model the companion signal and the starlight directly in the detector images, which mitigates the effects of NIRSpec's spatial undersampling. We demonstrate a sensitivity to planets that are 3e-6 fainter than their stars at 1'', or 3e-5 at 0.3''. Then, we implement a reference star point spread function (PSF) subtraction and a spectral extraction that does not require spatially and spectrally regularly sampled spectral cubes. This allows us to extract a moderate resolution (R~2,700) spectrum of the faint T-dwarf companion HD 19467 B from 2.9-5.2 um with signal-to-noise ratio (S/N)~10 per resolution element. Across this wavelength range, HD~19467~B has a flux ratio varying between 1e-5-1e-4 and a separation relative to its star of 1.6''. A companion paper by Hoch et al. more deeply analyzes the atmospheric properties of this companion based on the extracted spectrum. Using the methods developed here, NIRSpec's sensitivity may enable direct detection and spectral characterization of relatively old (~1 Gyr), cool (~250 K), and closely separated (~3-5 au) exoplanets that are less massive than Jupiter. △ Less

Submitted 31 May, 2024; v1 submitted 15 October, 2023; originally announced October 2023.

Comments: Accepted to AJ. The data analysis scripts for this work are published https://github.com/jruffio/HD_19467_B (https://zenodo.org/doi/10.5281/zenodo.11391740). The main revisions of the manuscript are listed in the change history section of the readme

JWST-TST DREAMS: Quartz Clouds in the Atmosphere of WASP-17b

Authors: David Grant, Nikole K. Lewis, Hannah R. Wakeford, Natasha E. Batalha, Ana Glidden, Jayesh Goyal, Elijah Mullens, Ryan J. MacDonald, Erin M. May, Sara Seager, Kevin B. Stevenson, Jeff A. Valenti, Channon Visscher, Lili Alderson, Natalie H. Allen, Caleb I. Cañas, Knicole Colón, Mark Clampin, Néstor Espinoza, Amélie Gressier, Jingcheng Huang, Zifan Lin, Douglas Long, Dana R. Louie, Maria Peña-Guerrero , et al. (17 additional authors not shown)

Abstract: Clouds are prevalent in many of the exoplanet atmospheres that have been observed to date. For transiting exoplanets, we know if clouds are present because they mute spectral features and cause wavelength-dependent scattering. While the exact composition of these clouds is largely unknown, this information is vital to understanding the chemistry and energy budget of planetary atmospheres. In this… ▽ More Clouds are prevalent in many of the exoplanet atmospheres that have been observed to date. For transiting exoplanets, we know if clouds are present because they mute spectral features and cause wavelength-dependent scattering. While the exact composition of these clouds is largely unknown, this information is vital to understanding the chemistry and energy budget of planetary atmospheres. In this work, we observe one transit of the hot Jupiter WASP-17b with JWST's MIRI LRS and generate a transmission spectrum from 5-12 $\rmμ$m. These wavelengths allow us to probe absorption due to the vibrational modes of various predicted cloud species. Our transmission spectrum shows additional opacity centered at 8.6 $\rmμ$m, and detailed atmospheric modeling and retrievals identify this feature as SiO$_2$(s) (quartz) clouds. The SiO$_2$(s) clouds model is preferred at 3.5-4.2$σ$ versus a cloud-free model and at 2.6$σ$ versus a generic aerosol prescription. We find the SiO$_2$(s) clouds are comprised of small ${\sim}0.01$ $\rmμ$m particles, which extend to high altitudes in the atmosphere. The atmosphere also shows a depletion of H$_2$O, a finding consistent with the formation of high-temperature aerosols from oxygen-rich species. This work is part of a series of studies by our JWST Telescope Scientist Team (JWST-TST), in which we will use Guaranteed Time Observations to perform Deep Reconnaissance of Exoplanet Atmospheres through Multi-instrument Spectroscopy (DREAMS). △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 19 pages, 7 figures, accepted for publication in ApJL

Equatorial Waves and Superrotation in the Stratosphere of a Titan General Circulation Model

Authors: Neil T. Lewis, Nicholas A. Lombardo, Peter L. Read, Juan M. Lora

Abstract: We investigate the characteristics of equatorial waves associated with the maintenance of superrotation in the stratosphere of a Titan general circulation model. A variety of equatorial waves are present in the model atmosphere, including equatorial Kelvin waves, equatorial Rossby waves, and mixed Rossby-gravity waves. In the upper stratosphere, acceleration of superrotation is strongest around so… ▽ More We investigate the characteristics of equatorial waves associated with the maintenance of superrotation in the stratosphere of a Titan general circulation model. A variety of equatorial waves are present in the model atmosphere, including equatorial Kelvin waves, equatorial Rossby waves, and mixed Rossby-gravity waves. In the upper stratosphere, acceleration of superrotation is strongest around solstice and is due to interaction between equatorial Kelvin waves and Rossby-type waves in winter-hemisphere mid-latitudes. The existence of this 'Rossby-Kelvin'-type wave appears to depend on strong meridional shear of the background zonal wind that occurs in the upper stratosphere at times away from the equinoxes. In the lower stratosphere, acceleration of superrotation occurs throughout the year and is partially induced by equatorial Rossby waves, which we speculate are generated by quasigeostrophic barotropic instability. Acceleration of superrotation is generally due to waves with phase speeds close to the zonal velocity of the mean flow. Consequently, they have short vertical wavelengths which are close to the model's vertical grid scale, and therefore are likely to be not properly represented. We suggest this may be a common issue amongst Titan GCMs which should be addressed by future model development. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Comments: Accepted for publication in The Planetary Science Journal

HST/WFC3 Light Curve Supports a Terrestrial Composition for the Closest Exoplanet to Transit an M Dwarf

Authors: Emily K Pass, Jennifer G Winters, David Charbonneau, Aurelia Balkanski, Nikole Lewis, Maura Lally, Jacob L Bean, Ryan Cloutier, Jason D Eastman

Abstract: Previous studies of the exoplanet LTT 1445Ac concluded that the light curve from the Transiting Exoplanet Survey Satellite (TESS) was consistent with both grazing and non-grazing geometries. As a result, the radius and hence density of the planet remained unknown. To resolve this ambiguity, we observed the LTT 1445 system for six spacecraft orbits of the Hubble Space Telescope (HST) using WFC3/UVI… ▽ More Previous studies of the exoplanet LTT 1445Ac concluded that the light curve from the Transiting Exoplanet Survey Satellite (TESS) was consistent with both grazing and non-grazing geometries. As a result, the radius and hence density of the planet remained unknown. To resolve this ambiguity, we observed the LTT 1445 system for six spacecraft orbits of the Hubble Space Telescope (HST) using WFC3/UVIS imaging in spatial scan mode, including one partial transit of LTT 1445Ac. This imaging produces resolved light curves of each of the three stars in the LTT 1445 system. We confirm that the planet transits LTT 1445A and that LTT 1445C is the source of the rotational modulation seen in the TESS light curve, and we refine the estimate of the dilution factor for the TESS data. We perform a joint fit to the TESS and HST observations, finding that the transit of LTT 1445Ac is not grazing with 97% confidence. We measure a planetary radius of 1.07$_{-0.07}^{+0.10}$ R$_\oplus$. Combined with previous radial velocity observations, our analysis yields a planetary mass of $1.37\pm0.19$ M$_\oplus$ and a planetary density of 5.9$_{-1.5}^{+1.8}$ g cm$^{-3}$. LTT 1445Ac is likely an Earth analog with respect to its mass and radius, albeit with a higher instellation, and is therefore an exciting target for future atmospheric studies. △ Less

Submitted 21 August, 2023; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: Accepted for publication in AJ. 10 pages, 7 figures, 3 tables

A Shallow Water Model Exploration of Atmospheric Circulation on Sub-Neptunes: Effects of Radiative Forcing and Rotation Period

Authors: Ekaterina Landgren, Alice Nadeau, Nikole Lewis, Tiffany Kataria, Peter Hitchcock

Abstract: Sub-Neptune type exoplanets are abundant in our galaxy yet have no solar system analogs. They exist in a broad range of stellar forcing and rotational regimes that are distinctly different from solar system planets and more commonly studied hot Jupiters. Here we present simulations that explore global atmospheric circulation of sub-Neptunes generated with a two-dimensional shallow-water model, SWA… ▽ More Sub-Neptune type exoplanets are abundant in our galaxy yet have no solar system analogs. They exist in a broad range of stellar forcing and rotational regimes that are distinctly different from solar system planets and more commonly studied hot Jupiters. Here we present simulations that explore global atmospheric circulation of sub-Neptunes generated with a two-dimensional shallow-water model, SWAMPE. We explore the circulation regimes of synchronously rotating sub-Neptunes with a focus on the interaction of planetary rotation rate and radiative timescale in a variety of stellar insolations. In highly irradiated, short-timescale regimes, our models exhibit high day-night geopotential contrasts. As the timescales become longer, the geopotential contrasts and longitudinal variability decrease, while temporal variability increases. The transition from day-to-night flow to jet-dominated flow is primarily driven by the radiative timescale. Strong- and medium-forcing regimes exhibit transitions between day-to-night flow and jet-dominated flow at similar points in the parameter space. Weak-forcing regime differs due to comparatively stronger rotational effects. Planetary rotation period dominates in determining equator-to-pole geopotential contrast. Our simulations exhibit higher time variability when either radiative timescale or rotation period is long. △ Less

Submitted 30 May, 2023; originally announced May 2023.

Comments: Accepted for publication in PSJ on 5/12/2023. 25 pages, 1 table, 11 figures

Mass Derivation of planets K2-21b and K2-21c from Transit Timing Variations

Authors: Maryame El Moutamid, Kevin B. Stevenson, Billy Quarles, Nikole K. Lewis, Erik Petigura Daniel Fabrycky, Jacob L. Bean, Diana Dragomir, Kristin S. Sotzenvand Michael W. Werner

Abstract: While various indirect methods are used to detect exoplanets, one of the most effective and accurate methods is the transit method, which measures the brightness of a given star for periodic dips when an exoplanet is passing in front of the parent star. For systems with multiple transiting planets, the gravitational perturbations between planets affect their transit times. The difference in transi… ▽ More While various indirect methods are used to detect exoplanets, one of the most effective and accurate methods is the transit method, which measures the brightness of a given star for periodic dips when an exoplanet is passing in front of the parent star. For systems with multiple transiting planets, the gravitational perturbations between planets affect their transit times. The difference in transit times allows a measurement of the planet masses and orbital eccentricities. These parameters help speculating on the formation, evolution and stability of the system. Using Transit Timing Variations (TTVs), we measure the masses and eccentricities of two planets orbiting K2-21, a relatively bright K7 dwarf star. These two planets exhibit measurable TTVs, have orbital periods of about 9.32 days and 15.50 days, respectively, and a period ratio of about 1.66, which is relatively near to the 5:3 mean motion resonance. We report that the inner and outer planets in the K2-21 system have properties consistent with the presence of a hydrogen and helium dominated atmospheres, as we estimate their masses to be 1.59^{+0.52}_{-0.44} M_E and 3.88^{+1.22}_{-1.07} M_E and densities of 0.22^{+0.05}_{-0.04} rho_E and 0.34^{+0.08}_{-0.06} rho_E, respectively (M_E and rho_E are the mass and density of Earth, respectively). Our results show that the inner planet is less dense than the outer planet; one more counter-intuitive exoplanetary system such as Kepler-105, LTT 1445, TOI-175 and Kepler-279 systems. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Schedule optimization for transiting exoplanet observations with NASA's Pandora SmallSat mission

Authors: Trevor O. Foote, Thomas Barclay, Christina L. Hedges, Nikole K. Lewis, Elisa V. Quintana, Benjamin V. Rackham

Abstract: Pandora is an upcoming NASA SmallSat mission that will observe transiting exoplanets to study their atmospheres and the variability of their host stars. Efficient mission planning is critical for maximizing the science achieved with the year-long primary mission. To this end, we have developed a scheduler based on a metaheuristic algorithm that is focused on tackling the unique challenges of time-… ▽ More Pandora is an upcoming NASA SmallSat mission that will observe transiting exoplanets to study their atmospheres and the variability of their host stars. Efficient mission planning is critical for maximizing the science achieved with the year-long primary mission. To this end, we have developed a scheduler based on a metaheuristic algorithm that is focused on tackling the unique challenges of time-constrained observing missions, like Pandora. Our scheduling algorithm combines a minimum transit requirement metric, which ensures we meet observational requirements, with a `quality' metric that considers three factors to determine the scientific quality of each observation window around an exoplanet transit (defined as a visit). These three factors are: observing efficiency during a visit, the amount of the transit captured by the telescope during a visit, and how much of the transit captured is contaminated by a coincidental passing of the observatory through the South Atlantic Anomaly. The importance of each of these factors can be adjusted based on the needs or preferences of the science team. Utilizing this schedule optimizer, we develop and compare a few schedules with differing factor weights for the Pandora SmallSat mission, illustrating trade-offs that should be considered between the three quality factors. We also find that under all scenarios probed, Pandora will not only be able to achieve its observational requirements using the planets on the notional target list but will do so with significant time remaining for ancillary science. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: 35 pages, 7 figures. Submitted to JATIS, SPIE. Python code is available at: https://github.com/PandoraMission/pandora-scheduler

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

JWST-TST Proper Motions: I. High-Precision NIRISS Calibration and Large Magellanic Cloud Kinematics

Authors: M. Libralato, A. Bellini, R. P. van der Marel, J. Anderson, S. T. Sohn, L. L. Watkins, L. Alderson, N. Allen, M. Clampin, A. Glidden, J. Goyal, K. Hoch, J. Huang, J. Kammerer, N. K. Lewis, Z. Lin, D. Long, D. Louie, R. J. MacDonald, M. Mountain, M. Peña-Guerrero, M. D. Perrin, L. Pueyo, I. Rebollido, E. Rickman , et al. (5 additional authors not shown)

Abstract: We develop and disseminate effective point-spread functions and geometric-distortion solutions for high-precision astrometry and photometry with the JWST NIRISS instrument. We correct field dependencies and detector effects, and assess the quality and the temporal stability of the calibrations. As a scientific application and validation, we study the proper motion (PM) kinematics of stars in the J… ▽ More We develop and disseminate effective point-spread functions and geometric-distortion solutions for high-precision astrometry and photometry with the JWST NIRISS instrument. We correct field dependencies and detector effects, and assess the quality and the temporal stability of the calibrations. As a scientific application and validation, we study the proper motion (PM) kinematics of stars in the JWST calibration field near the Large Magellanic Cloud (LMC) center, comparing to a first-epoch Hubble Space Telescope (HST) archival catalog with a 16-yr baseline. For stars with G~20, the median PM uncertainty is ~13 $μ$as yr$^{-1}$ (3.1 km s$^{-1}$), better than Gaia DR3 typically achieves for its very best-measured stars. We kinematically detect the known star cluster OGLE-CL LMC 407, measure its absolute PM for the first time, and show how this differs from other LMC populations. The inferred cluster dispersion sets an upper limit of 24 $μ$as yr$^{-1}$ (5.6 km s$^{-1}$) on systematic uncertainties. Red-giant-branch stars have a velocity dispersion of 33.8 $\pm$ 0.6 km s$^{-1}$, while younger blue populations have a narrower velocity distribution, but with a significant kinematical substructure. We discuss how this relates to the larger velocity dispersions inferred from Gaia DR3. These results establish JWST as capable of state-of-the-art astrometry, building on the extensive legacy of HST. This is the first paper in a series by our JWST Telescope Scientist Team (TST), in which we will use Guaranteed Time Observations to study the PM kinematics of various stellar systems in the Local Group. △ Less

Submitted 25 April, 2023; v1 submitted 28 February, 2023; originally announced March 2023.

Comments: 28 pages, 25 figures, 2 tables. Accepted for publication in ApJ. The effective point-spread-function models, the geometric-distortion solutions and a preliminary version of the code are available at the links provided in the manuscript

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

A JWST NIRSpec Phase Curve for WASP-121b: Dayside Emission Strongest Eastward of the Substellar Point and Nightside Conditions Conducive to Cloud Formation

Authors: Thomas Mikal-Evans, David K. Sing, Jiayin Dong, Daniel Foreman-Mackey, Tiffany Kataria, Joanna K. Barstow, Jayesh M. Goyal, Nikole K. Lewis, Joshua D. Lothringer, Nathan J. Mayne, Hannah R. Wakeford, Duncan A. Christie, Zafar Rustamkulov

Abstract: We present the first exoplanet phase curve measurement made with the JWST NIRSpec instrument, highlighting the exceptional stability of this newly-commissioned observatory for exoplanet climate studies. The target, WASP-121b, is an ultrahot Jupiter with an orbital period of 30.6 hr. We analyze two broadband light curves generated for the NRS1 and NRS2 detectors, covering wavelength ranges of 2.70-… ▽ More We present the first exoplanet phase curve measurement made with the JWST NIRSpec instrument, highlighting the exceptional stability of this newly-commissioned observatory for exoplanet climate studies. The target, WASP-121b, is an ultrahot Jupiter with an orbital period of 30.6 hr. We analyze two broadband light curves generated for the NRS1 and NRS2 detectors, covering wavelength ranges of 2.70-3.72 micron and 3.82-5.15 micron, respectively. Both light curves exhibit minimal systematics, with approximately linear drifts in the baseline flux level of 30 ppm/hr (NRS1) and 10 ppm/hr (NRS2). Assuming a simple brightness map for the planet described by a low-order spherical harmonic dipole, our light curve fits suggest that the phase curve peaks coincide with orbital phases $3.36 \pm 0.11$ deg (NRS1) and $2.66 \pm 0.12$ deg (NRS2) prior to mid-eclipse. This is consistent with the strongest dayside emission emanating from eastward of the substellar point. We measure planet-to-star emission ratios of $3,924 \pm 7$ ppm (NRS1) and $4,924 \pm 9$ ppm (NRS2) for the dayside hemisphere, and $136 \pm 8$ ppm (NRS1) and $630 \pm 10$ ppm (NRS2) for the nightside hemisphere. The latter nightside emission ratios translate to planetary brightness temperatures of $926 \pm 12$ K (NRS1) and $1,122 \pm 10$ K (NRS2), which are low enough for a wide range of refractory condensates to form, including enstatite and forsterite. A nightside cloud deck may be blocking emission from deeper, hotter layers of the atmosphere, potentially helping to explain why cloud-free 3D general circulation model simulations systematically over-predict the nightside emission for WASP-121b. △ Less

Submitted 16 February, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: Accepted for publication in Astrophysical Journal Letters on December 29, 2022

Optical Properties of Organic Haze Analogues in Water-rich Exoplanet Atmospheres Observable with JWST

Authors: Chao He, Michael Radke, Sarah E. Moran, Sarah M. Horst, Nikole K. Lewis, Julianne I. Moses, Mark S. Marley, Natasha E. Batalha, Eliza M. -R. Kempton, Caroline V. Morley, Jeff A. Valenti, Veronique Vuitton

Abstract: JWST has begun its scientific mission, which includes the atmospheric characterization of transiting exoplanets. Some of the first exoplanets to be observed by JWST have equilibrium temperatures below 1000 K, which is a regime where photochemical hazes are expected to form. The optical properties of these hazes, which controls how they interact with light, are critical for interpreting exoplanet o… ▽ More JWST has begun its scientific mission, which includes the atmospheric characterization of transiting exoplanets. Some of the first exoplanets to be observed by JWST have equilibrium temperatures below 1000 K, which is a regime where photochemical hazes are expected to form. The optical properties of these hazes, which controls how they interact with light, are critical for interpreting exoplanet observations, but relevant experimental data are not available. Here we measure the density and optical properties of organic haze analogues generated in water-rich exoplanet atmosphere experiments. We report optical constants (0.4 to 28.6 μm) of organic haze analogues for current and future observational and modeling efforts covering the entire wavelength range of JWST instrumentation and a large part of Hubble. We use these optical constants to generate hazy model atmospheric spectra. The synthetic spectra show that differences in haze optical constants have a detectable effect on the spectra, impacting our interpretation of exoplanet observations. This study emphasizes the need to investigate the optical properties of hazes formed in different exoplanet atmospheres, and establishes a practical procedure to determine such properties. △ Less

Submitted 28 November, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 4 figures, 1 Table, Published in Nature Astronomy

Journal ref: Nature Astronomy (2023)

In Search of the Edge: A Bayesian Exploration of the Detectability of Red Edges in Exoplanet Reflection Spectra

Authors: Jonathan Gomez Barrientos, Ryan J. MacDonald, Nikole K. Lewis, Lisa Kaltenegger

Abstract: Reflection spectroscopy holds great promise for characterizing the atmospheres and surfaces of potentially habitable terrestrial exoplanets. The surface of the modern Earth exhibits a sharp albedo change near 750 nm caused by vegetation - the red edge - which would leave a strong spectral signature if present on an exoplanet. However, the retrieval of wavelength-dependent surface properties from r… ▽ More Reflection spectroscopy holds great promise for characterizing the atmospheres and surfaces of potentially habitable terrestrial exoplanets. The surface of the modern Earth exhibits a sharp albedo change near 750 nm caused by vegetation - the red edge - which would leave a strong spectral signature if present on an exoplanet. However, the retrieval of wavelength-dependent surface properties from reflection spectra has seen relatively little study. Here, we propose a new surface albedo parameterization capable of retrieving the wavelength location of a priori unknown 'edge-like' features. We demonstrate that a wavelength-dependent surface albedo model achieves higher accuracy in retrieving atmospheric composition. Wavelength-dependent surfaces are also generally preferred over a uniform albedo model when retrieving simulated reflection spectra for a modern Earth analog, even for moderate signal-to-noise ratios (S/N = 10) and Earth-like clouds. Further, the location of the modern Earth's red edge can be robustly and precisely constrained (within 70 nm for S/N = 10). Our results suggest that future space-based direct imaging missions have the potential to infer surface compositions for rocky exoplanets, including spectral edges similar to those caused by life on the modern Earth. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: 23 pages, 17 figures. Accepted for publication in ApJ

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

The Climate and Compositional Variation of the Highly Eccentric Planet HD 80606 b -- the rise and fall of carbon monoxide and elemental sulfur

Authors: Shang-Min Tsai, Maria Steinrueck, Vivien Parmentier, Nikole Lewis, Raymond Pierrehumbert

Abstract: The gas giant HD 80606 b has a highly eccentric orbit (e $\sim$ 0.93). The variation due to the rapid shift of stellar irradiation provides a unique opportunity to probe the physical and chemical timescales and to study the interplay between climate dynamics and atmospheric chemistry. In this work, we present integrated models to study the atmospheric responses and the underlying physical and chem… ▽ More The gas giant HD 80606 b has a highly eccentric orbit (e $\sim$ 0.93). The variation due to the rapid shift of stellar irradiation provides a unique opportunity to probe the physical and chemical timescales and to study the interplay between climate dynamics and atmospheric chemistry. In this work, we present integrated models to study the atmospheric responses and the underlying physical and chemical mechanisms of HD 80606 b. We first run three-dimensional general circulation models (GCMs) to establish the atmospheric thermal and dynamical structures for different atmospheric metallicities and internal heat. Based on the GCM output, we then adopted a 1D time-dependent photochemical model to investigate the compositional variation along the eccentric orbit. The transition of the circulation patterns of HD 80606 b matched the dynamics regimes in previous works. Our photochemical models show that efficient vertical mixing leads to deep quench levels of the major carbon and nitrogen species and the quenching behavior does not change throughout the eccentric orbit. Instead, photolysis is the main driver of the time-dependent chemistry. While CH$_4$ dominates over CO through most of the orbits, a transient state of [CO]/[CH$_4$}] $>$ 1 after periastron is confirmed for all metallicity and internal heat cases. The upcoming JWST Cycle 1 GO program will be able to track this real-time CH$_4$--CO conversion and infer the chemical timescale. Furthermore, sulfur species initiated by sudden heating and photochemical forcing exhibit both short-term and long-term cycles, opening an interesting avenue for detecting sulfur on exoplanets. △ Less

Submitted 16 March, 2023; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 18 pages, 18 figures, published in MNRAS

A Comparative L-dwarf Sample Exploring the Interplay Between Atmospheric Assumptions and Data Properties

Authors: Eileen C. Gonzales, Ben Burningham, Jacqueline K. Faherty, Nikole K. Lewis, Channon Visscher, Mark Marley

Abstract: Comparisons of atmospheric retrievals can reveal powerful insights on the strengths and limitations of our data and modeling tools. In this paper, we examine a sample of 5 similar effective temperature (Teff) or spectral type L dwarfs to compare their pressure-temperature (P-T) profiles. Additionally, we explore the impact of an object's metallicity and the observations' signal-to-noise (SNR) on t… ▽ More Comparisons of atmospheric retrievals can reveal powerful insights on the strengths and limitations of our data and modeling tools. In this paper, we examine a sample of 5 similar effective temperature (Teff) or spectral type L dwarfs to compare their pressure-temperature (P-T) profiles. Additionally, we explore the impact of an object's metallicity and the observations' signal-to-noise (SNR) on the parameters we can retrieve. We present the first atmospheric retrievals: 2MASS J15261405$+$2043414, 2MASS J05395200$-$0059019, 2MASS J15394189$-$0520428, and GD 165B increasing the small but growing number of L-dwarfs retrieved. When compared to atmospheric retrievals of SDSS J141624.08+134826.7, a low-metallicity d/sdL7 primary in a wide L+T binary, we find similar Teff sources have similar P-T profiles with metallicity differences impacting the relative offset between their P-T profiles in the photosphere. We also find that for near-infrared spectra, when the SNR is $\gtrsim80$ we are in a regime where model uncertainties dominate over data measurement uncertainties. As such, SNR does not play a role in the retrieval's ability to distinguish between a cloud-free and cloudless model, but may impact the confidence of the retrieved parameters. Lastly, we also discuss how to break cloud model degeneracies and the impact of extraneous gases in a retrieval model. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, 3 figures, 7 tables

ACCESS: Tentative detection of H$_2$O in the ground-based optical transmission spectrum of the low-density hot Saturn HATS-5b

Authors: Natalie H. Allen, Néstor Espinoza, Andrés Jordán, Mercedes López-Morales, Dániel Apai, Benjamin V. Rackham, James Kirk, David J. Osip, Ian C. Weaver, Chima McGruder, Kevin Ortiz Ceballos, Henrique Reggiani, Rafael Brahm, Florian Rodler, Nikole K Lewis, Jonathan Fraine

Abstract: We present a precise ground-based optical transmission spectrum of the hot-Saturn HATS-5b ($T_{eq} =1025$ K), obtained as part of the ACCESS survey with the IMACS multi-object spectrograph mounted on the Magellan/Baade Telescope. Our spectra cover the 0.5 to 0.9 micron region, and are the product of 5 individual transits observed between 2014 and 2018. We introduce the usage of additional second-o… ▽ More We present a precise ground-based optical transmission spectrum of the hot-Saturn HATS-5b ($T_{eq} =1025$ K), obtained as part of the ACCESS survey with the IMACS multi-object spectrograph mounted on the Magellan/Baade Telescope. Our spectra cover the 0.5 to 0.9 micron region, and are the product of 5 individual transits observed between 2014 and 2018. We introduce the usage of additional second-order light in our analyses which allows us to extract an "extra" transit light curve, improving the overall precision of our combined transit spectrum. We find that the favored atmospheric model for this transmission spectrum is a solar-metallicity atmosphere with sub-solar C/O, whose features are dominated by H$_2$O and with a depleted abundance of Na and K. If confirmed, this would point to a "clear" atmosphere at the pressure levels probed by transmission spectroscopy for HATS-5b. Our best-fit atmospheric model predicts a rich near-IR spectrum, which makes this exoplanet an excellent target for future follow-up observations with the James Webb Space Telescope, both to confirm this H$_2$O detection and to superbly constrain the atmosphere's parameters. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 31 pages, 16 figures, 7 tables. Accepted for publication in the Astronomical Journal

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

Bistability of the atmospheric circulation on TRAPPIST-1e

Authors: Denis E. Sergeev, Neil T. Lewis, F. Hugo Lambert, Nathan J. Mayne, Ian A. Boutle, James Manners, Krisztian Kohary

Abstract: Using a 3D general circulation model, we demonstrate that a confirmed rocky exoplanet and a primary observational target, TRAPPIST-1e presents an interesting case of climate bistability. We find that the atmospheric circulation on TRAPPIST-1e can exist in two distinct regimes for a 1~bar nitrogen-dominated atmosphere. One is characterized by a single strong equatorial prograde jet and a large day-… ▽ More Using a 3D general circulation model, we demonstrate that a confirmed rocky exoplanet and a primary observational target, TRAPPIST-1e presents an interesting case of climate bistability. We find that the atmospheric circulation on TRAPPIST-1e can exist in two distinct regimes for a 1~bar nitrogen-dominated atmosphere. One is characterized by a single strong equatorial prograde jet and a large day-night temperature difference; the other is characterized by a pair of mid-latitude prograde jets and a relatively small day-night contrast. The circulation regime appears to be highly sensitive to the model setup, including initial and surface boundary conditions, as well as physical parameterizations of convection and cloud radiative effects. We focus on the emergence of the atmospheric circulation during the early stages of simulations and show that the regime bistability is associated with a delicate balance between the zonally asymmetric heating, mean overturning circulation, and mid-latitude baroclinic instability. The relative strength of these processes places the GCM simulations on different branches of the evolution of atmospheric dynamics. The resulting steady states of the two regimes have consistent differences in the amount of water content and clouds, affecting the water absorption bands as well as the continuum level in the transmission spectrum, although they are too small to be detected with current technology. Nevertheless, this regime bistability affects the surface temperature, especially on the night side of the planet, and presents an interesting case for understanding atmospheric dynamics and highlights uncertainty in 3D GCM results, motivating more multi-model studies. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 31 pages, 14 figures, accepted to the Planetary Science Journal

Temperature Structures Associated with Different Components of the Atmospheric Circulation on Tidally Locked Exoplanets

Authors: Neil T. Lewis, Mark Hammond

Abstract: Observations of time-resolved thermal emission from tidally locked exoplanets can tell us about their atmospheric temperature structure. Telescopes such as JWST and ARIEL will improve the quality and availability of these measurements. This motivates an improved understanding of the processes that determine atmospheric temperature structure, particularly atmospheric circulation. The circulation is… ▽ More Observations of time-resolved thermal emission from tidally locked exoplanets can tell us about their atmospheric temperature structure. Telescopes such as JWST and ARIEL will improve the quality and availability of these measurements. This motivates an improved understanding of the processes that determine atmospheric temperature structure, particularly atmospheric circulation. The circulation is important in determining atmospheric temperatures, not only through its ability to transport heat, but also because any circulation pattern needs to be balanced by horizontal pressure contrasts, therefore implying a particular temperature structure. In this work, we show how the global temperature field on a tidally locked planet can be decomposed into contributions that are balanced by different components of the atmospheric circulation. These are the superrotating jet, stationary Rossby waves, and the divergent circulation. To achieve this, we partition the geopotential field into components balanced by the divergent circulation and the rotational circulation, with the latter comprising the jet and Rossby waves. The partitioned geopotential then implies a corresponding partitioning of the temperature via the hydrostatic relation. We apply these diagnostics to idealised general circulation model simulations, to show how the separate rotational and divergent circulations together make up the total three-dimensional atmospheric temperature structure. We also show how each component contributes distinct signatures to the thermal phase curve of a tidally locked planet. We conclude that this decomposition is a physically meaningful separation of the temperature field that explains its global structure, and can be used to fit observations of thermal emission. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: Submitted to AAS journals. Comments welcome and appreciated

A New Analysis of 8 Spitzer Phase Curves and Hot Jupiter Population Trends: Qatar-1b, Qatar-2b, WASP-52b, WASP-34b, and WASP-140b

Authors: Erin May, Kevin Stevenson, Jacob Bean, Taylor Bell, Nicolas Cowan, Lisa Dang, Jean-Michel Desert, Jonathan Fortney, Dylan Keating, Eliza Kempton, Thaddeus Komacek, Nikole Lewis, Megan Mansfield, Caroline Morley, Vivien Parmentier, Emily Rauscher, Mark Swain, Robert Zellem, Adam Showman

Abstract: With over 30 phase curves observed during the warm Spitzer mission, the complete data set provides a wealth of information relating to trends and three-dimensional properties of hot Jupiter atmospheres. In this work we present a comparative study of seven new Spitzer phase curves for four planets with equilibrium temperatures of T$_{eq}\sim$ 1300K: Qatar-2b, WASP-52b, WASP-34b, and WASP-140b, as w… ▽ More With over 30 phase curves observed during the warm Spitzer mission, the complete data set provides a wealth of information relating to trends and three-dimensional properties of hot Jupiter atmospheres. In this work we present a comparative study of seven new Spitzer phase curves for four planets with equilibrium temperatures of T$_{eq}\sim$ 1300K: Qatar-2b, WASP-52b, WASP-34b, and WASP-140b, as well as the reanalysis of the 4.5 $\micron$ Qatar-1b phase curve due to the similar equilibrium temperature. In total, five 4.5 $\micron$ phase curves and three 3.6 $\micron$ phase curves are analyzed here with a uniform approach. Using these new results, in combination with literature values for the entire population of published Spitzer phase curves of hot Jupiters, we present evidence for a linear trend of increasing hot spot offset with increasing orbital period, as well as observational evidence for two classes of planets in apparent redistribution vs. equilibrium temperature parameter space, and tentative evidence for a dependence of hot spot offset on planetary surface gravity in our $\sim$ 1300 K sample. We do not find trends in apparent heat redistribution with orbital period or gravity. Non-uniformity in literature Spitzer data analysis techniques precludes a definitive determination of the sources or lack of trends. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 22 pages, 13 figures, 6 tables. Accepted for publication in AAS journals

HST/WFC3 transmission spectroscopy of the cold rocky planet TRAPPIST-1h

Authors: L. J. Garcia, S. E. Moran, B. V. Rackham, H. R. Wakeford, M. Gillon, J. de Wit, N. K. Lewis

Abstract: TRAPPIST-1 is a nearby ultra-cool dwarf star transited by seven rocky planets. We observed three transits of its outermost planet, TRAPPIST-1h, using the G141 grism of the Wide Field Camera 3 instrument aboard the Hubble Space Telescope to place constraints on its potentially cold atmosphere. In order to deal with the effect of stellar contamination, we model TRAPPIST-1 active regions as portions… ▽ More TRAPPIST-1 is a nearby ultra-cool dwarf star transited by seven rocky planets. We observed three transits of its outermost planet, TRAPPIST-1h, using the G141 grism of the Wide Field Camera 3 instrument aboard the Hubble Space Telescope to place constraints on its potentially cold atmosphere. In order to deal with the effect of stellar contamination, we model TRAPPIST-1 active regions as portions of a cooler and a hotter photosphere, and generate multi-temperature models that we compare to the out-of-transit spectrum of the star. Using the inferred spot parameters, we produce corrected transmission spectra for planet h under five transit configurations and compare these data to planetary atmospheric transmission models using the forward model CHIMERA. Our analysis reveals that TRAPPIST-1h is unlikely to host an aerosol-free H/He-dominated atmosphere. While the current data precision limits the constraints we can put on the planetary atmosphere, we find that the likeliest scenario is that of a flat, featureless transmission spectrum in the WFC3/G141 bandpass due to a high mean molecular weight atmosphere (>1000x solar), no atmosphere, or an opaque aerosol layer, all in absence of stellar contamination. This work outlines the limitations of modeling active photospheric regions with theoretical stellar spectra, and those brought by our lack of knowledge of the photospheric structure of ultracool dwarf stars. Further characterization of the planetary atmosphere of TRAPPIST-1h would require higher precision measurements over wider wavelengths, which will be possible with the James Webb Space Telescope. △ Less

Submitted 30 March, 2022; v1 submitted 25 March, 2022; originally announced March 2022.

Journal ref: A&A 665, A19 (2022)

A Comprehensive Analysis of WASP-17b's Transmission Spectrum from Space-Based Observations

Authors: L. Alderson, H. R. Wakeford, R. J. MacDonald, N. K. Lewis, E. M. May, D. Grant, D. K. Sing, K. B. Stevenson, J. Fowler, J. Goyal, N. E. Batalha, T. Kataria

Abstract: Due to its 1770 K equilibrium temperature, WASP-17b, a 1.99 $R_\mathrm{Jup}$, 0.486 $M_\mathrm{Jup}$ exoplanet, sits at the critical juncture between hot and ultra-hot Jupiters. We present its 0.3-5 $μm$ transmission spectrum, with newly obtained with Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) measurements, and, taking advantage of improved analysis techniques, reanalysed HST Space Te… ▽ More Due to its 1770 K equilibrium temperature, WASP-17b, a 1.99 $R_\mathrm{Jup}$, 0.486 $M_\mathrm{Jup}$ exoplanet, sits at the critical juncture between hot and ultra-hot Jupiters. We present its 0.3-5 $μm$ transmission spectrum, with newly obtained with Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) measurements, and, taking advantage of improved analysis techniques, reanalysed HST Space Telescope Imaging Spectrograph (STIS) and Spitzer Space Telescope Infrared Array Camera (IRAC) observations. We achieve a median precision of 132 ppm, with a mean of 272 ppm across the whole spectrum. We additionally make use of Transiting Exoplanet Survey Satellite (TESS) and ground-based transit observations to refine the orbital period of WASP-17b. To interpret the observed atmosphere, we make use of free and equilibrium chemistry retrievals using the POSEIDON and ATMO retrieval codes respectively. We detect absorption due to H$_2$O at $>7 σ$, and find evidence of absorption due to CO$_2$ at $>3 σ$. We see no evidence of previously detected Na and K absorption. Across an extensive suite of retrieval configurations, we find the data favours a bimodal solution with high or low metallicity modes as a result of poor constraints in the optical, and demonstrate the importance of using multiple statistics for model selection. Future James Webb Space Telescope (JWST) GTO observations, combined with the presented transmission spectrum, will enable precise constraints on WASP-17b's atmosphere. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 25 pages, 13 figures, 10 tables. Accepted for publication in MNRAS

Diurnal variations in the stratosphere of the ultrahot giant exoplanet WASP-121b

Authors: Thomas Mikal-Evans, David K. Sing, Joanna K. Barstow, Tiffany Kataria, Jayesh Goyal, Nikole Lewis, Jake Taylor, Nathan. J. Mayne, Tansu Daylan, Hannah R. Wakeford, Mark S. Marley, Jessica J. Spake

Abstract: The temperature profile of a planetary atmosphere is a key diagnostic of radiative and dynamical processes governing the absorption, redistribution, and emission of energy. Observations have revealed dayside stratospheres that either cool or warm with altitude for a small number of gas giant exoplanets, while other dayside stratospheres are consistent with constant temperatures. Here we report spe… ▽ More The temperature profile of a planetary atmosphere is a key diagnostic of radiative and dynamical processes governing the absorption, redistribution, and emission of energy. Observations have revealed dayside stratospheres that either cool or warm with altitude for a small number of gas giant exoplanets, while other dayside stratospheres are consistent with constant temperatures. Here we report spectroscopic phase curve measurements for the gas giant WASP-121b, which constrain stratospheric temperatures throughout the diurnal cycle. Variations measured for a water vapour spectral feature reveal a temperature profile that transitions from warming with altitude on the dayside hemisphere to cooling with altitude on the nightside hemisphere. The data are well explained by models assuming chemical equilibrium, with water molecules thermally dissociating at low pressures on the dayside and recombining on the nightside. Nightside temperatures are low enough for perovskite (CaTiO3) to condense, which could deplete titanium from the gas phase and explain recent non-detections at the day-night terminator. Nightside temperatures are also consistent with the condensation of refractory species such as magnesium, iron, and vanadium. Detections of these metals at the day-night terminator suggest, however, that if they do form nightside clouds, cold trapping does not efficiently remove them from the upper atmosphere. Horizontal winds and vertical mixing could keep these refractory condensates aloft in the upper atmosphere of the nightside hemisphere until they are recirculated to the hotter dayside hemisphere and vaporised. △ Less

Submitted 20 February, 2022; originally announced February 2022.

Comments: Published in Nature Astronomy (publisher version is Open Access)

Why is it So Hot in Here? Exploring Population Trends in $\textit{Spitzer}$ Thermal Emission Observations of Hot Jupiters using Planet-Specific Self-Consistent Atmospheric Models

Authors: Jayesh M Goyal, Nikole K Lewis, Hannah R Wakeford, Ryan J MacDonald, Nathan J Mayne

Abstract: Thermal emission has now been observed from many dozens of exoplanet atmospheres, opening the gateway to population-level characterization. Here, we provide theoretical explanations for observed trends in $\textit{Spitzer}$ IRAC channel 1 (3.6 $μm$) and channel 2 (4.5 $μm$) photometric eclipse depths (EDs) across a population of 34 hot Jupiters. We apply planet-specific, self-consistent atmospheri… ▽ More Thermal emission has now been observed from many dozens of exoplanet atmospheres, opening the gateway to population-level characterization. Here, we provide theoretical explanations for observed trends in $\textit{Spitzer}$ IRAC channel 1 (3.6 $μm$) and channel 2 (4.5 $μm$) photometric eclipse depths (EDs) across a population of 34 hot Jupiters. We apply planet-specific, self-consistent atmospheric models, spanning a range of recirculation factors, metallicities, and C/O ratios, to probe the information content of $\textit{Spitzer}$ secondary eclipse observations across the hot-Jupiter population. We show that most hot Jupiters are inconsistent with blackbodies from $\textit{Spitzer}$ observations alone. We demonstrate that the majority of hot Jupiters are consistent with low energy redistribution between the dayside and nightside (hotter dayside than expected with efficient recirculation). We also see that high equilibrium temperature planets (T$_{eq}$ $\ge$ 1800 K) favor inefficient recirculation in comparison to the low temperature planets. Our planet-specific models do not reveal any definitive population trends in metallicity and C/O ratio with current data precision, but more than 59 % of our sample size is consistent with the C/O ratio $\leq$ 1 and 35 % are consistent with whole range (0.35 $\leq$ C/O $\leq$ 1.5). We also find that for most of the planets in our sample, 3.6 and 4.5 $μm$ model EDs lie within $\pm$1 $σ$ of the observed EDs. Intriguingly, few hot Jupiters exhibit greater thermal emission than predicted by the hottest atmospheric models (lowest recirculation) in our grid. Future spectroscopic observations of thermal emission from hot Jupiters with the James Webb Space Telescope will be necessary to robustly identify population trends in chemical compositions with its increased spectral resolution, range and data precision. △ Less

Submitted 6 January, 2022; originally announced January 2022.

Comments: 24 pages and 10 Figures. Published in the Astrophysical Journal (ApJ). Supplementary figures available in the online version of the journal

Journal ref: The Astrophysical Journal, Volume 923, Number 2, Page 242, 2021 December 20

Sensitivity of the Roman Coronagraph Instrument to Exozodiacal Dust

Authors: Ewan S Douglas, John Debes, Bertrand Mennesson, Bijan Nemati, Jaren Ashcraft, Bin Ren, Karl Stapelfeldt, Dmitry Savransky, Nikole K. Lewis, Bruce Macintosh

Abstract: Exozodiacal dust, warm debris from comets and asteroids in and near the habitable zone of stellar systems, reveals the physical processes that shape planetary systems. Scattered light from this dust is also a source of background flux which must be overcome by future missions to image Earthlike planets. This study quantifies the sensitivity of the Nancy Grace Roman Space Telescope Coronagraph to l… ▽ More Exozodiacal dust, warm debris from comets and asteroids in and near the habitable zone of stellar systems, reveals the physical processes that shape planetary systems. Scattered light from this dust is also a source of background flux which must be overcome by future missions to image Earthlike planets. This study quantifies the sensitivity of the Nancy Grace Roman Space Telescope Coronagraph to light scattered by exozodi, the zodiacal dust around other stars. Using a sample of 149 nearby stars, previously selected for optimum detection of habitable exoplanets by space observatories, we find the maximum number of exozodiacal disks with observable \textit{inner} habitable zone boundaries is six and the number of observable outer habitable boundaries is 74. One zodi was defined as the visible-light surface brightness of 22 $m_{\rm V}\ $arcsec$^{-2}$ around a solar-mass star, approximating the scattered light brightness in visible light at the Earth-equivalent insolation. In the speckle limited case, where the signal-to-noise ratio is limited by speckle temporal stability rather than shot noise, the median $5σ$ sensitivity to habitable zone exozodi is 12 zodi per resolution element. This estimate is calculated at the inner-working angle of the coronagraph, for the current best estimate performance, neglecting margins on the uncertainty in instrument performance and including a post-processing speckle suppression factor. For an log-norm distribution of exozodi levels with a median exozodi of 3$\times$ the solar zodi, we find that the Roman Coronagraph would be able to make 5$σ$ detections of exozodiacal disks in scattered light from 13 systems with a 95\% confidence interval spanning 7-20 systems. This sensitivity allows Roman Coronagraph to complement ground-based measurements of exozodiacal thermal emission and constrain dust albedos. △ Less

Submitted 23 December, 2021; originally announced December 2021.

Comments: 30 pages, 8 figures, 2 tables, accepted to PASP

TRIDENT: A Rapid 3D Radiative Transfer Model for Exoplanet Transmission Spectra

Authors: Ryan J. MacDonald, Nikole K. Lewis

Abstract: Transmission spectroscopy is one of the premier methods used to probe the temperature, composition, and cloud properties of exoplanet atmospheres. Recent studies have demonstrated that the multidimensional nature of exoplanet atmospheres -- due to non-uniformities across the day-night transition and between the morning and evening terminators -- can strongly influence transmission spectra. However… ▽ More Transmission spectroscopy is one of the premier methods used to probe the temperature, composition, and cloud properties of exoplanet atmospheres. Recent studies have demonstrated that the multidimensional nature of exoplanet atmospheres -- due to non-uniformities across the day-night transition and between the morning and evening terminators -- can strongly influence transmission spectra. However, the computational demands of 3D radiative transfer techniques have precluded their usage within atmospheric retrievals. Here we introduce TRIDENT, a new 3D radiative transfer model which rapidly computes transmission spectra of exoplanet atmospheres with day-night, morning-evening, and vertical variations in temperature, chemical abundances, and cloud properties. We also derive a general equation for transmission spectra, accounting for 3D atmospheres, refraction, multiple scattering, ingress/egress, grazing transits, stellar heterogeneities, and nightside thermal emission. After introducing TRIDENT's linear algebra-based approach to 3D radiative transfer, we propose new parametric prescriptions for 3D temperature and abundance profiles and 3D clouds. We show that multidimensional transmission spectra exhibit two significant observational signatures: (i) day-night composition gradients alter the relative amplitudes of absorption features; and (ii) morning-evening composition gradients distort the peak-to-wing contrast of absorption features. Finally, we demonstrate that these signatures of multidimensional atmospheres incur residuals > 100 ppm compared to 1D models, rendering them potentially detectable with JWST. TRIDENT's rapid radiative transfer, coupled with parametric multidimensional atmospheres, unlocks the final barrier to 3D atmospheric retrievals. △ Less

Submitted 31 December, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: 41 pages, 15 figures. Added appendix with TRIDENT runtimes. Accepted for publication in ApJ

Thermal Phase Curves of XO-3b: an Eccentric Hot Jupiter at the Deuterium Burning Limit

Authors: Lisa Dang, Taylor J. Bell, Nicolas B. Cowan, Daniel Thorngren, Tiffany Kataria, Heather A. Knutson, Nikole K. Lewis, Keivan G. Stassun, Jonathan J. Fortney, Eric Agol, Gregory P. Laughlin, Adam Burrows, Karen A. Collins, Drake Deming, Diana Jovmir, Jonathan Langton, Sara Rastegar, Adam P. Showman

Abstract: We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $μ$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $μ$m and $1610 \pm 70$ ppm at 4.5 $μ$m show no evidence of the previously reported dayside temperature inversion. We also empirically derive the mass and radius of XO-3b and its host star using Gaia DR3's parallax measureme… ▽ More We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $μ$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $μ$m and $1610 \pm 70$ ppm at 4.5 $μ$m show no evidence of the previously reported dayside temperature inversion. We also empirically derive the mass and radius of XO-3b and its host star using Gaia DR3's parallax measurement and find a planetary mass $M_p=11.79 \pm 0.98 ~M_{\rm{Jup}}$ and radius $R_p=1.295 \pm 0.066 ~R_{\rm{Jup}}$. We compare our \textit{Spitzer} observations with multiple atmospheric models to constrain the radiative and advective properties of XO-3b. While the decorrelated 4.5 $μ$m observations are pristine, the 3.6 $μ$m phase curve remains polluted with detector systematics due to larger amplitude intrapixel sensitivity variations in this channel. We focus our analysis on the more reliable 4.5 $μ$m phase curve and fit an energy balance model with solid body rotation to estimate the zonal wind speed and the pressure of the bottom of the mixed layer. Our energy balance model fit suggests an eastward equatorial wind speed of $3.13 ^{+0.26} _{-0.83}$ km/s, an atmospheric mixed layer down to $2.40 ^{+0.92} _{-0.16}$ bar, and Bond albedo of $0.106 ^{+0.008} _{-0.106}$. We assume that the wind speed and mixed layer depth are constant throughout the orbit. We compare our observations with a 1D planet-averaged model predictions at apoapse and periapse and 3D general circulation model (GCM) predictions for XO-3b. We also investigate the inflated radius of XO-3b and find that it would require an unusually large amount of internal heating to explain the observed planetary radius. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 23 pages, 14 figures, accepted for publication in AJ

Hiding in plain sight: observing planet-starspot crossings with the James Webb Space Telescope

Authors: Giovanni Bruno, Nikole K. Lewis, Jeff A. Valenti, Isabella Pagano, Tom J. Wilson, Everett Schlawin, Joshua Lothringer, Antonino F. Lanza, Jonathan Fraine, Gaetano Scandariato, Giuseppina Micela, Gianluca Cracchiolo

Abstract: Transiting exoplanets orbiting active stars frequently occult starspots and faculae on the visible stellar disc. Such occultations are often rejected from spectrophotometric transits, as it is assumed they do not contain relevant information for the study of exoplanet atmopsheres. However, they can provide useful constraints to retrieve the temperature of active features and their effect on transm… ▽ More Transiting exoplanets orbiting active stars frequently occult starspots and faculae on the visible stellar disc. Such occultations are often rejected from spectrophotometric transits, as it is assumed they do not contain relevant information for the study of exoplanet atmopsheres. However, they can provide useful constraints to retrieve the temperature of active features and their effect on transmission spectra. We analyse the capabilities of the James Webb Space Telescope in the determination of the spectra of occulted starspots, despite its lack of optical wavelength instruments on board. Focusing on K and M spectral types, we simulate starspots with different temperatures and in different locations of the stellar disc, and find that starspot temperatures can be determined to within a few hundred kelvins using NIRSpec/Prism and the proposed NIRCam/F150W2$+$F322W2's broad wavelength capabilities. Our results are particularly promising in the case of K and M dwarfs of mag$_K \leq 12.5$ with large temperature contrasts. △ Less

Submitted 18 November, 2021; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 16 pages, 13 figures, accepted for publication on MNRAS, updated references in Sections 1, 3.1 and 5, updated affiliation details, added link to uploaded material in Data Availability section

Gemini/GMOS Transmission Spectroscopy of the Grazing Planet Candidate WD 1856+534 b

Authors: Siyi Xu, Hannah Diamond-Lowe, Ryan J. MacDonald, Andrew Vanderburg, Simon Blouin, P. Dufour, Peter Gao, Laura Kreidberg, S. K. Leggett, Andrew W. Mann, Caroline V. Morley, Andrew W. Stephens, Christopher E. O'Connor, Pa Chia Thao, Nikole K. Lewis

Abstract: WD 1856+534 b is a Jupiter-sized, cool giant planet candidate transiting the white dwarf WD 1856+534. Here, we report an optical transmission spectrum of WD 1856+534 b obtained from ten transits using the Gemini Multi-Object Spectrograph. This system is challenging to observe due to the faintness of the host star and the short transit duration. Nevertheless, our phase-folded white light curve reac… ▽ More WD 1856+534 b is a Jupiter-sized, cool giant planet candidate transiting the white dwarf WD 1856+534. Here, we report an optical transmission spectrum of WD 1856+534 b obtained from ten transits using the Gemini Multi-Object Spectrograph. This system is challenging to observe due to the faintness of the host star and the short transit duration. Nevertheless, our phase-folded white light curve reached a precision of 0.12 %. WD 1856+534 b provides a unique transit configuration compared to other known exoplanets: the planet is $8\times$ larger than its star and occults over half of the stellar disc during mid-transit. Consequently, many standard modeling assumptions do not hold. We introduce the concept of a `limb darkening corrected, time-averaged transmission spectrum' and propose that this is more suitable than $(R_{\mathrm{p}, λ} / R_{\mathrm{s}})^2$ for comparisons to atmospheric models for planets with grazing transits. We also present a modified radiative transfer prescription. Though the transmission spectrum shows no prominent absorption features, it is sufficiently precise to constrain the mass of WD 1856+534 b to be > 0.84 M$_\mathrm{J}$ (to $2 \, σ$ confidence), assuming a clear atmosphere and a Jovian composition. High-altitude cloud decks can allow lower masses. WD 1856+534 b could have formed either as a result of common envelope evolution or migration under the Kozai-Lidov mechanism. Further studies of WD 1856+534 b, alongside new dedicated searches for substellar objects around white dwarfs, will shed further light on the mysteries of post-main sequence planetary systems. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 21 pages, 10 figures, accepted for publication in AJ

A comprehensive revisit of select Galileo/NIMS observations of Europa

Authors: Ishan Mishra, Nikole Lewis, Jonathan Lunine, Kevin P. Hand, Paul Helfenstein, R. W. Carlson, Ryan J. MacDonald

Abstract: The Galileo Near Infrared Mapping Spectrometer (NIMS) collected spectra of Europa in the 0.7-5.2 $μ$m wavelength region, which have been critical to improving our understanding of the surface composition of this moon. However, most of the work done to get constraints on abundances of species like water ice, hydrated sulfuric acid, hydrated salts and oxides have used proxy methods, such as absorpti… ▽ More The Galileo Near Infrared Mapping Spectrometer (NIMS) collected spectra of Europa in the 0.7-5.2 $μ$m wavelength region, which have been critical to improving our understanding of the surface composition of this moon. However, most of the work done to get constraints on abundances of species like water ice, hydrated sulfuric acid, hydrated salts and oxides have used proxy methods, such as absorption strength of spectral features or fitting a linear mixture of laboratory generated spectra. Such techniques neglect the effect of parameters degenerate with the abundances, such as the average grain-size of particles, or the porosity of the regolith. In this work we revisit three Galileo NIMS spectra, collected from observations of the trailing hemisphere of Europa, and use a Bayesian inference framework, with the Hapke reflectance model, to reassess Europa's surface composition. Our framework has several quantitative improvements relative to prior analyses: (1) simultaneous inclusion of amorphous and crystalline water ice, sulfuric-acid-octahydrate (SAO), CO$_2$, and SO$_2$; (2) physical parameters like regolith porosity and radiation-induced band-center shift; and (3) tools to quantify confidence in the presence of each species included in the model, constrain their parameters, and explore solution degeneracies. We find that SAO strongly dominates the composition in the spectra considered in this study, while both forms of water ice are detected at varying confidence levels. We find no evidence of either CO$_2$ or SO$_2$ in any of the spectra; we further show through a theoretical analysis that it is highly unlikely that these species are detectable in any 1-2.5 $μ$m Galileo NIMS data. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Comments: 30 pages, 12 figures. Published in the Planetary Science Journal

Journal ref: Planet. Sci. J., 2, 183 (2021)

The first retrieval of a substellar subdwarf: A cloud-free SDSS J125637.13-022452.4

Authors: Eileen C. Gonzales, Ben Burningham, Jacqueline K. Faherty, Channon Visscher, Mark Marley, Roxana Lupu, Richard Freedman, Nikole K. Lewis

Abstract: We present the first retrieval analysis of a substellar subdwarf, SDSS J125637.13-022452.4 (SDSS J1256-0224), using the Brewster retrieval code base. We find SDSS J1256-0224 is best fit by a cloud-free model with an ion (neutral H, H-, and electron) abundance corresponding to ion [Fe/H]=-1.5. However, this model is indistinguishable from a cloud-free model with ion [Fe/H]=-2.0 and a cloud-free mod… ▽ More We present the first retrieval analysis of a substellar subdwarf, SDSS J125637.13-022452.4 (SDSS J1256-0224), using the Brewster retrieval code base. We find SDSS J1256-0224 is best fit by a cloud-free model with an ion (neutral H, H-, and electron) abundance corresponding to ion [Fe/H]=-1.5. However, this model is indistinguishable from a cloud-free model with ion [Fe/H]=-2.0 and a cloud-free model with ion Fe/H]=-1.5 assuming a subsolar carbon-to-oxygen ratio. We are able to constrain abundances for water, FeH, and CrH, with an inability to constrain any carbon-bearing species likely due to the low-metallicity of SDSS J1256-0224. We also present an updated spectral energy distribution (SED) and semi-empirical fundamental parameters. Our retrieval- and SED-based fundamental parameters agree with the Baraffe low-metallicity evolutionary models. From examining our "rejected" models (those with $Δ$BIC>45), we find that we are able to retrieve gas abundances consistent with those of our best-fitting model. We find the cloud in these poorer fitting "cloudy" models is either pushed to the bottom of the atmosphere or made optically thin. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: 19 pages, 5 figures. arXiv admin note: text overlap with arXiv:2010.01224

On the Utility of Transmission Color Analysis I: Differentiating Super-Earths and Sub-Neptunes

Authors: Kristin S. Sotzen, Kevin B. Stevenson, Erin M. May, Natasha E. Batalha, Noam R. Izenberg, Sarah M. Horst, Calley L. Tinsman, Carey M. Lisse, Nikole K. Lewis, Jayesh M. Goyal, Joseph J. Linden, Kathleen E. Mandt

Abstract: The majority of exoplanets found to date have been discovered via the transit method, and transmission spectroscopy represents the primary method of studying these distant worlds. Currently, in-depth atmospheric characterization of transiting exoplanets entails the use of spectrographs on large telescopes, requiring significant observing time to study each planet. Previous studies have demonstrate… ▽ More The majority of exoplanets found to date have been discovered via the transit method, and transmission spectroscopy represents the primary method of studying these distant worlds. Currently, in-depth atmospheric characterization of transiting exoplanets entails the use of spectrographs on large telescopes, requiring significant observing time to study each planet. Previous studies have demonstrated trends for solar system worlds using color-color photometry of reflectance spectra, as well as trends within transmission spectra for hot Jupiters. Building on these concepts, we have investigated the use of transmission color photometric analysis for efficient, coarse categorization of exoplanets and for assessing the nature of these worlds, with a focus on resolving the bulk composition degeneracy to aid in discriminating super-Earths and sub-Neptunes. We present our methodology and first results, including spectrum models, model comparison frameworks, and wave band selection criteria. We present our results for different transmission "color" metrics, filter selection methods, and numbers of filters. Assuming noise-free spectra of isothermal atmospheres in chemical equilibrium, with our pipeline, we are able to constrain atmospheric mean molecular weight in order to distinguish between super-Earth and sub-Neptune atmospheres with >90$\%$ overall accuracy using as few as two specific low-resolution filter combinations. We also found that increasing the number of filters does not substantially impact this performance. This method could allow for broad characterization of large numbers of planets much more efficiently than current methods permit, enabling population and system-level studies. Additionally, data collected via this method could inform follow-up observing time by large telescopes for more detailed studies of worlds of interest. △ Less

Submitted 6 September, 2021; originally announced September 2021.

The Pandora SmallSat: Multiwavelength Characterization of Exoplanets and their Host Stars

Authors: Elisa V. Quintana, Knicole D. Colón, Gregory Mosby, Joshua E. Schlieder, Pete Supsinskas, Jordan Karburn, Jessie L. Dotson, Thomas P. Greene, Christina Hedges, Dániel Apai, Thomas Barclay, Jessie L. Christiansen, Néstor Espinoza, Susan E. Mullally, Emily A. Gilbert, Kelsey Hoffman, Veselin B. Kostov, Nikole K. Lewis, Trevor O. Foote, James Mason, Allison Youngblood, Brett M. Morris, Elisabeth R. Newton, Joshua Pepper, Benjamin V. Rackham , et al. (2 additional authors not shown)

Abstract: Pandora is a SmallSat mission designed to study the atmospheres of exoplanets, and was selected as part of NASA's Astrophysics Pioneers Program. Transmission spectroscopy of transiting exoplanets provides our best opportunity to identify the makeup of planetary atmospheres in the coming decade. Stellar brightness variations due to star spots, however, can impact these measurements and contaminate… ▽ More Pandora is a SmallSat mission designed to study the atmospheres of exoplanets, and was selected as part of NASA's Astrophysics Pioneers Program. Transmission spectroscopy of transiting exoplanets provides our best opportunity to identify the makeup of planetary atmospheres in the coming decade. Stellar brightness variations due to star spots, however, can impact these measurements and contaminate the observed spectra. Pandora's goal is to disentangle star and planet signals in transmission spectra to reliably determine exoplanet atmosphere compositions. Pandora will collect long-duration photometric observations with a visible-light channel and simultaneous spectra with a near-IR channel. The broad-wavelength coverage will provide constraints on the spot and faculae covering fractions of low-mass exoplanet host stars and the impact of these active regions on exoplanetary transmission spectra. Pandora will subsequently identify exoplanets with hydrogen- or water-dominated atmospheres, and robustly determine which planets are covered by clouds and hazes. Pandora will observe at least 20 exoplanets with sizes ranging from Earth-size to Jupiter-size and host stars spanning mid-K to late-M spectral types. The project is made possible by leveraging investments in other projects, including an all-aluminum 0.45-meter Cassegrain telescope design, and a NIR sensor chip assembly from the James Webb Space Telescope. The mission will last five years from initial formulation to closeout, with one-year of science operations. Launch is planned for the mid-2020s as a secondary payload in Sun-synchronous low-Earth orbit. By design, Pandora has a diverse team, with over half of the mission leadership roles filled by early career scientists and engineers, demonstrating the high value of SmallSats for developing the next generation of space mission leaders. △ Less

Submitted 19 August, 2021; v1 submitted 13 August, 2021; originally announced August 2021.

Comments: Proceedings of the Small Satellite Conference, Science/Mission Payloads, SSC21-VI-02 (2021)

The Emission Spectrum of the Hot Jupiter WASP-79b from HST/WFC3

Authors: Trevor O. Foote, Nikole K. Lewis, Brian M. Kilpatrick, Jayesh M. Goyal, Giovanni Bruno, Hannah R. Wakeford, Nina Robbins-Blanch, Tiffany Kataria, Ryan J. MacDonald, Mercedes López-Morales, David K. Sing, Thomas Mikal-Evans, Vincent Bourrier, Gregory Henry, Lars A. Buchhave

Abstract: Here we present a thermal emission spectrum of WASP-79b, obtained via Hubble Space Telescope Wide Field Camera 3 G141 observations as part of the PanCET program. As we did not observe the ingress or egress of WASP-79b's secondary eclipse, we consider two scenarios: a fixed mid-eclipse time based on the expected occurrence time, and a mid-eclipse time as a free parameter. In both scenarios, we can… ▽ More Here we present a thermal emission spectrum of WASP-79b, obtained via Hubble Space Telescope Wide Field Camera 3 G141 observations as part of the PanCET program. As we did not observe the ingress or egress of WASP-79b's secondary eclipse, we consider two scenarios: a fixed mid-eclipse time based on the expected occurrence time, and a mid-eclipse time as a free parameter. In both scenarios, we can measure thermal emission from WASP-79b from 1.1 to 1.7 $μ$m at 2.4$σ$ confidence consistent with a 1900 K brightness temperature for the planet. We combine our observations with Spitzer dayside photometry (3.6 and 4.5 $μ$m) and compare these observations to a grid of atmospheric forward models that span a range of metallicities, carbon-to-oxygen ratios, and recirculation factors. Given the strength of the planetary emission and the precision of our measurements, we found a wide range of forward models to be consistent with our data. The best-match equilibrium model suggests that WASP-79b's dayside has a solar metallicity and carbon-to-oxygen ratio, alongside a recirculation factor of 0.75. Models including significant H- opacity provide the best match to WASP-79b's emission spectrum near 1.58 $μ$m. However, models featuring high-temperature cloud species-formed via vigorous vertical mixing and low sedimentation efficiencies-with little day-to-night energy transport also match WASP-79b's emission spectrum. Given the broad range of equilibrium chemistry, disequilibrium chemistry, and cloudy atmospheric models consistent with our observations of WASP-79b's dayside emission, further observations will be necessary to constrain WASP-79b's dayside atmospheric properties. △ Less

Submitted 14 December, 2021; v1 submitted 29 July, 2021; originally announced July 2021.

Comments: 13 pages, 8 figures. Accepted for publication in AJ

Haze Evolution in Temperate Exoplanet Atmospheres Through Surface Energy Measurements

Authors: Xinting Yu, Chao He, Xi Zhang, Sarah M. Hörst, Austin H. Dymont, Patricia McGuiggan, Julianne I. Moses, Nikole K. Lewis, Jonathan J. Fortney, Peter Gao, Eliza M. -R. Kempton, Sarah E. Moran, Caroline V. Morley, Diana Powell, Jeff A. Valenti, Véronique Vuitton

Abstract: Photochemical hazes are important opacity sources in temperate exoplanet atmospheres, hindering current observations from characterizing exoplanet atmospheric compositions. The haziness of an atmosphere is determined by the balance between haze production and removal. However, the material-dependent removal physics of the haze particles is currently unknown under exoplanetary conditions. Here we p… ▽ More Photochemical hazes are important opacity sources in temperate exoplanet atmospheres, hindering current observations from characterizing exoplanet atmospheric compositions. The haziness of an atmosphere is determined by the balance between haze production and removal. However, the material-dependent removal physics of the haze particles is currently unknown under exoplanetary conditions. Here we provide experimentally-measured surface energies for a grid of temperate exoplanet hazes to characterize haze removal in exoplanetary atmospheres. We found large variations of surface energies for hazes produced under different energy sources, atmospheric compositions, and temperatures. The surface energies of the hazes were found to be the lowest around 400 K for the cold plasma samples, leading to the lowest removal rates. We show a suggestive correlation between haze surface energy and atmospheric haziness with planetary equilibrium temperature. We hypothesize that habitable zone exoplanets could be less hazy, as they would possess high-surface-energy hazes which can be removed efficiently. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 6 figures, 2 tables, 23 pages

Spitzer phase curve observations and circulation models of the inflated ultra-hot Jupiter WASP-76b

Authors: Erin M. May, Thaddeus D. Komacek, Kevin B. Stevenson, Eliza M. -R. Kempton, Jacob L. Bean, Matej Malik, Jegug Ih, Megan Mansfield, Arjun B. Savel, Drake Deming, Jean-Michel Desert, Y. Katherina Feng, Jonathan J. Fortney, Tiffany Kataria, Nikole Lewis, Caroline Morley, Emily Rauscher, Adam Showman

Abstract: The large radii of many hot Jupiters can only be matched by models that have hot interior adiabats, and recent theoretical work has shown that the interior evolution of hot Jupiters has a significant impact on their atmospheric structure. Due to its inflated radius, low gravity, and ultra-hot equilibrium temperature, WASP-76b is an ideal case study for the impact of internal evolution on observabl… ▽ More The large radii of many hot Jupiters can only be matched by models that have hot interior adiabats, and recent theoretical work has shown that the interior evolution of hot Jupiters has a significant impact on their atmospheric structure. Due to its inflated radius, low gravity, and ultra-hot equilibrium temperature, WASP-76b is an ideal case study for the impact of internal evolution on observable properties. Hot interiors should most strongly affect the non-irradiated side of the planet, and thus full phase curve observations are critical to ascertain the effect of the interior on the atmospheres of hot Jupiters. In this work, we present the first Spitzer phase curve observations of WASP-76b. We find that WASP-76b has an ultra-hot day side and relatively cold nightside with brightness temperatures of $2471 \pm 27~\mathrm{K}$/$1518 \pm 61~\mathrm{K}$ at $3.6~\micron$ and $2699 \pm 32~\mathrm{K}$/$1259 \pm 44~\mathrm{K}$ at $4.5~\micron$, respectively. These results provide evidence for a dayside thermal inversion. Both channels exhibit small phase offsets of $0.68 \pm 0.48^{\circ}$ at $3.6~\micron$ and $0.67 \pm 0.2^{\circ}$ at $4.5~μ\mathrm{m}$. We compare our observations to a suite of general circulation models that consider two end-members of interior temperature along with a broad range of frictional drag strengths. Strong frictional drag is necessary to match the small phase offsets and cold nightside temperatures observed. From our suite of cloud-free GCMs, we find that only cases with a cold interior can reproduce the cold nightsides and large phase curve amplitude at $4.5~\micron$, hinting that the hot interior adiabat of WASP-76b does not significantly impact its atmospheric dynamics or that clouds blanket its nightside. △ Less

Submitted 18 August, 2021; v1 submitted 7 July, 2021; originally announced July 2021.

Comments: 24 pages, 10 Figures, 5 Tables. Accepted to AJ. Co-First Authors

L 98-59: a Benchmark System of Small Planets for Future Atmospheric Characterization

Authors: Daria Pidhorodetska, Sarah E. Moran, Edward W. Schwieterman, Thomas Barclay, Thomas J. Fauchez, Nikole K. Lewis, Elisa V. Quintana, Geronimo L. Villanueva, Shawn D. Domagal-Goldman, Joshua E. Schlieder, Emily A. Gilbert, Stephen R. Kane, Veselin B. Kostov

Abstract: L 98-59 is an M3V dwarf star that hosts three small (R < 1.6 Earth radii) planets. The host star is bright (K = 7.1) and nearby (10.6 pc), making the system a prime target for follow-up characterization with the Hubble Space Telescope (HST) and the upcoming James Webb Space Telescope (JWST). Herein, we use simulated transmission spectroscopy to evaluate the detectability of spectral features with… ▽ More L 98-59 is an M3V dwarf star that hosts three small (R < 1.6 Earth radii) planets. The host star is bright (K = 7.1) and nearby (10.6 pc), making the system a prime target for follow-up characterization with the Hubble Space Telescope (HST) and the upcoming James Webb Space Telescope (JWST). Herein, we use simulated transmission spectroscopy to evaluate the detectability of spectral features with HST and JWST assuming diverse atmospheric scenarios (e.g., atmospheres dominated by H2, H2O, CO2, or O2). We find that H2O and CH4 present in a low mean-molecular weight atmosphere could be detected with HST in 1 transit for the two outermost planets, while H2O in a clear steam atmosphere could be detected in 6 transits or fewer with HST for all three planets. We predict that observations using JWST/NIRISS would be capable of detecting a clear steam atmosphere in 1 transit for each planet, and H2O absorption in a hazy steam atmosphere in 2 transits or less. In a clear, desiccated atmosphere, O2 absorption may be detectable for all three planets with NIRISS. If the L 98-59 planets possess a clear, Venus-like atmosphere, NIRSpec could detect CO2 within 26 transits for each planet, but the presence of H2SO4 clouds would significantly suppress CO2 absorption. The L 98-59 system is an excellent laboratory for comparative planetary studies of transiting multiplanet systems, and observations of the system via HST and JWST would present a unique opportunity to test the accuracy of the models presented in this study. △ Less

Submitted 1 June, 2021; originally announced June 2021.

Comments: Submitted to AJ

ACCESS & LRG-BEASTS: a precise new optical transmission spectrum of the ultrahot Jupiter WASP-103b

Authors: James Kirk, Ben Rackham, Ryan MacDonald, Mercedes López-Morales, Néstor Espinoza, Monika Lendl, Jamie Wilson, David J. Osip, Peter J. Wheatley, Ian Skillen, Dániel Apai, Alex Bixel, Neale P. Gibson, Andrés Jordan, Nikole K. Lewis, Tom Louden, Chima D. McGruder, Nikolay Nikolov, Florian Rodler, Ian C. Weaver

Abstract: We present a new ground-based optical transmission spectrum of the ultrahot Jupiter WASP-103b ($T_{eq} = 2484$K). Our transmission spectrum is the result of combining five new transits from the ACCESS survey and two new transits from the LRG-BEASTS survey with a reanalysis of three archival Gemini/GMOS transits and one VLT/FORS2 transit. Our combined 11-transit transmission spectrum covers a wavel… ▽ More We present a new ground-based optical transmission spectrum of the ultrahot Jupiter WASP-103b ($T_{eq} = 2484$K). Our transmission spectrum is the result of combining five new transits from the ACCESS survey and two new transits from the LRG-BEASTS survey with a reanalysis of three archival Gemini/GMOS transits and one VLT/FORS2 transit. Our combined 11-transit transmission spectrum covers a wavelength range of 3900--9450A with a median uncertainty in the transit depth of 148 parts-per-million, which is less than one atmospheric scale height of the planet. In our retrieval analysis of WASP-103b's combined optical and infrared transmission spectrum, we find strong evidence for unocculted bright regions ($4.3σ$) and weak evidence for H$_2$O ($1.9σ$), HCN ($1.7σ$), and TiO ($2.1σ$), which could be responsible for WASP-103b's observed temperature inversion. Our optical transmission spectrum shows significant structure that is in excellent agreement with the extensively studied ultrahot Jupiter WASP-121b, for which the presence of VO has been inferred. For WASP-103b, we find that VO can only provide a reasonable fit to the data if its abundance is implausibly high and we do not account for stellar activity. Our results highlight the precision that can be achieved by ground-based observations and the impacts that stellar activity from F-type stars can have on the interpretation of exoplanet transmission spectra. △ Less

Submitted 30 April, 2021; originally announced May 2021.

Comments: 33 pages, 17 figures, 7 tables. Accepted for publication in AJ

HST PanCET Program: A Complete Near-UV to Infrared Transmission Spectrum for the Hot Jupiter WASP-79b

Authors: Alexander D. Rathcke, Ryan J. MacDonald, Joanna K. Barstow, Jayesh M. Goyal, Mercedes Lopez-Morales, João M. Mendonça, Jorge Sanz-Forcada, Gregory W. Henry, David K. Sing, Munazza K. Alam, Nikole K. Lewis, Katy L. Chubb, Jake Taylor, Nikolay Nikolov, Lars A. Buchhave

Abstract: We present a new optical transmission spectrum of the hot Jupiter WASP-79b. We observed three transits with the STIS instrument mounted on HST, spanning 0.3 - 1.0 um. Combining these transits with previous observations, we construct a complete 0.3 - 5.0 um transmission spectrum of WASP-79b. Both HST and ground-based observations show decreasing transit depths towards blue wavelengths, contrary to… ▽ More We present a new optical transmission spectrum of the hot Jupiter WASP-79b. We observed three transits with the STIS instrument mounted on HST, spanning 0.3 - 1.0 um. Combining these transits with previous observations, we construct a complete 0.3 - 5.0 um transmission spectrum of WASP-79b. Both HST and ground-based observations show decreasing transit depths towards blue wavelengths, contrary to expectations from Rayleigh scattering or hazes. We infer atmospheric and stellar properties from the full near-UV to infrared transmission spectrum of WASP-79b using three independent retrieval codes, all of which yield consistent results. Our retrievals confirm previous detections of H$_{2}$O (at 4.0$σ$ confidence), while providing moderate evidence of H$^{-}$ bound-free opacity (3.3$σ$) and strong evidence of stellar contamination from unocculted faculae (4.7$σ$). The retrieved H$_{2}$O abundance ($\sim$ 1$\%$) suggests a super-stellar atmospheric metallicity, though stellar or sub-stellar abundances remain consistent with present observations (O/H = 0.3 - 34$\times$ stellar). All three retrieval codes obtain a precise H$^{-}$ abundance constraint: log(X$_{\rm{H^{-}}}$) $\approx$ -8.0 $\pm$ 0.7. The potential presence of H$^{-}$ suggests that JWST observations may be sensitive to ionic chemistry in the atmosphere of WASP-79b. The inferred faculae are $\sim$ 500 K hotter than the stellar photosphere, covering $\sim$ 15$\%$ of the stellar surface. Our analysis underscores the importance of observing UV - optical transmission spectra in order to disentangle the influence of unocculted stellar heterogeneities from planetary transmission spectra. △ Less

Submitted 21 April, 2021; originally announced April 2021.

Comments: 27 pages, 10 figures. Resubmitted to AJ after referee report

ACCESS: An optical transmission spectrum of the high-gravity, hot Jupiter HAT-P-23b

Authors: Ian C. Weaver, Mercedes López-Morales, Munazza K. Alam, Néstor Espinoza, Benjamin V. Rackham, Jayesh M. Goyal, Ryan J. MacDonald, Nikole K. Lewis, Dániel Apai, Alex Bixel, Andrés Jordán, James Kirk, Chima McGruder, David J. Osip

Abstract: We present a new ground-based visible transmission spectrum of the high-gravity, hot Jupiter HAT-P-23b, obtained as part of the ACCESS project. We derive the spectrum from five transits observed between 2016 and 2018, with combined wavelength coverage between 5200 Å - 9269 Å in 200 Å bins, and with a median precision of 247 ppm per bin. HAT-P-23b's relatively high surface gravity (g ~ 30 m/s^2), c… ▽ More We present a new ground-based visible transmission spectrum of the high-gravity, hot Jupiter HAT-P-23b, obtained as part of the ACCESS project. We derive the spectrum from five transits observed between 2016 and 2018, with combined wavelength coverage between 5200 Å - 9269 Å in 200 Å bins, and with a median precision of 247 ppm per bin. HAT-P-23b's relatively high surface gravity (g ~ 30 m/s^2), combined with updated stellar and planetary parameters from Gaia DR2, gives a 5-scale-height signal of 384 ppm for a hydrogen-dominated atmosphere. Bayesian models favor a clear atmosphere for the planet with the tentative presence of TiO, after simultaneously modeling stellar contamination, using spots parameter constraints from photometry. If confirmed, HAT-P-23b would be the first example of a high-gravity gas giant with a clear atmosphere observed in transmission at optical/NIR wavelengths; therefore, we recommend expanding observations to the UV and IR to confirm our results and further characterize this planet. This result demonstrates how combining transmission spectroscopy of exoplanet atmospheres with long-term photometric monitoring of the host stars can help disentangle the exoplanet and stellar activity signals. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 28 pages, 18 Figures, accepted for publication in AJ. arXiv admin note: text overlap with arXiv:1911.03358